4.1. Particle mass flow rates

The particle mass flow rate for particles exiting a hopper can be estimated from the well-established Beverloo equation (Beverloo, Leniger & Van de Velde Reference Beverloo, Leniger and Van de Velde1961), given by

(4.1)\begin{equation} \dot M_p=C_1 \rho_b g^{0.5} D^{1.5}_{eff}, \end{equation}

where $D_{eff}=D-C_2 \skew6\bar{d}_p$ is the effective hopper outlet thickness, which takes into account the influence of a thin layer of slow (or non-) moving particles near to the wall of the hopper exit, resulting in a reduction in the effective hopper outlet size (Beverloo et al. Reference Beverloo, Leniger and Van de Velde1961; Mankoc et al. Reference Mankoc, Janda, Arevalo, Pastor, Zuriguel, Garcimartín and Maza2007; Janda et al. Reference Janda, Zuriguel and Maza2012). The degree to which the effective hopper outlet size is reduced by the presence of particles is taken into account through the coefficient $C_2$, which is likely influenced by the particle and hopper surface properties (Beverloo et al. Reference Beverloo, Leniger and Van de Velde1961; Mankoc et al. Reference Mankoc, Janda, Arevalo, Pastor, Zuriguel, Garcimartín and Maza2007). The coefficient $C_1$ in (4.1) is a lumped parameter which partially takes into account the flow resistance due to particle–particle and particle–wall friction, with the value of $C_1$ expected to be approximately inversely proportional to the friction within the flow.

Re-arranging (4.1) results in the normalised particle mass flow rate per metre of hopper width as follows:

(4.2)\begin{equation} \tilde{m}=\frac{\dot M_p}{\rho_b g^{0.5} D^{1.5}}= C_1 \left(1-\frac{C_2 \skew6\bar{d}_p}{D} \right)^{1.5}. \end{equation}

Since $C_1$ may be influenced by particle size, (4.2) can be fitted for any given particle size with the measured experimental data to obtain $C_1$ as a function of the reciprocal of particle diameter, $1/\skew6\bar{d}_p$. The results are presented in figure 5(a). These show that the value of $C_1$ is approximately linearly proportional to the reciprocal of particle diameter, increasing with $\skew6\bar{d}_p$. Since $\tilde {m}$ is linearly proportional to $C_1$, this implies that, all else being equal, the particle mass flow rate will increase with $\skew6\bar{d}_p$. This is attributed to the decrease in surface area (per unit volume) as $\skew6\bar{d}_p$ is increased, which in turn reduces particle–particle friction. A linear curve fit of the data in figure 5(a) results in $C_1 \approx -0.0416/\skew6\bar{d}_p +1.91$ (note that the first coefficient on the right-hand side of the equation has dimensions of mm), which suggests that, for the cases where the particles are large, $C_1 \rightarrow 1.9$. Nevertheless, the values of $1.63 \leq C_1 \leq 1.84$ found here are significantly different from the $C_1 = 1.033$ obtained by Kim et al. (Reference Kim, Siegel, Kolb, Rangaswamy and Moujaes2009) and Ho et al. (Reference Ho, Christian, Romano, Yellowhair, Siegel, Savoldi and Zanino2017) using similar particles, and from the value of $C_1 = 0.583$ measured by Beverloo et al. (Reference Beverloo, Leniger and Van de Velde1961) using agricultural solids. The higher value of $C_1$ found in the present experiment is attributed to differences in hopper angle. For example, the present experiment utilises wedged-shaped hoppers with a wall inclination half-angle of $5^{\circ }$, while the hoppers used in Beverloo et al. (Reference Beverloo, Leniger and Van de Velde1961), Kim et al. (Reference Kim, Siegel, Kolb, Rangaswamy and Moujaes2009) and Ho et al. (Reference Ho, Christian, Romano, Yellowhair, Siegel, Savoldi and Zanino2017) had the opening on a flat bottom, with a hopper wall inclination half-angle of $90^{\circ }$. It is known from the previous studies that a smaller hopper angle leads to an increase in the particle mass flow rates (Brown Reference Brown1961; Anand et al. Reference Anand, Curtis, Wassgren, Hancock and Ketterhagen2008). Taking the effect of hopper angle into account, a greater value of $C_1$ found in the present measurements is therefore consistent with previous work.

Figure 5. The variation of the parameter $C_1$ as a function of the inverse of particle diameter, together with the measured mass flow rate parameter $(\tilde {m}/C_1 )^{2/3}$ as a function of normalised particle size. Here, $C_1$ and $\tilde {m}$ are from (4.2). The black dashed line in both figures represents a linear curve fit of the data.

Figure 5(b) presents the term $(\tilde {m}/C_1 )^{2/3}$ as a function of the normalised particle diameter, $\skew6\bar{d}_p/D$, for all investigated hopper sizes and particle diameters. The axes in this figure were chosen on the basis that the data should follow a linear profile according to (4.2). The results show that the data for all particle diameters and hopper sizes collapse onto a single linear profile. A linear curve fit of the data results in a profile following the formula $(\tilde {m}/C_1 )^{2/3}=1-1.64 \skew6\bar{d}_p/D$. That is, the value of $C_2$ found in the present experiments is approximately constant at $C_2=1.64$. The constant $C_2$ is slightly higher than the value 1.4 reported by Beverloo et al. (Reference Beverloo, Leniger and Van de Velde1961), Kim et al. (Reference Kim, Siegel, Kolb, Rangaswamy and Moujaes2009) and Ho et al. (Reference Ho, Christian, Romano, Yellowhair, Siegel, Savoldi and Zanino2017). The reason for this small discrepancy is unclear, although it could be due to differences in the particle properties (e.g. the particle–particle friction) between the present experiments and the previous studies. Nevertheless, the results show that with a constant value of $C_2$, and with $C_1$ varying with $\skew6\bar{d}_p$ according to figure 5(a), the measured mass flow rates are close to that estimated from the linear curve fit formula to within 5.8 %. This result provides evidence that $C_1$ is highly dependent on particle size, while $C_2$ is relatively insensitive to particle size, for any given particle type.

The results presented here also demonstrate that the Beverloo equation particularly suits the scaling assessments of flows in a hopper because it is an analytical equation that directly accounts for the dominant dimensionless parameters, whilst lumping the second-order terms to incorporate, albeit imperfectly, the influence of highly complex and nonlinear processes, into the two empirical parameters $C_1$ and $C_2$. These parameters simplify with reasonable accuracy the complex particle dynamics in the granular (highly dense) flow regime inside the hopper and include the influences of particle size distribution, surface roughness and shape on the particle–particle and particle–wall interactions. The challenges of advancing the understanding in this environment mean that, despite many attempts to improve on it for specific scenarios, no widely acknowledged method has yet emerged.

4.2. Particle curtain expansion

