Numerical Simulation of Turbulent Flow and Heat Transfer in a Three-Dimensional Channel Coupled with Flow Through Porous Structures

Transport of Turbulence Across Permeable Interface in a Turbulent Channel Flow: Interface-Resolved Direct Numerical Simulation

Turbulence transportation across permeable interfaces is investigated using direct numerical simulation, and the connection between the turbulent surface flow and the pore flow is explored. The porous media domain is constructed with an in-line arranged circular cylinder array. The effects of Reynolds number and porosity are also investigated by comparing cases with two Reynolds numbers ($$Re\approx 3000,6000$$ R e ≈ 3000 , 6000 ) and two porosities ($$\varphi =0.5,0.8$$ φ = 0.5 , 0.8 ). It was found that the change of porosity leads to the variation of flow motions near the interface region, which further affect turbulence transportation below the interface. The turbulent kinetic energy (TKE) budget shows that turbulent diffusion and pressure transportation work as energy sink and source alternatively, which suggests a possible route for turbulence transferring into porous region. Further analysis on the spectral TKE budget reveals the role of modes of different wavelengths. A major finding is that mean convection not only affects the distribution of TKE in spatial space, but also in scale space. The permeability of the wall also have an major impact on the occurrence ratio between blow and suction events as well as their corresponding flow structures, which can be related to the change of the Kármán constant of the mean velocity profile.

A Hybrid-Dimensional Coupled Pore-Network/Free-Flow Model Including Pore-Scale Slip and Its Application to a Micromodel Experiment

Modeling coupled systems of free flow adjacent to a porous medium by means of fully resolved Navier–Stokes equations is limited by the immense computational cost and is thus only feasible for relatively small domains. Coupled, hybrid-dimensional models can be much more efficient by simplifying the porous domain, e.g., in terms of a pore-network model. In this work, we present a coupled pore-network/free-flow model taking into account pore-scale slip at the local interfaces between free flow and the pores. We consider two-dimensional and three-dimensional setups and show that our proposed slip condition can significantly increase the coupled model’s accuracy: compared to fully resolved equidimensional numerical reference solutions, the normalized errors for velocity are reduced by a factor of more than five, depending on the flow configuration. A pore-scale slip parameter $$\beta _{{{{\rm pore}}}}$$ β pore required by the slip condition was determined numerically in a preprocessing step. We found a linear scaling behavior of $$\beta _{{{{\rm pore}}}}$$ β pore with the size of the interface pore body for three-dimensional and two-dimensional domains. The slip condition can thus be applied without incurring any run-time cost. In the last section of this work, we used the coupled model to recalculate a microfluidic experiment where we additionally exploited the flat structure of the micromodel which permits the use of a quasi-3D free-flow model. The extended coupled model is accurate and efficient.

Obstacles, Interfacial Forms, and Turbulence: A Numerical Analysis of Soil–Water Evaporation Across Different Interfaces

Exchange processes between a turbulent free flow and a porous media flow are sensitive to the flow dynamics in both flow regimes, as well as to the interface that separates them. Resolving these complex exchange processes across irregular interfaces is key in understanding many natural and engineered systems. With soil–water evaporation as the natural application of interest, the coupled behavior and exchange between flow regimes are investigated numerically, considering a turbulent free flow as well as interfacial forms and obstacles. Interfacial forms and obstacles will alter the flow conditions at the interface, creating flow structures that either enhance or reduce exchange rates based on their velocity conditions and their mixing with the main flow. To evaluate how these interfacial forms change the exchange rates, interfacial conditions are isolated and investigated numerically. First, different flow speeds are compared for a flat surface. Second, a porous obstacle of varied height is introduced at the interface, and the effects the flow structures that develop have on the interface are analyzed. The flow parameters of this obstacle are then varied and the interfacial exchange rates investigated. Next, to evaluate the interaction of flow structures between obstacles, a second obstacle is introduced, separated by a varied distance. Finally, the shape of these obstacles is modified to create different wave forms. Each of these interfacial forms and obstacles is shown to create different flow structures adjacent to the surface which alter the mass, momentum, and energy conditions at the interface. These changes will enhance the exchange rate in locations where higher velocity gradients and more mixing with the main flow develop, but will reduce the exchange rate in locations where low velocity gradients and limited mixing with the main flow occur.

Stochastic Analysis of the Gas Flow at the Gas Diffusion Layer/Channel Interface of a High-Temperature Polymer Electrolyte Fuel Cell

Gas diffusion layers (GDLs) play a significant role in the efficient operation of high-temperature polymer electrolyte fuel cells. They connect the electrodes to the gas channels of the bipolar plate by porous material with a meso-scale geometric structure. The electrodes must be sufficiently supplied by gases from the channels to operate fuel cells efficiently. Furthermore, reaction products must be transported in the other direction. The gas transport is simulated in the through-plane direction of the GDL, and its microstructure created by a stochastic model is equivalent to the structure of real GDL material. Continuum approaches in cell-scale simulations have model parameters for porous regions that can be taken from effective properties calculated from the meso-scale simulation results, as one feature of multi-scale simulations. Another significant issue in multi-scale simulations is the interface between two regions. The focus is on the gas flow at the interface between GDL and the gas channel, which is analyzed using statistical methods. Quantitative relationships between functionality and microstructure can be detected. With this approach, virtual GDL materials can possibly be designed with improved transport properties. The evaluation of the surface flow with stochastic methods offers substantiated benefits that are suitable for connecting the meso-scale to larger spatial scales.