TY - JOUR

T1 - Viscoelastic flow simulations in random porous media

AU - De, S.

AU - Kuipers, J. A. M.

AU - Peters, E. A. J. F.

AU - Padding, J. T.

PY - 2017

Y1 - 2017

N2 - We investigate creeping flow of a viscoelastic fluid through a three dimensional random porous medium using computational fluid dynamics. The simulations are performed using a finite volume methodology with a staggered grid. The no slip boundary condition on the fluid-solid interface is implemented using a second order finite volume immersed boundary (FVM-IBM) methodology [1]. The viscoelastic fluid is modeled using a FENE-P type model. The simulations reveal a transition from a laminar regime to a nonstationary regime with increasing viscoelasticity. We find an increased flow resistance with increase in Deborah number even though shear rheology is shear thinning nature of the fluid. By choosing a length scale based on the permeability of the porous media, a Deborah number can be defined, such that a universal curve for the flow transition is obtained. A study of the flow topology shows how in such disordered porous media shear, extensional and rotational contributions to the flow evolve with increased viscoelasticity. We correlate the flow topology with the dissipation function distribution across the porous domain, and find that most of the mechanical energy is dissipated in shear dominated regimes instead, even at high viscoelasticity.

AB - We investigate creeping flow of a viscoelastic fluid through a three dimensional random porous medium using computational fluid dynamics. The simulations are performed using a finite volume methodology with a staggered grid. The no slip boundary condition on the fluid-solid interface is implemented using a second order finite volume immersed boundary (FVM-IBM) methodology [1]. The viscoelastic fluid is modeled using a FENE-P type model. The simulations reveal a transition from a laminar regime to a nonstationary regime with increasing viscoelasticity. We find an increased flow resistance with increase in Deborah number even though shear rheology is shear thinning nature of the fluid. By choosing a length scale based on the permeability of the porous media, a Deborah number can be defined, such that a universal curve for the flow transition is obtained. A study of the flow topology shows how in such disordered porous media shear, extensional and rotational contributions to the flow evolve with increased viscoelasticity. We correlate the flow topology with the dissipation function distribution across the porous domain, and find that most of the mechanical energy is dissipated in shear dominated regimes instead, even at high viscoelasticity.

UR - http://resolver.tudelft.nl/uuid:f9e8c39b-e4e9-4915-8c15-58dd5d2c642d

U2 - 10.1016/j.jnnfm.2017.08.010

DO - 10.1016/j.jnnfm.2017.08.010

M3 - Article

VL - 248

SP - 50

EP - 61

JO - Journal of Non-Newtonian Fluid Mechanics

JF - Journal of Non-Newtonian Fluid Mechanics

SN - 0377-0257

ER -