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
What We Can Learn From DNS of Turbulence in Porous Media: Modeling Turbulent Flow in Composite Porous/Fluid Domains
Our recently published papers reporting results of Direct Numerical Simulation (DNS) of forced convection flows in porous media suggest that in a porous medium the size of turbulent structures is restricted by the pore scale. Since the turbulent kinetic energy is predominantly contained within large eddies, this suggests that turbulent flow in a porous medium may carry less energy that its counterpart in a clear fluid domain. We use this insight to develop a practical model of turbulent flow in composite porous/fluid domains. In such domains, most of the flow is expected to occur in the clear fluid region; therefore, in most cases the flow in the porous region either remains laminar or starts its transition to turbulence even if the flow in the clear fluid region is fully turbulent. This conclusion is confirmed by comparing appropriate Reynolds numbers with their critical values. Therefore, for most cases, using the Forchheimer term in the momentum equation and the thermal dispersion term in the energy equation may result in a sufficiently good model for the porous region. However, what may really affect turbulent convection in composite domains is the roughness of the porous/fluid interface. If particles or fibers that constitute the porous medium (and the pores) are relatively large, the impact of the roughness on convection heat transfer in composite porous/fluid domains may be much more significant than the impact of possible turbulence in the porous region. We use the above considerations to develop a practical model of turbulent flow in a composite porous/fluid domain, concentrating on the effect of interface roughness on turbulence.
more »
« less
- Award ID(s):
- 1642262
- PAR ID:
- 10057659
- Date Published:
- Journal Name:
- Proceedings of the 13th International Conference on Heat Transfer, Fluid Mechanics and Thermodynamics HEFAT2017, July 17-19, 2017, Portoroz, Slovenia
- Page Range / eLocation ID:
- pp. 932-939
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
One of the most controversial topics in the field of convection in porous media is the issue of macroscopic turbulence. It remains unclear whether it can occur in porous media. It is difficult to carry out velocity measurements within porous media, as they are typically optically opaque. At the same time, it is now possible to conduct a definitive direct numerical simulation (DNS) study of this phenomenon. We examine the processes that take place in porous media at large Reynolds numbers, attempting to accurately describe them and analyze whether they can be labeled as true turbulence. In contrast to existing work on turbulence in porous media, which relies on certain turbulence models, DNS allows one to understand the phenomenon in all its complexity by directly resolving all the scales of motion. Our results suggest that the size of the pores determines the maximum size of the turbulent eddies. If the size of turbulent eddies cannot exceed the size of the pores, then turbulent phenomena in porous media differ from turbulence in clear fluids. Indeed, this size limitation must have an impact on the energy cascade, for in clear fluids the turbulent kinetic energy is predominantly contained within large eddies.more » « less
-
A ubiquitous arrangement in nature is a free-flowing fluid coupled to a porous medium, for example a river or lake lying above a porous bed. Depending on the environmental conditions, thermal convection can occur and may be confined to the clear fluid region, forming shallow convection cells, or it can penetrate into the porous medium, forming deep cells. Here, we combine three complementary approaches – linear stability analysis, fully nonlinear numerical simulations and a coarse-grained model – to determine the circumstances that lead to each configuration. The coarse-grained model yields an explicit formula for the transition between deep and shallow convection in the physically relevant limit of small Darcy number. Near the onset of convection, all three of the approaches agree, validating the predictive capability of the explicit formula. The numerical simulations extend these results into the strongly nonlinear regime, revealing novel hybrid configurations in which the flow exhibits a dynamic shift from shallow to deep convection. This hybrid shallow-to-deep convection begins with small, random initial data, progresses through a metastable shallow state and arrives at the preferred steady state of deep convection. We construct a phase diagram that incorporates information from all three approaches and depicts the regions in parameter space that give rise to each convective state.more » « less
-
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
-
A large eddy simulation (LES) is performed for a turbulent open channel flow over a porous sediment bed at permeability Reynolds number of ReK∼2.56 (Reτ = 270) representative of aquatic systems. A continuum approach based on the upscaled, volume-averaged Navier−Stokes (VaNS) equations is used by defining smoothly varying porosity across the sediment water interface (SWI) and modeling the drag force in the porous bed using a modified Ergun equation with Forchheimer corrections for inertial terms. The results from the continuum approach are compared with a pore-resolved direct numerical simulation (PR-DNS) in which turbulent flow over a randomly packed sediment bed of monodispersed particles is investigated [Karra et al.,J. Fluid Mech. 971, A23 (2023)] A spatially varying porosity profile generated from the pore-resolved DNS is used in the continuum approach. Mean flow, Reynolds stress statistics, and net momentum exchange between the freestream and the porous bed are compared between the two studies, showing reasonably good agreement. Small deviations within the transitional region between the sediment bed and the freestream as compared to the PR-DNS results are attributed to the local protrusions of particles in a randomly packed bed that are absent in the continuum approach but are present in the PR-DNS. A better representation of the effective permeability in the top transition layer that accounts for roughness effect of exposed particles is necessary. The continuum approach significantly reduces the computational cost, thereby making it suitable to study hyporheic exchange of mass and momentum in large scale aquatic domains with combined influence of bedform and bed roughness.more » « less
An official website of the United States government