Figure 6 presents a representative instantaneous image of the side view of the free-falling particle curtain. As can be seen, the particle curtain thickness is small relative to its streamwise length, with the particles falling predominantly in the direction of gravity, consistent with previous work (Ho et al. Reference Ho, Christian, Romano, Yellowhair, Siegel, Savoldi and Zanino2017). Visual inspection of the particle curtain flow also reveals that there are three distinct flow regions which characterise the particle curtain flow in the direction of particle fall, namely a near-field expansion, a neck zone and an intermediate-field expansion. The ‘near-field expansion’ occurs near to the hopper exit (i.e. approximately within the first few hopper outlet thicknesses from the exit), and is characterised by strong particle dispersion laterally from the centreline, as can be seen from the zoomed-in inset on the right-hand side of figure 6. Downstream from the near-field expansion region, the particle curtain appears to reduce in thickness with streamwise distance. This region is called the ‘neck zone’, and may extend as far downstream as $x/D \approx 100$. Further downstream from this neck zone, the curtain can be seen to expand again with downstream distance. In this ‘intermediate-field expansion’, the curtain can be seen to exhibit large-scale flow unsteadiness that spans the entire width of the curtain. This unsteadiness is consistent with previous work, being driven by the shear force between the curtain and the surrounding flow, which generates such roll-up features (Möbius Reference Möbius2006; Amarouchene, Boudet & Kellay Reference Amarouchene, Boudet and Kellay2008; Lau & Nathan Reference Lau and Nathan2016; Chu et al. Reference Chu, Wang, Zheng, Yu and Pan2020). It should be noted that the term ‘intermediate field’ is used here, as distinct from the more commonly used ‘far field’, because the latter typically refers to the region of the flow where particles approach their settling velocity (Sedaghatizadeh et al. Reference Sedaghatizadeh, Arjomandi, Lau and Nathan2022), which is not necessarily true in the present intermediate field.

Figure 6. An instantaneous side-view image of the particle curtain, with an illustration of the three characteristic regions of the falling particle curtain, namely the near field, the neck zone and the intermediate field. In this case, $D = 5$ mm and $\skew6\bar{d}_p = 163\,\mathrm {\mu }$m. Also shown is the dilute-particle region where particles at the edge of the curtain interact with the surrounding flow, resulting in a region of low particle volume fraction relative to the core of the curtain.

Visualisation of the particle curtain flow also shows that the lateral distribution of particles can be characterised by a local ‘core’ region of densely seeded particles close to the centreline where the flow opacity approaches unity (the black regions in the figure), whose position fluctuates with time in the downstream region of the curtain. Surrounding this ‘core’ region, a dilute-particle region can be seen where the opacity is noticeably lower than unity. From a practical perspective, the lateral dispersion of particles away from the main curtain is also important, since it is likely to be a key contributor to the generation of ‘dust’ that has potential to escape from devices such as FPRs to generate both a challenge to reduce pollution and a loss of efficiency.

Figure 7 presents instantaneous side-view images of the particle curtain issuing from a single hopper for seven different particle sizes. Here, the hopper outlet thickness is fixed at $D = 5$ mm. The three particle curtain regions can be generally identified for all particle sizes. The largest particle size can be seen to exhibit the overall greatest curtain thickness in the upper half of the curtain. However, the flow features including the neck zone and the unsteady fluctuation of the curtain in the intermediate field become more prominent with a decrease in particle size. Particles with different size distributions, shapes (e.g. oval or cylindrical) and surface roughnesses will inevitably influence the particle flow behaviour via particle–wall, particle–particle and particle–fluid interactions, causing differences in the values. However, the controlling scaling relationship will still apply, albeit with somewhat different values.

Figure 7. Instantaneous side-view images of the particle curtain $0 \lesssim x \lesssim 940$ mm for seven different particle sizes ($\skew6\bar{d}_p$). Here, the hopper outlet thickness ($D$) is constant at 5 mm. The particle flow is from top to bottom following the gravitational direction ($g$).

Figure 8 presents the time-averaged measurements of the axial evolution of the particle curtain thickness ($\varDelta _{0.5}$) as a function of falling distance for seven particle sizes ($\skew6\bar{d}_p = 163$, 192, 216, 257, 363, 399 and $500\,\mathrm {\mu }$m) and for two hopper outlet thicknesses ($D = 3$ and 8 mm). Note that reliable data are not available immediately below the hopper exit, $0 \leq x/D \lesssim 1$, due to the parallax error, i.e. the edge of the hopper closer to the backlight obscures light through a small region of the flow in front of it near to the exit. Nevertheless, the results show that, for all of these particles and hopper outlet thicknesses, there is a strong increase in $\varDelta _{0.5}$ with downstream distance in the near-field expansion region ($0 \lesssim x \lesssim 75$ mm), due to strong particle dispersion away from the centreline, consistent with the results presented in figure 6. The presence of this near-field expansion was also observed in similar studies of particles falling from a round hopper with $4\,{\rm mm}\leq D \leq 20$ mm and $40\,\mathrm {\mu }{\rm m}\lesssim d_p \lesssim 1000\,\mathrm {\mu }$m (Cooper & Arnold Reference Cooper and Arnold1995; Ogata et al. Reference Ogata, Funatsu and Tomita2001; Prado et al. Reference Prado, Amarouchene and Kellay2013). The reason for this phenomenon is yet to be confirmed, but we postulate that it is due to particle collisions near to the hopper exit. More specifically, particle dispersion away from the centreline could be attributed to collisions between the faster moving particles near to the hopper centreline and the slower moving particles near to the hopper wall.

Figure 8. Axial variation of the free-falling particle curtain thickness ($\varDelta _{0.5}$), defined as the full width at half-maximum, measured for seven different particle diameters for cases of (a) $D = 3$ mm and (b) $D = 8$ mm. Note that there are two labels for the $y$-axes, one for the non-normalised curtain thickness ($\varDelta _{0.5}$), and the other for the curtain thickness normalised by hopper outlet thickness ($\varDelta _{0.5}/D$). The circles and x symbols denote the locations of the local minima and maxima of $\varDelta _{0.5}$ for each case.

Interestingly, the initial rate of expansion (i.e. the gradient of the curves in figure 8) in the near field is approximately the same for all particle sizes and for any given hopper outlet thickness. This suggests that the mechanism that causes the particles to disperse laterally away from the hopper centreline in the near field is relatively insensitive to particle size. Furthermore, the larger particles disperse further away (laterally) from the curtain centreline in the near field than do the smaller particles for all hopper outlet thicknesses. This is attributed to the greater inertia of the larger particles. That is, larger particles will be transported further away from the centreline than smaller particles, due to their greater lateral momentum relative to the drag force.

Further downstream from the initial expansion region, for all cases except for the largest particles ($\skew6\bar{d}_p = 500\,\mathrm {\mu }$m) issuing from the smallest hopper ($D = 3$ mm), the curtain thickness exhibits a local maximum (a peak). The location of this local maximum, $x_{pk}$ is approximately within the region $5 \lesssim x_{pk}/D \lesssim 26$ for all cases, with this location shifting further downstream as the particle size is increased. This is consistent with the larger particles having greater momentum (relative to drag) than the smaller particles. The prominence of the local maximum appears to decrease with an increase in particle diameter and a decrease in hopper outlet thickness, i.e. for greater values of $\skew6\bar{d}_p/D$. For the case where $\skew6\bar{d}_p/D$ is the greatest, i.e. $\skew6\bar{d}_p/D =1/6$, which occurs with the largest particles ($\skew6\bar{d}_p = 500\,\mathrm {\mu }$m) and the smallest hopper ($D = 3$ mm), the peak disappears completely, as previously noted. That is, this case does not exhibit a neck zone.

