Improved PISO algorithms for modeling density varying flow in conjugate fluid–porous domains

Abstract

Two modified segregated PISO algorithms are proposed, which are constructed to avoid the development of spurious oscillations in the computed flow near sharp interfaces of conjugate fluid-porous domains. The new collocated finite volume algorithms modify the Rhie-Chow interpolation to maintain a correct pressure-velocity coupling when large discontinuous momentum sources associated with jumps in the local permeability and porosity are present. The Re-Distributed Resistivity (RDR) algorithm is based on spreading flow resistivity over the grid cells neighboring a discontinuity in material properties of the porous medium. The Face Consistent Pressure (FCP) approach derives an auxiliary pressure value at the fluid-porous interface using momentum balance around the interface. Such derived pressure correction is designed to avoid spurious oscillations as would otherwise arise with a strictly central discretization. The proposed algorithms are successfully compared against published data for the velocity and pressure for two reference cases of viscous flow. The robustness of the proposed algorithms is additionally demonstrated for strongly reduced viscosity, i.e., higher Reynolds number flows and low Darcy numbers, i.e., low permeability of the porous regions in the domain, for which solutions without unphysical oscillations are computed. Both RDR and FCP are found to accurately represent porous media flow near discontinuities in material properties on structured grids.

Polymer Flow Through Porous Media: Numerical Prediction of the Contribution of Slip to the Apparent Viscosity

“…Since the grids are collocated (pressure and velocity are calculated at the same point), a Rhie-Chow type interpolation (Rhie and Chow 1983) is employed during the pressure correction steps in OpenFOAM, by estimating the pressure gradient on the cell faces with the pressure values of the neighboring cell centers. This Rhie-Chow method prevents the solution from becoming unstable, by enforcing the pressure/velocity coupling in the course of the SIMPLE algorithm [see a detailed description in Nordlund et al (2016)]. …”

Analysis and Prediction of Methane Invasion Distance Considering Real Ground Boundary

“…The finite volume method of ANSYS FLUENT software in pressure-based solver mode was used for the study. For the coupling method of velocity and pressure in simulation calculation, simple algorithm was used for steady simulation calculation and PISO algorithm was used for transient simulation calculation. , To ensure the accuracy of calculation, the second-order upwind discretization format was adopted. For steady simulation, the calculation residual was less than 10 –6 and the calculation was stopped to ensure the absolute convergence.…”

Neams Th Crab

“…Other modifications to the Rhie Chow interpolation, proposed by Zhang [18] and Nordlund [19], may also be considered. One proposed modification is to use a baffle/interface approach to connect two subdomain solutions with a simple equation for a set pressure drop, as laid out by [20].…”