The onset of fluid-elastic instability in cylinder arrays is usually thought to depend primarily on the mean flow velocity, the Scruton number and the natural frequency of the cylinders. Currently, there is considerable evidence from experimental measurements and computational fluid dynamic (CFD) simulations that the Reynolds number is also an important parameter. However, the available data are not sufficient to understand or quantify this effect. In this study we use a high resolution pseudo-spectral scheme to solve 2-D penalized Navier-Stokes equations in order to accurately model turbulent flow past cylinder array. To uncover the Reynolds number effect we perform simulations that vary Reynolds number independent of flow velocity at a fixed Scruton number, and then analyze the cylinder responses. The computational complexity of our algorithm is a function of Reynolds number. Therefore, we developed a high performance parallel code which allows us to simulate high Reynolds numbers at a reasonable computational cost. The simulations reveal that increasing Reynolds number has a strong de-stabilizing effect for staggered arrays. On the other hand, for the in-line array case Reynolds number still affects the instability threshold, but the effect is not monotonic with increasing Reynolds number. In addition, our findings suggest that geometry is also an important factor since at low Reynolds numbers critical flow velocity in the staggered array is considerably higher than the in-line case. This study helps to better predict how the onset of fluid-elastic instability depends on Reynolds number and reduces uncertainties in the experimental data which usually do not consider the effect of Reynolds number. / Thesis / Master of Science (MSc)

This thesis reports on an experimental study of slow viscoelastic flow through models of fibrous porous media. The models were square arrays of parallel cylinders, with solid volume fractions or ‘solidities’ of 2.5%, 5.0%, and 10%. An initial study using a Newtonian fluid provided a baseline for comparison with results for two Boger fluids, so that the effects of fluid elasticity could be determined. Boger fluids are elastic fluids that have near constant viscosities and can be used in experiments without having to account for shear-thinning effects. The experimental approach involved measurements of pressure loss through the three arrays and interior velocity measurements using particle image velocimetry (PIV). For the Newtonian flows, pressure loss measurements were in good agreement with the analytical predictions of Sangani and Acrivos (1982). PIV measurements showed velocity profiles which were symmetrical and independent of flow rate. Pressure loss measurements for the Boger fluid flows revealed that the onset of elastic effects occurred at a Deborah number of approximately 0.5, for both fluids and the three arrays. Flow resistance data collapsed for the two Boger fluids, and increased with solidity. For all three models, the flow resistance increased monotonically with Deborah number, reaching values up to four times the Newtonian resistance for the 10% model. PIV measurements showed that the transverse velocity profiles for the Newtonian and Boger fluids were the same at Deborah numbers below the elastic onset. Above onset, the profiles became skewed. The skewness, like the flow resistance, was observed to increase with both Deborah number and solidity. In the wake regions between cylinders in a column, periodic flow structures formed in the spanwise direction. The structures were staggered from column to column, consistent with the skewing. As either Deborah number or solidity increased, the flow structures became increasingly three-dimensional, and the stagger became more symmetric. An analysis of fluid stresses reveals that the elastic flow resistance is attributed to additional normal stresses caused by shearing, and not by extension.

This thesis reports on an experimental study of slow viscoelastic flow through models of fibrous porous media. The models were square arrays of parallel cylinders, with solid volume fractions or ‘solidities’ of 2.5%, 5.0%, and 10%. An initial study using a Newtonian fluid provided a baseline for comparison with results for two Boger fluids, so that the effects of fluid elasticity could be determined. Boger fluids are elastic fluids that have near constant viscosities and can be used in experiments without having to account for shear-thinning effects. The experimental approach involved measurements of pressure loss through the three arrays and interior velocity measurements using particle image velocimetry (PIV). For the Newtonian flows, pressure loss measurements were in good agreement with the analytical predictions of Sangani and Acrivos (1982). PIV measurements showed velocity profiles which were symmetrical and independent of flow rate. Pressure loss measurements for the Boger fluid flows revealed that the onset of elastic effects occurred at a Deborah number of approximately 0.5, for both fluids and the three arrays. Flow resistance data collapsed for the two Boger fluids, and increased with solidity. For all three models, the flow resistance increased monotonically with Deborah number, reaching values up to four times the Newtonian resistance for the 10% model. PIV measurements showed that the transverse velocity profiles for the Newtonian and Boger fluids were the same at Deborah numbers below the elastic onset. Above onset, the profiles became skewed. The skewness, like the flow resistance, was observed to increase with both Deborah number and solidity. In the wake regions between cylinders in a column, periodic flow structures formed in the spanwise direction. The structures were staggered from column to column, consistent with the skewing. As either Deborah number or solidity increased, the flow structures became increasingly three-dimensional, and the stagger became more symmetric. An analysis of fluid stresses reveals that the elastic flow resistance is attributed to additional normal stresses caused by shearing, and not by extension.