Where a local maximum occurs, it is followed by a region in which the curtain thickness decreases with downstream distance, before increasing again. This is the well-known ‘necking’ effect (Hardalupas, Taylor & Whitelaw Reference Hardalupas, Taylor and Whitelaw1989; Lau & Nathan Reference Lau and Nathan2014, Reference Lau and Nathan2016; Ho et al. Reference Ho, Christian, Romano, Yellowhair, Siegel, Savoldi and Zanino2017), where the streamwise acceleration of particles becomes significant, leading to a decreased local pressure within the curtain. This, in turn, results in particles (and fluid) being drawn towards the centreline. The curtain thickness typically continues to decrease for a few hundred millimetres downstream until $40D \lesssim x \lesssim 80D$ where $\varDelta _{0.5}/D$ reaches a local minimum (a trough), before it begins to expand again. The latter marks the beginning of the intermediate-field expansion region. That is, downstream of the local minimum, the curtain thickness expands, presumably due to the entrainment of surrounding gas into the curtain, consistent with the well-known dynamics of particle-laden jets (Longmire & Eaton Reference Longmire and Eaton1992; Lau & Nathan Reference Lau and Nathan2016). It should be noted that, in the intermediate field, the gravity-driven acceleration of the particles is expected to be relatively small, with the drag forces dominating over the gravity-driven forces such that the curtain approximates a jet. Near to the end of the measurement region, the expansion of the curtain approaches a linear function, with the rate of expansion (i.e. the gradient of the curves in figure 8) increasing as the particle diameter is decreased. This is because smaller particles, with their lower inertia, are expected to have a stronger interaction with the ambient air, causing smaller particles to spread more rapidly than the larger particles (Lau & Nathan Reference Lau and Nathan2016).

The results also show that the maximum curtain thickness in the near-field expansion region, $\varDelta _{0.5, pk}/D$, typically decreases with an increase in the hopper outlet thickness. For example, in the cases of $\skew6\bar{d}_p = 163\,\mathrm {\mu }$m particles, $\varDelta _{0.5, pk}/D \approx 1.85$ for the hopper of $D = 3$ mm, decreasing to $\varDelta _{0.5, pk}/D \approx 1.61$ for the hopper of $D = 8$ mm. This suggests that a normalisation of the curtain thickness by the hopper outlet thickness does not result in the self-similarity profiles in flow thickness that are observed in jets (Mi, Nobes & Nathan Reference Mi, Nobes and Nathan2001). This difference is not surprising, but implies that the mechanisms controlling the initial curtain expansion does not scale with $D$ alone.

Figure 9 presents an example of the curtain side view illustrating the boundaries of the dilute-particle region and the curtain core defined by 5 % and 50 % of the peak opacity, i.e. $\varDelta _{0.05}$ and $\varDelta _{0.5}$. The dilute-particle region is clearly differentiated from the curtain core by a much lower volume fraction. The radial profiles do not vary very much with axial distance, which implies that the radial movement of the particles is weak. The 5 % contour is well above the noise threshold at $|y/D| > 4$, giving confidence in the choice of boundary. These boundaries are also reliable for $|y/D| < 1$.

Figure 9. The contours of the half-width, $\varDelta _{0.5}$ (red line) and the contours of 5 % of the peak opacity $\varDelta _{0.05}$ (green line), defining the boundaries of the dilute-particle region, as presented in (a), the axial evolution and (b) the radial profile at $x/D = 40$. Here, $D = 5$ mm and $\skew6\bar{d}_p = 270\,\mathrm {\mu }$m.

Figure 10 presents the time-averaged measurements of the axial evolution of total thickness of the dilute-particle region, accounting for both sides of the curtain core region for seven particle sizes ($\skew6\bar{d}_p = 163$, 192, 216, 257, 363, 399 and $500\,\mathrm {\mu }$m) and for two hopper outlet thicknesses ($D = 3$ and 8 mm). This dilute-particle region is defined by a difference of curtain thicknesses between 5 % and 50 % of the peak opacity at the lateral plane, that is $\varDelta _{0.05} - \varDelta _{0.5}$. The results show that $\varDelta _{0.05} - \varDelta _{0.5}$, for $x/D < 10$, increases from a non-zero value near to the hopper outlet to a value some two times the hopper outlet thickness, which is of the same magnitude as that for $\varDelta _{0.5}$, as shown in figure 8. Furthermore, the largest particles expand the most. This indicates that the dilute-particle region originates from the near-exit region, where inter-particle collisions and particle–hopper interactions generate significant particle lateral momentum, causing the initial dispersion, as previously discussed. That the outer layer of the particles is not entrained back into the curtain core implies that the initial expansion moves them into a region of negligible radial pressure gradient, including that induced by particle acceleration in the neck zone of the curtain core (figure 8). That is, these particles are not re-entrained into the curtain.

Figure 10. Axial evolution of the dilute-particle region characterised by the difference of curtain thicknesses $\varDelta _{0.05} - \varDelta _{0.5}$, normalised by hopper outlet thickness $D$, for cases of (a) $D = 3$ mm and (b) $D = 8$ mm with seven different particle sizes ($\skew6\bar{d}_p$).

For $x/D > 130$ ($D = 3$ mm) and 60 ($D = 8$ mm), there is a second noticeable increase of the regional thickness for the smallest particles. This coincides with the region of intermediate-field expansion as identified in figure 8, with a higher expansion rate for the smaller $\skew6\bar{d}_p$. The consistency in particle expansion between these two lateral regions implies a common controlling mechanism – namely that both are in the drag-dominated regime where smaller particles have a higher drag per unit volume ratio, therefore expanding faster. In contrast, the larger particles, $\skew6\bar{d}_p = 257 \sim 500\,\mathrm {\mu }$m, do not show the same expansion. This suggests the persistence of the upstream mechanisms throughout the downstream region.

Figure 11 presents a schematic diagram of the side view of the falling particle curtain, together with the notation used to characterise the near- and intermediate-field expansions in the flow. The magnitudes of the local ‘peak’ and ‘trough’ features of the curtain thickness are denoted by $\varDelta _{0.5,pk}$ and $\varDelta _{0.5,tr}$, respectively, with their corresponding streamwise locations measured from the hopper exit represented by $x_{pk}$ and $x_{tr}$, respectively. The gradients of expansion in the near field and the intermediate field are denoted separately by $K_{near}$ and $K_{int}$, with both defined as the differentiation of the full curtain thickness over a specific streamwise location defined as appropriate below.

Figure 11. Schematic diagram of the curtain side view, together with the notation used to describe the key features of the local peak (subscript ‘pk’) and trough (subscript ‘tr’) that characterise the curtain thickness $\varDelta _{0.5}$. The gradients of the expansion in the near field and the intermediate field are denoted by $K_{near}$ and $K_{int}$. The ‘dilute-particle’ region with particles dispersed around the denser curtain ‘core’ region, as denoted by $\varDelta _{0.05} - \varDelta _{0.5}$, is also shown.

This image also reveals more details about the source of the particles that are found at low number density in the outer ‘dilute-particle’ region that surrounds the dense particle curtain, denoted by $\varDelta _{0.05} - \varDelta _{0.5}$. The generation of this region of more dispersed particles is important because they are likely contributors to the generation of fine particle emissions that can readily migrate from a curtain by being entrained into secondary flows within a falling particle receiver or other types of equipment. This image suggests that these particles originate from the near-field expansion region, where the collision-induced force disperses them beyond the range at which they can be re-entrained into the jet.

