skip to main content

Explore Research Products in the PAR
It may take a few hours for recently added research products to appear in PAR search results.


Title: Symmetry breaking of turbulent flow in porous media composed of periodically arranged solid obstacles
The focus of this paper is a numerical simulation study of the flow dynamics in a periodic porous medium to analyse the physics of a symmetry-breaking phenomenon, which causes a deviation in the direction of the macroscale flow from that of the applied pressure gradient. The phenomenon is prominent in the range of porosity from 0.43 to 0.72 for circular solid obstacles. It is the result of the flow instabilities formed when the surface forces on the solid obstacles compete with the inertial force of the fluid flow in the turbulent regime. We report the origin and mechanism of the symmetry-breaking phenomenon in periodic porous media. Large-eddy simulation (LES) is used to simulate turbulent flow in a homogeneous porous medium consisting of a periodic, square lattice arrangement of cylindrical solid obstacles. Direct numerical simulation is used to simulate the transient stages during symmetry breakdown and also to validate the LES method. Quantitative and qualitative observations are made from the following approaches: (1) macroscale momentum budget and (2) two- and three-dimensional flow visualization. The phenomenon draws its roots from the amplification of a flow instability that emerges from the vortex shedding process. The symmetry-breaking phenomenon is a pitchfork bifurcation that can exhibit multiple modes depending on the local vortex shedding process. The phenomenon is observed to be sensitive to the porosity, solid obstacle shape and Reynolds number. It is a source of macroscale turbulence anisotropy in porous media for symmetric solid-obstacle geometries. In the macroscale, the principal axis of the Reynolds stress tensor is not aligned with any of the geometric axes of symmetry, nor with the direction of flow. Thus, symmetry breaking in porous media involves unresolved flow physics that should be taken into consideration while modelling flow inhomogeneity in the macroscale.  more » « less
Award ID(s):
2042834
PAR ID:
10299875
Author(s) / Creator(s):
; ; ;
Date Published:
Journal Name:
Journal of Fluid Mechanics
Volume:
929
ISSN:
0022-1120
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ; ; ; ( , Proceeding of 5-6th Thermal and Fluids Engineering Conference (TFEC))
    null (Ed.)
    Turbulence modeling in porous media can be greatly improved by combining high-resolution numerical methods with modern data-driven techniques. The development of accurate macroscale models (length scale greater than the pore size) will enable real-time systemic simulations of porous media flow. We consider the case of turbulent flow in homogeneous porous media, typically encountered in engineered porous media (heat exchangers, metamaterials, combustors, etc.). The underlying microscale flow field is inhomogeneous and determined by the geometry of the porous medium. Neural Networks are able to resolve the geometry-dependence and the non-linearity of porous media turbulent flow. We are proposing to separate the macroscale model into individual blocks that predict a unique aspect of the microscale flow, such as microscale spatial flow distribution and vortex dynamics. In the present work, we determine the feasibility of the prediction of the Reynolds-averaged microscale flow patterns by using Convolutional Neural Networks (CNN). The porous medium is represented by using a square lattice arrangement of circular cylinder solid obstacles. The pore-scale Reynolds number of the flow is 300. The porosity of the porous medium is varied from 0.45 to 0.92 with 60 steps. The microscale flow field is simulated by using Large Eddy Simulation (LES) with a compact sixth-order finite difference method. We demonstrate satisfactory prediction of the microscale flow field using the CNN with a global error less than 10%. We vary the number of training samples to study the deterioration of the model accuracy. The CNN model offers a O(106) speedup over LES with only 10% loss in accuracy. 
    more » « less
  2. ; ; ( , Proceeding of 5-6th Thermal and Fluids Engineering Conference (TFEC))
    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
  3. ; ; ( , Transport in Porous Media)
    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
  4. ; ( , Physics of Fluids)
    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
  5. ; ( , AIAA Journal)
    A novel formulation of the vortex particle method (VPM) is developed for large-eddy simulation (LES) in a meshless scheme that is numerically stable. A new set of VPM governing equations are derived from the LES-filtered Navier–Stokes equations. The new equations reinforce the conservation of angular momentum by resizing vortex elements subject to vortex stretching. In addition to the VPM reformulation, a new anisotropic dynamic model of subfilter-scale (SFS) vortex stretching is developed. This SFS model is well suited for turbulent flows with coherent vortical structures, where the predominant cascade mechanism is vortex stretching. The mean and fluctuating components of turbulent flow and Reynolds stresses are validated through the simulation of a turbulent round jet. The computational efficiency of the scheme is showcased in the simulation of an aircraft rotor in hover, showing our meshless LES to be 100 times faster than a mesh-based LES with similar fidelity. The implementation of our meshless LES scheme is released as open-source software, called FLOWVPM 
    more » « less
  • Citation Formats
  • MLA
  • APA
  • Chicago
  • BibTeX
  • Save / Share this Record
  • Send to Email