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
Using Direct Numerical Simulations for Investigating Physics of Turbulence in Porous Media
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
- Award ID(s):
- 1642262
- PAR ID:
- 10057667
- Date Published:
- Journal Name:
- Proceedings of the ASME 2017 Fluids Engineering Division Summer Meeting, FEDSM2017, July 31-August 3, 2017, Waikoloa, Hawai'i, USA
- Page Range / eLocation ID:
- FEDSM2017-69041
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Natural convection in porous media is a fundamental process for the long-term storage of CO 2 in deep saline aquifers. Typically, details of mass transfer in porous media are inferred from the numerical solution of the volume-averaged Darcy–Oberbeck–Boussinesq (DOB) equations, even though these equations do not account for the microscopic properties of a porous medium. According to the DOB equations, natural convection in a porous medium is uniquely determined by the Rayleigh number. However, in contrast with experiments, DOB simulations yield a linear scaling of the Sherwood number with the Rayleigh number ( $Ra$ ) for high values of $Ra$ ( $Ra\gg 1300$ ). Here, we perform direct numerical simulations (DNS), fully resolving the flow field within the pores. We show that the boundary layer thickness is determined by the pore size instead of the Rayleigh number, as previously assumed. The mega- and proto-plume sizes increase with the pore size. Our DNS results exhibit a nonlinear scaling of the Sherwood number at high porosity, and for the same Rayleigh number, higher Sherwood numbers are predicted by DNS at lower porosities. It can be concluded that the scaling of the Sherwood number depends on the porosity and the pore-scale parameters, which is consistent with experimental studies.more » « less
-
null (Ed.)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
-
Recent numerical predictions of turbulent boundary layers subject to very strong Favorable Pressure Gradient (FPG) with high spatial/temporal resolution, i.e. Direct Numerical Simulation (DNS), have shown a meaningful weakening of the Reynolds shear stresses with a lengthy logarithmic behavior [1,2]. In the present study, assessment of the Shear Stress Transport and Spalart-Allmaras turbulence models (hence- forth SST and SA, respectively) in Reynolds-averaged Navier-Stokes (RANS) simulations is performed. The main objective is to evaluate the ability of popular turbulence models in capturing the characteristic features present during the quasi-laminarization phenomenon in highly accelerating turbulent boundary layers. A favorable pressure gradient is prescribed by a top converging surface (sink flow) with an approximately constant acceleration parameter of K = 4 . 0 ×10 −6 . Validation of RANS results is carried out by means of a large DNS dataset [1]. Generally speaking, the SA turbulence model has demonstrated the best compromise between accuracy and quick adaptation to the turbulent inflow conditions. Turbulence models properly captured the increasing trend of the freestream and friction velocity in highly accelerated flows; however, they fail to reproduce the decreasing behavior of the skin friction coefficient, which is typical in early stages of the quasi-laminarization process. Both models have shown deficient predictions of the decreasing and logarithmic behavior of Reynolds shear stresses as well as significantly overpredicted the production of Turbulent Kinetic Energy (TKE) in turbulent boundary layers subject to very strong FPG. https://doi.org/10.1016/j.compfluid.2020.104494more » « less
-
more » « less
Abstract In very stable boundary layers (VSBL), a “cocktail” of submeso motions routinely result in elevated mean wind speed maxima above the ground, acting as a new source of turbulence generation. This new source of turbulent kinetic energy enhances turbulent mixing and causes mean wind profile distortion (WPD). As a results, this transient distortion in the wind profile adjusts the classical log‐law. Addressing how WPD‐induced turbulence regulates flow structures, turbulent fluxes, and transitions in stability regimes across layers remains a challenge. Eddy covariance data measured at four levels on a 62‐m tower are employed to address these questions. It is shown that the WPD initiates large turbulent eddies that penetrate downward, leading to enhanced vertical mixing and comparable turbulent transport efficiencies across layers. As a consequence, turbulence intensity and fluxes are increased. As the WPD is intensified, turbulent fluxes and turbulent flux transport caused by large eddies are also enhanced, leading to a transition from very stable to weakly stable regimes. Due to the influence of WPD‐induced large eddies, the large‐eddy turbulent Prandtl number does not deviate appreciably from unity and the partitioning between turbulent kinetic and potential energies is linearly related to the gradient Richardson number.
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Conference Paper:
- The DOI is not currently available.
