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 -