Figure 12 presents the gradient of the near-field expansion ($K_{near}$) and the intermediate-field expansion ($K_{int}$) as a function of normalised particle diameter, $\skew6\bar{d}_p/D$, for all examined particle diameters and hopper outlet thicknesses. These values were obtained by curve fitting a linear function to the data in the near-field expansion region ($0 \lesssim x/D \lesssim 13$), and the downstream region of the intermediate field ($76 \lesssim x/D \lesssim 300$) where the evolution of $\varDelta _{0.5}/D$ with axial distance is approximately linear. For all cases, the linear curve fits have a coefficient of determination, $R^2 \approx 95\,\%$. The results from figure 12(a) show that $K_{near}$ is typically an order of magnitude larger than $K_{int}$ from figure 12(b) for all cases. The near-field expansion gradients found here, $0.0796 \leq K_{near} \leq 0.1703$, translate to a near-field expansion half-angle in the range $2.3^{\circ } \leq \alpha \leq 4.8^{\circ }$. The near-field expansion rate was also found to increase with hopper outlet thickness for almost all normalised particle diameters ($\skew6\bar{d}_p/D$), while only being weakly dependent on the particle diameter, consistent with the measurements (figure 8). Although particle initial momentum is based on particle inertia, which depends on particle diameter, this finding suggests that particle lateral momentum is gained mostly from the increase in $D$. This increases the fraction of particles within the flow that are unaffected by drag at the wall, thereby increasing the particle axial momentum and, in turn, the lateral momentum that arises from inter-particle collisions in the near field. Such collisions occur because the particles near to the hopper internal wall fall more slowly than those closer to the centreplane. In addition, a higher particle volume fraction at the exit was found to be associated with a larger $D$ (Janda et al. Reference Janda, Zuriguel and Maza2012; Wang et al. Reference Wang, Yang, Zhang, Cai, Lin, Chen and Yang2020), which also increases the near-field expansion of the particle curtains by enhancing particle–particle collisions.

Figure 12. The dependence of the linear gradients of (a) the near-field expansion, $K_{near}$ and (b) the intermediate-field expansion, $K_{int}$ on the normalised particle diameter ($\skew6\bar{d}_p/D$) for a series of hopper outlet thicknesses ($D = 3$, 4, 5 and 8 mm). Gradients are obtained through curve fitting of the curtain thickness axial profiles (see figures 8 and 11).

Interestingly, the rate of expansion in the intermediate field shows the opposite trend to that in the near field, with the intermediate-field expansion rate increasing as both the hopper outlet thickness and the particle diameter are decreased, as shown in figure 12(b). This implies a different controlling mechanism of the particle curtain dynamics in this region compared with that in the near field. The general trend of decreased values of $K_{int}$ with an increase in $D$ is consistent with the Beverloo equation, which shows that the mass flow rate, and hence momentum, scales with $D^{1.5}$. That is, the influence of particle inertia is greater for larger hopper outlet thickness. This influence may also be compounding, since a lower spreading rate will also lead to higher local volume fractions. Previous work has shown that, where the volume fraction is sufficiently high, the presence of downstream particles can reduce the drag of the upstream particles due to the effect of aerodynamic wake (Cooper & Arnold Reference Cooper and Arnold1995; Wang et al. Reference Wang, Ren, Zhao, Chu, Cao, Yang, Duan, Fan and Qu2016; Sedaghatizadeh et al. Reference Sedaghatizadeh, Arjomandi, Lau and Nathan2022). Here, we call it the ‘particle group wake’ (or the particle ‘group effect’) to highlight its bulk effect on the particle stream. This, in turn, will further reduce the expansion rate due to the reduced entrainment into the curtain. This trend has also been found to occur further downstream for a particle stream issuing from a round hopper, although the flow regimes are significantly different from those identified in the present study (Wang et al. Reference Wang, Ren, Zhao, Chu, Cao, Yang, Duan, Fan and Qu2016). In addition, the results also show that $K_{int}$ increases with a decrease in $\skew6\bar{d}_p$, with $K_{int}$ exhibiting a stronger dependence on particle size than $K_{near}$. This highlights the importance of particle–air interaction through the shear force between the particles and the surrounding fluid, which dominates the intermediate-field expansion in a jet-like flow regime, consistent with the results presented in figure 8.

Figure 13 presents the relative axial streamwise positions of the local peaks and troughs in the curtain thickness profiles, $x_{pk}/D$ or $x_{tr}/D$, and their relative magnitudes, $\varDelta _{0.5,pk}/D$ or $\varDelta _{0.5,tr}/D$, as identified from figure 8 and illustrated in figure 11, as a function of normalised particle size ($\skew6\bar{d}_p/D$). In particular, figure 13(a) shows that the value of $x_{pk}/D$ for any given hopper outlet thickness increases monotonically with the particle size. This is consistent with the momentum of particles relative to the drag force being greater for larger particles, as previously discussed. The influence of particle diameter on $x_{pk}/D$ can be illustrated using a simplified model comprising a single spherical particle in a quiescent environment imparted with an initial radial (or lateral) velocity in the Stokes regime. Under these conditions, the equation of motion for the particle is

(4.3)\begin{equation} \frac{\partial V_p}{\partial t} + \frac{V_p}{\tau_p}=0, \end{equation}

where $V_p$ is the particle velocity in the radial direction (i.e. the $y$-direction). Solving this equation for particle radial displacement results in $y_p (t)=\tau _p V_0 [1-\exp ((-t)/\tau _p )]$, where $y_p$ is the displacement of the particle from the centreplane and $V_0$ is the particle initial velocity in the radial direction. Therefore, the radial distance required for the particle to reach a radial velocity component of zero is ${\propto }\tau _p$. Since $\tau _p\propto \skew6\bar{d}_p ^2$, this implies that $y_p$ and hence $x_{pk}/D$ are also positively correlated with $\skew6\bar{d}_p^2$ ignoring any differences in $K_{near}$ for different particle sizes shown in figures 8 and 12.

Figure 13. The influence of normalised particle diameter ($\skew6\bar{d}_p/D$) on the streamwise locations of the local peaks, $x_{pk}/D$, and troughs, $x_{tr}/D$, in the curtain thickness axial profiles (see also figures 8 and 11), together with the magnitudes of the peaks, $\varDelta _{0.5,pk}/D$ and troughs, $\varDelta _{0.5,tr}/D$, for a series of hopper outlet thicknesses ($D = 3$, 4 ,5 and 8 mm).

Figure 13(a) also shows that the profiles of $x_{tr}/D$ for the four different hopper outlet thicknesses collapse quite well to a single arch-shaped curve, with the maximum $x_{tr}/D$ appearing for $\skew6\bar{d}_p/D \approx 0.05$, although the flow inside the hopper was previously found to be independent of particle size for $\skew6\bar{d}_p/D < 0.05$ (Beverloo et al. Reference Beverloo, Leniger and Van de Velde1961). As the neck zone is formed by particle–air interaction, the local ‘trough’ of the particle curtain thickness established at the end of this region depends on the balance between particle inertia, which drives the near-field expansion, and the influence of the entrained surrounding air induced by particle acceleration that drives particles toward the curtain ‘core’. That is, the larger particles in the range $\skew6\bar{d}_p/D > 0.05$ have a higher ratio of inertia to drag, so that they tend to retain their original trajectory to disperse away from the curtain core, while the relatively smaller hopper outlet thickness within $\skew6\bar{d}_p/D > 0.05$ cannot generate sufficient pressure difference in the fluid phase between curtain core and surroundings to entrain these larger particles back to the curtain core. This results in lower values of $x_{tr}/D$ and of $\varDelta _{0.5,tr} /\varDelta _{0.5,pk}$, indicating that particles in the range $\skew6\bar{d}_p/D > 0.05$ become less responsive (less necking) to the pressure gradient as particle size increases. Indeed, no neck zone at all is formed for the extreme case of $D = 3$ mm and $\skew6\bar{d}_p = 500\,\mathrm {\mu }$m, as can be seen in figure 8(a). In contrast, the smallest particles exhibit the opposite trend due to their stronger dependence on aerodynamic drag, leading to the most significant necking effects. That is, $x_{tr}/D$ decreases with an increase in $\skew6\bar{d}_p/D$ for $\skew6\bar{d}_p/D > 0.05$.

