skip to main content


Title: Volumetric measurement of a Newtonian fluid flow through three-dimensional porous media using Lagrangian particle tracking (Shake-the-Box) technique
This work experimentally investigates the pressure-driven flow of a pure Newtonian fluid through three-dimensional (3D) porous media models. The porous media model consists of square arrays of rods that also could be interpreted as a periodic tandem rod arrangement. We employed a time-resolved three-dimensional particle tracking velocimetry (3D Shake-the-Box) technique for a range of Reynolds numbers [Formula: see text] to observe flow structures and vortex formation between the rods in porous media structures with different porosities of [Formula: see text] which corresponds to the spacing ratio of [Formula: see text], where L is the distance between the centers of the rods, and D is the diameter of the rods. For all the examined cases, we further analyzed the effect of the Reynolds number and the spacing ratio on the instantaneous and averaged patterns of velocity, vorticity, and the other flow parameters after obtaining the two-dimensional velocity fields using the bin-averaging method. We observed both symmetrical and asymmetrical patterns of structure and recirculation regions between the rods depending on the Reynolds number and spacing ratio. Increasing the Reynolds number reduced the symmetrical patterns of flow structures with respect to the centerline of the gap region, while the spacing ratio was randomly affecting the symmetry degree. Vortex shedding was considerable for the two examined high Reynolds numbers of Re = 444 and Re = 890 behind the upstream rod as the porosity increased. The backward movement of the reattachment point has been observed by increasing the Reynolds number.  more » « less
Award ID(s):
1854376
PAR ID:
10430729
Author(s) / Creator(s):
;
Date Published:
Journal Name:
Physics of Fluids
Volume:
35
Issue:
4
ISSN:
1070-6631
Page Range / eLocation ID:
043605
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. High-resolution large eddy simulations and complementary laboratory experiments using particle image velocimetry were performed to provide a detailed quantitative assessment of flow response to gaps in cylinder arrays. The base canopy consists of a dense array of emergent rigid cylinders placed in a regular staggered pattern. The gaps varied in length from [Formula: see text] to 24, in intervals of 4 d, where d is the diameter of the cylinders. The analysis was performed under subcritical conditions with Froude numbers [Formula: see text] and bulk Reynolds numbers [Formula: see text]. Results show that the gaps affect the flow statistics at the upstream and downstream proximity of the canopy. The affected zone was [Formula: see text] for the mean flow and [Formula: see text] for the second-order statistics. Dimensionless time-averaged streamwise velocity within the gap exhibited minor variability with gap spacing; however, in-plane turbulent kinetic energy, k, showed a consistent decay rate when normalized with that at [Formula: see text] from the beginning of the gap. The emergent canopy acts as a passive turbulence generator for the gap flow for practical purposes. The streamwise dependence of k follows an exponential trend within [Formula: see text] and transitions to a power-law at [Formula: see text]. The substantially lower maximum values of k within the gap compared to k within the canopy evidence a limitation of gap measurements representative of canopy flow statistics. We present a base framework for estimating representative in-canopy statistics from measurements in the gap. 
    more » « less
  2. In this paper, we report a novel symmetry-breaking phenomenon that occurs in turbulent convective flow in periodic porous media in the intermediate porosity flow regime for values of porosities in between 0.8 and 0.9. Large eddy simulation is used to numerically simulate the momentum and thermal transport inside the porous medium at the microscale level. The phenomenon is observed to occur for periodically repeating porous media consisting of an in-line arrangement of circular cylinder solid obstacles, such as typically found in heat exchangers. The transition from symmetric to asymmetric flow occurs in between the pore scale Reynolds numbers of 37 (laminar) and 100 (turbulent), and asymmetric flow patterns are reported for Reynolds numbers up to 1,000. A Hopf bifurcation resulting in unsteady oscillatory laminar flow marks the origin of a secondary flow instability arising from the interaction of the shear layers around the solid obstacle. In turbulent flow, stochastic phase difference in the vortex wake oscillations caused by the secondary flow instability results in asymmetrical velocity and temperature distributions in the pore space. Consequently, high and low velocity flow channels are formed in the pore space that leads to the asymmetrical velocity and pressure distributions. At the macroscale level, symmetry-breaking results in residual transverse drag force components acting on the solid obstacle surfaces. The vortex wake oscillations caused by the secondary flow instability promote attached flow on the solid obstacle surface, which improves the surface averaged heat flux from the solid obstacle surface. 
    more » « less
  3. This paper describes the numerical study of oscillating circular cylinders with rigid splitter plates of different lengths. These geometries may be used as disturbance generators for the study of unsteady airfoils and wings operating in highly vortical flowfields. It has been shown that cylinders undergoing forced rotational oscillations at their natural shedding frequency can produce wakes with minimal deviation in cycle-to-cycle vortex strength and position. Adding a splitter plate allows these deviations to be reduced even further. We present cases for oscillating cylinders having splitter-plate lengths up to [Formula: see text] at a Reynolds number of 7600. Frequencies are maintained at the natural shedding frequency, and a rotational amplitude of 45 deg is used. Numerical simulations are performed using a two-dimensional unsteady Reynolds-averaged Navier–Stokes (RANS) code. Results are presented in the form of vorticity contours and cycle-averaged velocity profiles, as well as the dominant frequencies of cylinder lift force and downstream velocity angles. The results show that splitter-plate lengths shorter than [Formula: see text] adversely affect the ability to generate a coherent vortex wake due to shear layer roll-up near the trailing edge of the plate. Splitter plates longer than [Formula: see text] produced a reverse von Kármán wake with consistent cycle-to-cycle vortex shedding.

     
    more » « less
  4. When fluids flow through straight channels sustained turbulence occurs only at high Reynolds numbers [typically [Formula: see text]]. It is difficult to mix multiple fluids flowing through a straight channel in the low Reynolds number laminar regime [[Formula: see text]] because in the absence of turbulence, mixing between the component fluids occurs primarily via the slow molecular diffusion process. This Letter reports a simple way to significantly enhance the low Reynolds number (in our case [Formula: see text]) passive microfluidic flow mixing in a straight microchannel by introducing asymmetric wetting boundary conditions on the floor of the channel. We show experimentally and numerically that by creating carefully chosen two-dimensional hydrophobic slip patterns on the floor of the channels, we can introduce stretching, folding, and/or recirculation in the flowing fluid volume, the essential elements to achieve mixing in the absence of turbulence. We also show that there are two distinctive pathways to produce homogeneous mixing in microchannels induced by the inhomogeneity of the boundary conditions. It can be achieved either by (1) introducing stretching, folding and twisting of fluid volumes, i.e., via a horse-shoe type transformation map, or (2) by creating chaotic advection, achieved through manipulation of the hydrophobic boundary patterns on the floor of the channels. We have also shown that by superposing stretching and folding with chaotic advection, mixing can be optimized in terms of significantly reducing mixing length, thereby opening up new design opportunities for simple yet efficient passive microfluidic reactors. 
    more » « less
  5. null (Ed.)
    The influence of microscale flow structures (smaller than the pore size) on turbulent heat transfer in porous media has not been yet investigated. The goal of this study is to determine the influence of the micro-vortices on convection heat transfer in turbulent porous media flow. Turbulent flow in a homogeneous porous medium was investigated using Large Eddy Simulation (LES) at a Reynolds number of 300. We observed that the convection heat transfer characteristics are dependent on whether the micro-vortices are attached or detached from the surface of the obstacle. There is a spectral correlation between the Nusselt number and the pressure instabilities due to vortex shedding. A secondary flow instability occurs due to high pressure regions forming periodically near the converging pathway between obstacles. This causes local adverse pressure gradient, affecting the flow velocity and convection heat transfer. This study has been performed for obstacles with shapes of square and circular cylinders at porosities of 0.50 and 0.87. Understanding the dominant modes that affect convection heat transfer can aid in finding an optimum geometry for the porous medium. 
    more » « less