skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Porous cylinder arrays for optimal wake and drag characteristics
The root systems of mangroves, a tree species found in intertidal tropical and subtropical coastal zones, provide a natural barrier that dissipates wave energy effectively and reduces sediment erosion. Here, we use a combination of experiments and numerical simulations to examine the wake and drag characteristics of porous arrays of cylinders, which serve as simplified models of mangrove root networks. Optimal arrangements of the arrays are obtained by coupling Navier–Stokes simulations with a multi-objective optimization algorithm, which seeks configurations that minimize wake enstrophy and maximize drag on the porous structure. These optimal configurations are investigated using particle image velocimetry, and the internal and external flows around the porous arrays are analysed using a combination of proper orthogonal decomposition and Lagrangian particle tracking. Large variations in drag and enstrophy are observed by varying the relative positions of the cylinders, which indicates that the geometrical arrangement of porous arrays plays a prominent role in determining wake and drag characteristics. A sensitivity analysis suggests that enstrophy is more sensitive than drag to specific cylinder placement, and depends on distinctive flow patterns that develop in the interior due to interactions among neighbouring cylinders. Arrays with higher drag involve a combination of larger projected frontal area and minimal flux through the interior, leading to increased wake enstrophy, which is unfavourable for particle deposition and erosion. Based on the analysis of characteristics associated with the optimal arrays, several manually designed arrays are tested, and they display the expected behaviour with regard to drag and enstrophy.  more » « less
Award ID(s):
2103536
PAR ID:
10504044
Author(s) / Creator(s):
; ; ;
Publisher / Repository:
Cambridge University Press
Date Published:
Journal Name:
Journal of Fluid Mechanics
Volume:
961
ISSN:
0022-1120
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. null (Ed.)
    It is well known that perforation of a flat plate reduces its drag when exposed to a flow. However, studies have shown an opposite effect in the case of cylinders. Such a counterintuitive result can have significant consequences on the momentum modelling often used for wind turbine performance predictions, where increased porosity is intrinsically linked to lower drag. Here, a study of the drag of various types of porous cylinders, bars and plates under steady laminar inflow is presented. It is shown that, for most cases, the drag decreases with increased porosity. Only special types of perforations can increase the drag on both cylinders and bars, either by enhancing the effect of the rear half of the models or by organizing the wake structures. These rare occurrences are not relevant to wind turbine modelling, which indicates that current momentum models exhibit the qualitatively correct behaviour. 
    more » « less
  2. Abstract The Beaufort Gyre (BG) is hypothesized to be partially equilibrated by those mesoscale eddies that form via baroclinic instabilities of its currents. However, our understanding of the eddy field’s dependence on the mean BG currents and the role of sea ice remains incomplete. This theoretical study explores the scales and vertical structures of eddies forming specifically due to baroclinic instabilities of interior BG flows. An idealized quasi-geostrophic model is used to show that flows driven only by the Ekman pumping contain no interior potential vorticity (PV) gradients and generate weak and large eddies, ℴ(200km) in size, with predominantly barotropic and first baroclinic mode energy. However, flows containing realistic interior PV gradients in the Pacific halocline layer generate significantly smaller eddies of about 50 km in size, with a distinct second baroclinic mode structure and a subsurface kinetic energy maximum. The dramatic change in eddy characteristics is shown to be caused by the stirring of interior PV gradients by large-scale barotropic eddies. The sea ice-ocean drag is identified as the dominant eddy dissipation mechanism, leading to realistic sub-surface maxima of eddy kinetic energy for drag coefficients higher than about 2×10 −3 . A scaling law is developed for the eddy potential enstrophy, demonstrating that it is directly proportional to the interior PV gradient and the square root of the barotropic eddy kinetic energy. This study proposes a possible formation mechanism of large BG eddies and points to the importance of accurate representation of the interior PV gradients and eddy dissipation by ice-ocean drag in BG simulations and theory. 
    more » « less
  3. null (Ed.)
    Abstract Mangrove swamps are extremely productive ecosystems providing many ecological services in coastal regions. The hydrodynamic interactions of mangrove roots and water flow have been proposed as a key element to mitigate erosion. Several studies reveal that precise prediction of the morphological evolution of coastal areas, in the face of global warming and the consequent sea-level rise, requires an understanding of interactions between root porosity (the fraction of the volume of void space over the total volume), water flows, and sediment transport. Water flows around the mangrove prop roots create a complex energetic process that mixes up sediments and generates a depositional region posterior to the roots. In this work, we investigated the boundary layer behind permeable arrays of cylinders (patch) that represent the mangrove roots to explore the impact of patch porosity on the onset of sediment transport. The flow measurements were performed in a vertical plane along the water depth downstream of the mangrove root models. A high-resolution Particle Image Velocimetry (PIV) was used in a flume to observe the impact of porosity on the mean flow, velocity derivatives, skin friction coefficient, and production of turbulent kinetic energy for Reynolds number of 2500 (based on patch diameter length-scale). Here, we proposed a predictive model for critical velocity for incipient motion that takes into account the mangrove roots porosity and the near-bed turbulence effect. It is found that the patch with the $$\phi =47\%$$ ϕ = 47 % porosity, has the maximum critical velocity over which the sediment transport initiates. We found the optimum porosity has the minimum sediment erosion and creates negative vorticity sources near the bed that increases the critical velocity. This signifies an optimum porosity for the onset of sediment transport consistent with the porosity of mangroves in nature. The phenomenological model is elucidated based on an analysis of the vorticity evolution equation for viscous incompressible flows. For the optimum porous patch, a sink of vorticity was formed which yielded to lower the near-bed turbulence and vorticity. The minimum velocity fluctuations were sufficient to initiate the boundary layer transition, however, the viscous dissipation dominated the turbulence production to obstruct the sediment transport. This work identified the pivotal role of mangrove root porosity in sediment transport in terms of velocity and its derivatives in wall-bounded flows. Our work also provides insight into the sediment transport and erosion processes that govern the evolution of the shapes of shorelines. 
    more » « less
  4. The role of solid obstacle surface roughness in turbulent convection in porous media is not well understood, even though it is frequently used for heat transfer enhancement in many applications. The focus of this paper is to systematically study the influence of solid obstacle surface roughness in porous media on the microscale flow physics and report its effect on macroscale drag and Nusselt number. The Reynolds-averaged flow field is numerically simulated using the realizable k-ε model for a flow through a periodic porous medium consisting of an in-line arrangement of square cylinders with square roughness particles on the cylinder surface. Two flow regimes are identified with respect to the surface roughness particle height—fine and coarse roughness regimes. The effect of the roughness particles in the fine roughness regime is limited to the near-wall boundary layer around the solid obstacle surface. In the coarse roughness regime, the roughness particles modify the microscale flow field in the entire pore space of the porous medium. In the fine roughness regime, the heat transfer from the rough solid obstacles to the fluid inside the porous medium is less than that from a smooth solid obstacle. In the coarse roughness regime, there is an enhancement in the heat transfer from the rough solid obstacle to the fluid inside the porous medium. Total drag reduction is also observed in the fine roughness regime for the smallest roughness particle height. The surface roughness particle spacing determines the fractional area of the solid obstacle surface covered by recirculating, reattached, and stagnating flow. As the roughness particle spacing increases, there are two competing factors for the heat transfer rate—increase due to more surface area covered by reattached flow and decrease due to fewer roughness particles on the solid obstacle surface. Decreasing the porosity and increasing the Reynolds number amplify the effect of the surface roughness on the microscale flow. The results suggest that heat transfer in porous media can be enhanced, if the increase in drag can be overcome. The results also show that the fine roughness regime, which is frequently encountered due to corrosion, is detrimental to the heat transfer performance of porous media. 
    more » « less
  5. null (Ed.)
    Machine learning techniques have received attention in fluid dynamics in terms of predicting, clustering and classifying complex flow physics. One application has been the classification or clustering of various wake structures that emanate from blu˙ bodies such as cylinders or flapping foils, creating a rich diversity of vortex formations specific to flow conditions, geometry, and/or kinematics of the body. When utilizing oscillating foils to harvest energy from tidal or river flows, it is critical to understand the intricate and nonlinear relationship between flapping kinematics and the downstream vortex wake structure for optimal siting and operation of arrays. This paper develops a classification model to obtain groups of kinematics that contain similar wake patterns within the energy harvesting regime. Data is obtained through simulations of 27 unique oscillating foil kinematics for a total of 13,650 samples of the wake vorticity field. Within these samples three groups are visually labeled based on the relative angle of attack. A machine learning approach combining a convolutional neural network (CNN) with long short-term memory (LSTM) units is utilized to automatically classify the wakes into the three groups. The average accuracy on five test data subsets is 80% when the three visually labeled groups are used for classification. After analyzing the test subset with lowest accuracy, an update on the group division boundaries is proposed. With this update, the algorithm achieves an average accuracy of 90%, demonstrating that the three groups are able to discern distinct wake structures within a range of energy harvesting kinematics. 
    more » « less