While particle accelerations were not directly measured in the present study, they can be estimated using the Stokes assumption (3.5). The calculated accelerations for $D = 3$ mm, with seven different values of $\skew6\bar{d}_p$ are shown in figure 14. As can be clearly seen, the particle acceleration decreases with an increase in $x$ and a decrease in $\skew6\bar{d}_p$. The negative correlation of $x_{tr}/D$ with particle size for $\skew6\bar{d}_p/D > 0.05$ can therefore be attributed to the increase in particle acceleration with particle diameter, owing to their stronger dependence on gravity than on drag, resulting in a greater acceleration closer to the exit plane.

Figure 14. The axial evolution of the particle axial acceleration calculated with the simplified analytical model in the Stokes regime (3.5) for seven different particle sizes and for $D = 3$ mm.

For $\skew6\bar{d}_p/D < 0.05$, the particle acceleration plausibly generates a lateral pressure gradient of sufficient magnitude for the drag of the entrained fluid to become significant relative to particle inertia. As a result, the larger particles of $\skew6\bar{d}_p/D < 0.05$ are no longer able to follow their original spreading trajectory but rather move towards to the curtain core. This can be directly observed from figure 8(b), where a neck zone is even generated for $\skew6\bar{d}_p = 500\,\mathrm {\mu } {\rm m}$. However, the process of contraction requires a longer distance and is also less prominent for the larger particles, due to their lower ratio of drag to inertia. Therefore, the neck zone in the range of $\skew6\bar{d}_p/D < 0.05$ can be deduced to be dominated by hopper-controlled particle acceleration, which induces significant pressure gradient.

Figure 13(b) presents the relative magnitudes of the thickness of the peaks (local maxima) and troughs (local minima), $\varDelta _{0.5,pk}/D$ or $\varDelta _{0.5,tr}/D$, as a function of $\skew6\bar{d}_p/D$ for $D = 3$, 4, 5 and 8 mm. The results show that, for all cases, $\varDelta _{0.5,pk}/D$ is significantly greater than the hopper outlet thickness, with $1.61 \leq \varDelta _{0.5,pk}/D \leq 2.54$ across all examined conditions. That is, the actual thickness of the particle curtain can be up to a factor of 2.5 greater than the hopper outlet thickness, even in the near field. Furthermore, the results show that the difference between $\varDelta _{0.5,pk}$ and $\varDelta _{0.5,tr}$ decreases with an increase in $\skew6\bar{d}_p/D$. That is, the near-field expansion and necking are more prominent with larger hopper outlet thickness and smaller particle sizes, consistent with figure 8. For this reason, an increase in $\skew6\bar{d}_p/D$ also causes the values of $\varDelta _{0.5,pk}$ to approach $\varDelta _{0.5,tr}$ and, for sufficiently large $\skew6\bar{d}_p/D$, the values of $\varDelta _{0.5,pk}$ and $\varDelta _{0.5,tr}$ may eventually converge, at which point the neck zone ceases to exist, as can also be seen in figure 8(a). Indeed, the measurements of $\varDelta _{0.5,pk}$ and $\varDelta _{0.5,tr}$, together with $x_{pk}/D$ and $x_{tr}/D$, for the smallest hopper outlet thickness of $D = 3$ mm and the largest particle size of $\skew6\bar{d}_p = 500\,\mathrm {\mu } {\rm m}$ are not presented in figure 13 because no peak or trough exists for this case.

The results also show that $\varDelta _{0.5,pk}/D$ increases approximately linearly with $\skew6\bar{d}_p/D$ and the profiles for all four different hoppers collapse well with these axes. The increase in $\varDelta _{0.5,pk}/D$ with particle diameter is attributed to the greater momentum (relative to drag) of the larger particles, as previously discussed. Together with the finding from figure 13(a), it can be deduced that both the location and the magnitude of the peaks at the near-field expansion scale with the particle size and hopper outlet thickness, potentially highlighting the influence of the particle–particle and particle–wall interactions as being the dominant mechanisms controlling the flow dynamics in the near-field expansion region. This finding of the dominant mechanisms is consistent with the trends identified previously for the influence of hopper on the near-field flow region using a round hopper outlet (Wang et al. Reference Wang, Ren, Zhao, Chu, Cao, Yang, Duan, Fan and Qu2016).

Figure 14 presents the evolution of particle axial acceleration ($a_{p,SD}$) calculated with the simplified analytical model in the Stokes regime (3.5). The smaller particles present a lower $a_{p,SD}$ throughout the examined falling distance, although the starting $a_{p,SD}$ is similar for the various $\skew6\bar{d}_p$. The smaller particles also undergo a faster decay of particle acceleration in the axial direction, decreasing to more than half of the starting $a_{p,SD}$ at the first 100$D$, due to their higher drag per unit volume ratio. In contrast, the largest particle $\skew6\bar{d}_p = 500\,\mathrm {\mu } {\rm m}$ continues preserving the high acceleration due to a greater influence of gravitational acceleration.

4.3. Particle curtain transmittance

Figure 15 presents seven representative instantaneous images of measured curtain transmittance, each with a different particle diameter, here all obtained with hopper $D = 5$ mm. A qualitative trend can be seen in which the curtain transmittance is low near to the exit and increases with downstream distance. This overall trend is consistent with previous investigations (Ogata et al. Reference Ogata, Funatsu and Tomita2001; Sedaghatizadeh et al. Reference Sedaghatizadeh, Arjomandi, Lau and Nathan2022), and is due to the particle velocity increasing with downstream distance (due to the acceleration of gravity), resulting in a decrease in local particle number density with streamwise distance and an increase in transmittance. However, the higher resolution of the present images provides additional detail over what has been reported previously. In particular, it can be seen that the distance to the streamwise location at which the opacity begins to fall decreases with an increase in particle diameter. Furthermore, the region downstream from this location exhibits significant instantaneous structure, which can be deduced to be driven by large-scale flow instabilities.

Figure 15. Instantaneous front-view images of the particle curtain issuing from a rectangular hopper of $D = 5$ mm, recorded at a falling distance between 0 and $\approx$900 mm from the hopper exit, for seven different particle diameters ($\skew6\bar{d}_p$). The positions of $x_{pk}$ and $x_{tr}$ (see also figure 11) are shown with red and blue lines, respectively. Here, the particle flow is from top to bottom following the direction of gravity.

The overall trend of a decrease in the distance at which transmittance begins to increase with particle diameter is consistent with the increase in particle acceleration (see figure 14). That is, the axial spacing between particles depends on acceleration, which increases with $\skew6\bar{d}_p$. This is also consistent with the previous measurements using hoppers with rectangular outlets, which show that particle velocity increases with $\skew6\bar{d}_p$ and with axial distance (Kim et al. Reference Kim, Siegel, Kolb, Rangaswamy and Moujaes2009; Sedaghatizadeh et al. Reference Sedaghatizadeh, Arjomandi, Lau and Nathan2022). In addition, the particle cross-sectional area per unit volume also decreases with an increase in $\skew6\bar{d}_p$, which will lead to a greater $T$ for a given hopper outlet thickness.

The position of the local maximum in curtain thickness at the end of the near-field expansion region, $x_{pk}/D$, and the location of the local minimum curtain thickness at the end of the neck zone, $x_{tr}/D$ (see also figure 11), are also shown in figure 15 (with the red and blue lines, respectively). It can be seen that both of these features occur well upstream from the point at which the opacity begins to fall noticeably, with no clear correlation between the locations of $x_{pk}/D$ and $x_{tr}/D$ and local values of light transmittance. The lack of a direct influence is because these features cause a difference in lateral spacing (along the $y$-axis) between particles, rather than axial spacing, while the transmittance of the particle curtain is mostly dependent on the axial spacing between particles.

