Smoothed Particle Hydrodynamics Simulations of Porous Medium Flow Using Ergun’s Fixed-Bed Equation

• Published in: Water (Basel) Vol. 15; no. 13; p. 2358
• Main Authors: Alvarado-Rodríguez, Carlos E., Díaz-Damacillo, Lamberto, Plaza, Eric, Sigalotti, Leonardo Di G.
• Format: Journal Article
• Language: English
• Published: Basel MDPI AG 01-07-2023
• Subjects: Accuracy; Arrays; Comparative analysis; Enhanced oil recovery; Ergun equation; Fluid dynamics; Fluid flow; Fluid mechanics; Fluidized beds; Hydrodynamics; Methods; numerical methods; packed beds; Permeability; pore-scale simulation; Porosity; Porous materials; Porous media; Pressure loss; REV-scale simulation; Reynolds number; Simulation; Simulation methods; Smooth particle hydrodynamics; Velocity; Velocity distribution
• ISSN: 2073-4441; 2073-4441
• DOI: 10.3390/w15132358

A popular equation that is often employed to represent the relationship between the pressure loss and the fluid flow in fluidized or packed granular beds is the Ergun equation, which is an extension of Darcy’s law. In this paper, the method of Smoothed Particle Hydrodynamics (SPH) is used to numerically study the flow field across a rectangular channel partially filled with a porous layer both at the Representative Elementary Volume (REV) scale using the Ergun equation and at the pore scale. Since the flow field can be estimated at the REV scale with a much lower cost compared to the pore scale, it is important to evaluate how accurately the pore-scale results can be reproduced at the REV scale. The comparison between both scales is made in terms of the velocity profiles at the outlet of the rectangular channel and the pressure losses across the clear and porous zones for three different arrays of solid grains at the pore scale. The results show that minimum differences in the flow structure and velocity profiles between the REV and the pore scale always occur at intermediate values of the porosity (ϕ=0.44 and 0.55). As the porosity increases, the differences between the REV and the pore scale also increase. The details of the pressure losses are affected by the geometry of the porous medium. In particular, we find that the pressure profiles at the REV scale match those at the pore scale almost independently of the porosity only when the grains are uniformly distributed in a non-staggered square array.