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.
Turbulent flow in a homogeneous porous medium was investigated through the use of numerical methods by employing the Reynolds Averaged Navier-Stokes (RANS) modeling technique. The focus of our research was to study how microscopic vortices in porous media flow influence the heat transfer from the solid obstacles comprising the porous medium to the fluid. A Representative Elementary Volume (REV) with 4 × 4 cylindrical obstacles and periodic boundary conditions was used to represent the infinite porous medium structure. Our hypothesis is that the rate of heat transfer between the obstacle surface and the fluid (qavg) is strongly influenced by the size of the contact area between the vortices and the solid obstacles in the porous medium (Avc). This is because vortices are regions with low velocity that form an insulating layer on the surface of the obstacles. Factors such as the porosity (φ), Pore Scale Reynolds number (Rep), and obstacle shape of the porous medium were investigated. All three of these factors have different influences on the contact area Avc, and, by extension, the overall heat transfer rate qavg. Under the same Pore Scale Reynolds number (Rep), our results suggest that a higher overall heat transfer rate is exhibited for smaller contact areas between the vortices and the obstacle surface. Although the size of the contact area, Avc, is affected by Rep, the direct influence of Rep on the overall heat transfer rate qavg is much stronger, and exceeds the effect of Avc on qavg. The Pore Scale Reynolds number, Rep, and the mean Nusselt number, Num, have a seemingly logarithmic relationship.
New insight into the contribution of the microscale vortex evolution to convection heat transfer in porous media is presented in this paper. The objective is to determine how the microscale vortices influence convection heat transfer in turbulent flow inside porous media. The microscale temperature distribution is analysed using flow visualization in two dimensions using streamlines and in three dimensions using the Q-criterion. The pertinent observations are supplemented with a comparison of surface skin friction and heat transfer using: (i) surface skin-friction lines and (ii) the joint probability density function of the pressure and skin-friction coefficients, along with the Nusselt number. The microscale flow phenomena observed are corroborated with the features of the frequency spectra of the drag coefficient and macroscale Nusselt number. The large eddy simulation technique is used in this study to investigate the flow field inside a periodic porous medium. The Reynolds numbers of the flow are 300 and 500. The porous medium consists of solid obstacles in the shape of square and circular cylinders. Two distinct flow regimes are represented by using the porosities of 0.50 and 0.87. The results show that the surface Nusselt number distribution is dependent on whether the micro-vortices are attached to or detached from the surface of the obstacle. The spectra of the macroscale Nusselt number and the pressure drag are similar, signifying a correlation between the dynamics of heat transfer and the microscale turbulent structures. Both vortex shedding and secondary flow instabilities are observed that significantly influence the Nusselt number. The fundamental insight gained in this paper can inform the development of more robust macroscale models of convection heat transfer in turbulent flow in porous media.
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.
hi@scite.ai
10624 S. Eastern Ave., Ste. A-614
Henderson, NV 89052, USA
Copyright © 2024 scite LLC. All rights reserved.
Made with 💙 for researchers
Part of the Research Solutions Family.