The instantaneous images presented in figure 15 also reveal the presence of highly localised regions of low curtain transmittance, with these regions typically consisting of filament-like or rope-like structures that have long axial coherence. These structures are typically aligned in the streamwise direction, with their span being of the same order as the local thickness of the curtain. This is consistent with the presence of particle ‘clusters’, which have been demonstrated to exist in a range of particle-laden jets (Lau & Nathan Reference Lau and Nathan2017; Lau, Frank & Nathan Reference Lau, Frank and Nathan2019). Here, the term ‘cluster’ is used following the pioneering work of Longmire & Eaton (Reference Longmire and Eaton1992) and Eaton & Fessler (Reference Eaton and Fessler1994) to refer to local groups of particles whose spacing is significantly less than average due to the role of large-scale flow features, but does not imply that any particles actually touch each other to form aggregates. These clusters are typically characterised by highly localised regions of high particle number density within the flow, and are typically caused by the interaction of the particles with the fluid flow (Maxey Reference Maxey1987; Squires & Eaton Reference Squires and Eaton1991; Longmire & Eaton Reference Longmire and Eaton1992; Eaton & Fessler Reference Eaton and Fessler1994; Capecelatro, Desjardins & Fox Reference Capecelatro, Desjardins and Fox2015). Whatever their origin, the effect of these clusters is to induce local, instantaneous non-uniformity in transmittance, that is to both reduce transmittance locally where the clusters are present, and to increase the transmittance on either side of the clusters. The presence of these clusters is particularly significant for applications such as FPRs, as the increased transmittance in non-clustered regions may result in solar radiation passing through the curtain, potentially damaging internal walls and surfaces. Additionally, the particles within the clusters may be shadowed from solar radiation due to their high opacity, resulting in lower average particle temperatures. These images also highlight the need to better understand the origins of these clusters, the need to be able to identify the axial distance at which they become significant and the need for methods to manage or account for their influence. Furthermore, the clear presence of these particles in almost all cases also implies that current steady-state models, both analytical and numerical, are likely to be unreliable in regions where clustering is significant.

Figure 16 presents the time-averaged curtain transmittance as a function of the falling distance for seven different particle sizes and for a fixed hopper outlet thickness $D = 3$ mm. The results show that the slope of the initial rise of transmittance depends strongly on the particle diameter. For the largest particle, this slope is high, so that $T$ increases strongly within the first 20$D$. In contrast, the initial slope is near to zero for the smallest particles. These trends are consistent with the visual observations shown in figure 15 and the accompanying discussion, with more details to support these deductions being derived from the analysis below. The slopes of the transmittance axial profiles reduce beyond a falling distance of approximately 50$D$, with monotonically consistent trends (almost parallel curves) in the relationship between $T$ and $x/D$ for all cases of different $\skew6\bar{d}_p$ at the end of the measurement window. This former trend can be explained by a much stronger decrease in particle acceleration with axial distance for smaller particles than larger particles, when approaching the terminal velocity (see figure 14). That is, the rate of acceleration decays much faster for smaller particles than larger ones, which reduces the rate of increase in particle vertical separation and the rate of increase in transmittance. However, the latter trend regarding the near-parallel and asymptotically stable profiles of $T$ between different $\skew6\bar{d}_p$ can be attributed to the relatively lower decay rate of particle acceleration. This trend also suggests a simpler regime for the curtain transmittance than for the curtain expansion under the effect of $\skew6\bar{d}_p$, which can be further discovered by a normalisation of the data.

Figure 16. Experimental measurements of time-averaged particle curtain transmittance as a function of particle falling distance, $x$, normalised by the hopper outlet thickness ($D$) for seven different particle diameters. In this case, the hopper outlet thickness was fixed at $D = 3$ mm.

Figure 17 presents the same particle curtain transmittance as plotted in figure 16, except that the data are now plotted as a function of the particle falling distance that is normalised by the hopper outlet thickness and particle Froude number based on the estimated particle terminal velocity (using the Stokes-drag model, as discussed in § 3.2), $Fr_t = U_{t,SD}/(gD)^{0.5}$. The data are only plotted for $D = 3$ mm because of the most dispersed particle condition, which best suits the adopted model. It can be seen that this normalisation provides a good collapse of the cases with different values of $\skew6\bar{d}_p$ onto a single curve. This highlights the dependence of the transmittance on the gravity-driven evolution of particle vertical velocity. It is also consistent with previous findings that the magnitude of particle streamwise velocity depends on $Fr_t$ (Sedaghatizadeh et al. Reference Sedaghatizadeh, Arjomandi, Lau and Nathan2022). The collapse with the current normalisation also provides evidence that the balance between the inertial and gravitational forces dominates the evolution of particle vertical velocity rather than any influences of the shear forces in the shear layer, at least in this set comparing different values of $\skew6\bar{d}_p$.

Figure 17. Measured particle curtain transmittance as a function of particle falling distance ($x$) normalised by the hopper outlet thickness ($D$) and the terminal Froude number ($Fr_t$) for seven particle sizes. The value of $Fr_t$ was calculated using the particle terminal velocity estimated from the analytical model in the Stokes regime, $U_{t,SD}$ (3.4).

Figure 18 presents the measured values of time-averaged curtain transmittance ($T$) as a function of the falling distance for four different hopper outlet thicknesses ($D = 3$, 4, 5 and 8 mm), but with a fixed particle size $\skew6\bar{d}_p = 363\,\mathrm {\mu }$m. This shows that, somewhat analogously to figure 16, the initial rate at which $T$ increases also varies greatly, from a high slope for $D = 3$ mm, to a near-zero initial slope for $D = 8$ mm. Taken together with figures 15 and 16, this shows that the axial distance to the point at which the transmittance starts to increase rapidly with axial distance increases with both an increase in $D$ and a decrease in $\skew6\bar{d}_p$. This is partially attributed to the higher particle volume fraction for using a greater $D$, as discussed before. That is, using a hopper with a larger $D$ increases the particle volume fraction significantly at the near-exit region which, in turn, causes a dramatic difference in the axial profile of curtain transmittance. However, this difference reduces with the axial distance as particles fall downwards, influenced by varying particle–fluid interactions depending on the size of $D$. The bulk-average particle–fluid interaction with curtain expansion is expected to be stronger for a smaller $D$ due to the reduced influence of the particle group effect. This results in a lower particle acceleration and a smaller particle vertical separation at a given downstream location, thereby leading to a faster reduction of the rising rate in $T$ for a smaller $D$. Nevertheless, the transmittance of a curtain with a larger $D$ does not exceed that for a smaller $D$ at a given downstream location. For example, $T = 0.66$ at $x = 900$ mm for $D = 3$ mm, which is some six times greater than $T = 0.11$ for $D = 8$ mm at the same location. This is because an increase in $D$ also increases the number of particles at any given falling distance, which decreases $T$. Noting also that, according to the Beverloo equation the mass flow rate through the hopper (per unit width) scales with $D^{1.5}$, from (3.8). That is, $\bar \theta _0\ \rho _p \ \forall _p\ \bar V_0 \propto D^{0.5}$. For any given particle size, this implies that an increase in $D$, causes an increase in the initial particle velocity, and/or the initial particle concentration at the exit. The latter results in a decrease in transmittance in the near field, as measured in figure 15. A complete summary of measured curtain thickness and transmittance, along with details of the hopper geometry, is provided in the Supplementary Material available at https://doi.org/10.1017/jfm.2024.896.

Figure 18. Experimental measurements of time-averaged particle curtain transmittance as a function of normalised particle falling distance, $x/D$, for four different hopper outlet thicknesses ($D = 3$, 4, 5 and 8 mm) and for a fixed particle size $\skew6\bar{d}_p = 363\,\mathrm {\mu }$m.

4.4. Insights into curtain evolution from analytical models

Figure 19 presents a comparison of the measured axial evolution of a particle curtain transmittance together with that estimated from the Beer–Lambert model (3.9), for both the zero-drag (3.2) and Stokes-drag (3.4) models. As noted in § 3, these integrated models represent two extreme scenarios, which were developed to support the interpretation of the measured results, most notably, the axial evolution of the particle dynamics and its correlation with curtain transmittance. The results are shown for three different values of particle sizes (i.e. 163, 257 and $500\,\mathrm {\mu }$m) and a single hopper outlet thickness of $D = 5$ mm. It can be seen that the transmittance calculated with the zero-drag model is always higher than that calculated with the Stokes-drag model. This is consistent with the understanding that the axial particle acceleration is greater for the case with no drag, leading to a greater axial separation between particles, and a higher void fraction. It can also be seen that the difference between the inferred transmittance profile for the Stokes-drag and the zero-drag models decreases with an increase in particle diameter. This can be explained by the drag force per unit weight of a particle decreasing with an increase in $\skew6\bar{d}_p$. That is, the influence of air drag decreases with an increase in $\skew6\bar{d}_p$. It can also be seen that both models tend to over-predict the transmittance in the near field and under-predict it at the end of the measured region in the intermediate field, particularly for the smallest particles. This is consistent with the particle dynamics in the near field being controlled by non-aerodynamic mechanisms of particle–particle interactions, as expected for a flow that originates in the granular flow regime. This granular regime decreases particle acceleration even beyond that for the Stokes regime. Furthermore, the measured transmittance at the end of the measurement region is greater than predicted with the zero-drag model, particularly for the smallest particles. This either implies that the particles accelerate faster than those with zero drag, which is implausible, or that the assumptions inherent in the current transmittance model, such as the assumption that particle scattering, particle shadowing and multiple scattering effects are negligible, and that the particle number density and velocity across the curtain at any given falling distance is uniform, begin to break down. Nevertheless, the comparison between different scenarios and associated mechanisms revealed herein also indicates that the particle dynamics predicted through the use of the Beverloo equation in these models is reasonable and consistent, thereby proving its validity in scaling assessments of flows as discussed for figure 5.

Figure 19. Comparison of the measured and modelled axial variation of transmittance for a particle curtain issuing from a hopper with different particle sizes ($\skew6\bar{d}_p = 163$, 257 and $500\,\mathrm {\mu }$m) and a fixed hopper outlet thickness of $D = 5$ mm. The figure includes the transmittance calculated using two models, namely the zero-drag model (§ 3.1) and the Stokes-drag model (§ 3.2).

It can also be seen from figure 19 that the zero-drag model predicts the measured transmittance quite well for the largest particle size, consistent with the analysis of Sedaghatizadeh et al. (Reference Sedaghatizadeh, Arjomandi, Lau and Nathan2022). This indicates that these larger particles approximate particles with zero drag, which is attributed to the decrease in drag per unit volume for larger $\skew6\bar{d}_p$ and the group effect whereby leading particles shield the upstream particles from the effects of the surrounding fluid, causing them to fall faster than if they were in a quiescent, unperturbed, environment.

Given that the calculated curtain transmittance from the zero-drag model matches the experimental data more closely than the Stoke-drag model, but demonstrated with cases comprising only three particle sizes and one hopper outlet thickness, it would be insightful to expand the range of cases to fully examine its dimensionless correlations between all particle sizes and all hopper outlet thicknesses. In addition, it is also important to understand that the closest match for the largest particle size is primarily due to the effect of particle size rather than the hopper outlet thickness. To provide clear examinations, a special approach of rearranging the integrated analytical model into a linear relationship incorporating the Beverloo equation, the Beer–Lambert law and the zero-drag scenario, has been taken as follows.

The transmittance model can be further simplified by assuming the particles have zero drag. Combining (3.2) (particle velocity under the zero-drag assumption) and (3.9), we obtain

(4.4)\begin{equation} T(x)=\exp \left[-\frac{3 \dot M_p}{2 \rho_p d_p (U_0^2+2gx)^{0.5}} \right]. \end{equation}

If the particle-phase density (not to be confused with the particle material density, $\rho _p$) at the hopper exit is assumed to be the same as that within the hopper (i.e. equal to the particle bulk density $\rho _b$), then a calculation of the mass flow rate through the hopper exit results in a particle initial velocity of $U_0= (\dot M_p)/(\rho _b D)$. Now, utilising the Beverloo equation (4.1), and assuming that the $d_p/D$ term is small, (4.4) becomes

(4.5)\begin{equation} \frac{C_1^2}{2}\left[\frac{1}{(\kappa \ln{T(x)})^2} -1 \right]=\frac{x}{D}, \end{equation}

where $\kappa = 2/3(\rho _p/\rho _b)(d_p/D)$ is a dimensionless ratio, accounting for density and size differences.

The experimental results of $C_1^2[1/(\kappa \ln {T(x)})^2 -1]/2$ as a function of $x/D$ are presented in figure 20, with panel (a) presenting results for a single hopper outlet thickness of $D = 5$ mm with seven different $\skew6\bar{d}_p$ and panel (b) presenting results for a single particle size $\skew6\bar{d}_p = 500\,\mathrm {\mu }$m with different $D$. The black dashed line in figure 20(a) represents a linear relationship with a gradient of unity, as expected from (4.5). It can be seen that the profile for the largest particle size, $\skew6\bar{d}_p = 500\,\mathrm {\mu }$m, matches the linear function quite closely, with the data departing from this line with a reduction in $\skew6\bar{d}_p$, particularly for $\skew6\bar{d}_p \leq 363\,\mathrm {\mu }$m. The departure of the experimental data from the trends predicted with (4.5) for sufficiently small particles, as shown in figure 20(a), is attributed to the simplifications and assumptions invoked in the derivation of (4.5), particularly in the assumption that the particle drag is zero. The relatively closely matched profile for $\skew6\bar{d}_p = 500\,\mathrm {\mu }$m can also be seen in figure 20(b) for other $D$, although some departure happens for $D = 3$ mm, potentially due to an increasing effect of non-zero drag caused by particle dispersion, which is consistent with the analysis of $K_{int}$ in figure 12(b). That is, this analytical model provides a reasonable estimate of transmittance for sufficiently large particles ($\skew6\bar{d}_p \geq 399\,\mathrm {\mu } {\rm m}$), with the hopper thickness as a secondary effect, notwithstanding the slight under-prediction in the near field and over-prediction in the intermediate field.

Figure 20. The measured particle curtain transmittance ($T$), normalised based on (4.5), as a function of normalised falling distance ($x/D$) for (a) seven different particle sizes with the hopper $D = 5$ mm and (b) four different hopper outlet thicknesses with the particle size $\skew6\bar{d}_p = 500\,\mathrm {\mu }$m. The black dashed line represents a gradient of unity, while the red dashed line indicates a different asymptotic curve fitted based on profiles of $\skew6\bar{d}_p = 163$–$257\,\mathrm {\mu }$m using a quadratic function which is included in the inset equation.

In addition, it can also be seen that this method of normalisation gives a pretty good collapse of the measurements of the curtains with smaller particles, but onto a different curve, as described by the function inserted in figure 20(a), which applies for the regime where particle drag is significant. These data have a lower slope than unity function in the near field and a higher slope for $x/D >110$. That the data collapse with this normalisation highlights the importance of initial curtain momentum in both regimes, since the normalisation utilising the Beverloo equation accounts for the particle mass flow rate and initial velocity. That the data collapse onto two curves demonstrates that curtains can be characterised broadly into the regimes in which aerodynamic drag is overall dominant or weak throughout these near and intermediate regions.

4.5. Summary of characterised flow regions for free-falling particle curtain flows

Figure 21 presents a schematic diagram of the key features of, and mechanisms for, the evolution of a free-falling particle curtain from a rectangular orifice at the base of a hopper. The ‘hopper-flow dominant regime’ identified by Sedaghatizadeh et al. (Reference Sedaghatizadeh, Arjomandi, Lau and Nathan2022) has been theoretically expanded into two regions to yield a total of four regions in the axial direction with the additional information also provided on the third region, i.e. ‘the similarity regime’, from the present new data. At the lateral plane, new information of an extra region of dispersed particles in addition to the curtain core region has also been provided. Overall, the curtain has been characterised into regions in both vertical and horizontal directions as follows:

1. (i) The near-field expansion region: this is characterised by rapid spreading immediately below the hopper exit and is deduced to be induced by particle collisions that result from the gradients in particle velocity close to the edge of the curtain due to wall friction within the hopper.

2. (ii) The neck zone: this zone is typically downstream from the near-field expansion region, and is characterised by a necking process whereby the curtain thickness decreases with axial distance. This process is attributed to the influence of particle acceleration due to gravity, which induces a low-pressure region along the axis of the curtain.

3. (iii) The intermediate-field expansion region: this region is characterised by an expansion of the curtain thickness downstream from the curtain neck (i.e. the point where the curtain thickness approaches a local minimum) due to the entrainment of air into the curtain. This region is characterised by large-scale flow unsteadiness that is deduced to be induced by a jet-like aerodynamics. The particle flow in this region transitions from the hopper-flow dominant condition to the regime where particles are approaching terminal velocity so that it theoretically matches with the similarity regime as identified by Sedaghatizadeh et al. (Reference Sedaghatizadeh, Arjomandi, Lau and Nathan2022).

4. (iv) The terminal velocity region: this ‘far-field’ region, while not investigated directly here, has been shown previously to correspond to the region in which particles decelerate towards the terminal velocity of an individual particle (Sedaghatizadeh et al. Reference Sedaghatizadeh, Arjomandi, Lau and Nathan2022).

5. (v) The curtain core: this central region of the curtain is densely loaded with particles, so that their wakes are interconnected to generate a strong particle group effect that reduces the bulk fluid drag. This region also provides a foundation for the axial regions defined above.

6. (vi) The dilute-particle region: this region laterally surrounds the curtain core, expanding continuously downstream from the near-exit region. It results from the initial particle lateral momentum, which is generated in turn by particle–wall and particle–particle interactions at the exit.

Figure 21. Schematic diagram with cartoon of key features of, and deduced mechanism for, a free-falling particle curtain discharged from a hopper with a rectangular outlet.

In the near-field expansion region, particles were found to expand immediately downstream from the exit plane with an expansion gradient of $0.0796 \leq K_{near} \leq 0.1703$, corresponding to an expansion half-angle of $2.3^{\circ } \leq \alpha \leq 4.8^{\circ }$. The near-field expansion rate increases with hopper outlet thickness and, weakly, with particle size. At the end of the near-field expansion region, a local peak in the curtain thickness was typically observed, with a local maximum curtain thickness falling in the range of $1.61 \leq \varDelta _{0.5,pk}/D \leq 2.54$. The relative axial location of the local peak, $4.96 \leq x_{pk}/D \leq 26.29$, was found to increase monotonically with particle size. Particle collisions were deduced to be the primary source of the particle lateral momentum and to be responsible for the near-field expansion, from which particle–hopper interactions can be expected to have a dominant influence on the flow behaviour. This needs to be further verified by measuring the particle velocity and the inter-particle collisions at the near-exit region.

In the neck zone, just downstream from the local maximum in the curtain thickness (i.e. at the downstream end of the near-field expansion region), the particle curtain was found to decrease in thickness, leading to a local ‘trough’ feature with the curtain thickness narrowing down to $0.78 \leq \varDelta _{0.5,tr}/D \leq 2.53$. The magnitude of $\varDelta _{0.5,tr}/D$ was found to increase with $\skew6\bar{d}_p/D$, while the difference in the magnitude between $\varDelta _{0.5, pk}/D$ and $\varDelta _{0.5, tr}/D$ was found to reduce with increasing $\skew6\bar{d}_p/D$. As a result, no obvious peak or trough feature was identified for the case of $\skew6\bar{d}_p = 500\,\mathrm {\mu }$m and $D = 3$ mm. Interestingly, the relative location of the local minimum in curtain thickness, $42.35 \leq x_{tr}/D \leq 73.72$, was found to form an ‘arch’-shaped curve with the maximum value of $x_{tr}/D$ occurring for $\skew6\bar{d}_p/D \approx 0.05$. The above findings highlight the complex nature of the flow in the neck zone, where the dominant influence is plausibly attributed to the momentum balance between the particle inertia (initial momentum) which drives the particle spreading continuously from the near-field expansion and the drag of the entrained air on the particles towards the curtain core due to the pressure gradients generated by particle acceleration. That is, particle inertia dominates the flow behaviour for the normalised particle size range of $\skew6\bar{d}_p/D > 0.05$, while the effect of surrounding fluid generated by the high particle mass flow rate is the key influence in the cases using a relatively large hopper outlet thickness, i.e. $\skew6\bar{d}_p/D < 0.05$. Further insight would require additional measurements of the particle velocity and the air velocity in the surrounding region close to the particles.

In the intermediate-field expansion region, downstream from the local minimum in curtain thickness (i.e. the neck), particles were found to expand with axial distance. In this region, large-scale flow unsteadiness that spans the entire width of the curtain was observed. The rate of expansion in the intermediate field, $0.00347 \leq K_{int} \leq 0.0226$, is approximately one order of magnitude smaller than that in the near field. Furthermore, the trend of the dependence of the expansion rate in the intermediate field on particle size and hopper outlet thickness was found to be opposite to that in the near field, with $K_{int}$ decreasing as $\skew6\bar{d}_p$ and $D$ are increased. This is consistent with the controlling mechanism being driven by aerodynamic drag, as proposed previously (Sedaghatizadeh et al. Reference Sedaghatizadeh, Arjomandi, Lau and Nathan2022). In addition, the influence of particle inertia increases with $D$ due to both the greater particle axial momentum and particle group effect reducing the fluid drag, which leads to a decrease of the $K_{int}$.

The radial particle distribution of the dilute region was found to have a rapid initial expansion comparable to that of the near-field expansion. This suggests that both are caused by similar mechanisms of particle–particle and particle–wall interactions. However, downstream from there, the dilute-particle region expands continuously, so that it does not exhibit the necking, which implies that it is not influenced by the transverse pressure gradients within the core that are induced by particle acceleration. The lateral expansion of the dilute-particle region continues to increase with a decrease in $\skew6\bar{d}_p$ through the intermediate field, which suggests the fluid drag controls its entire lateral evolution.