A coupled finite volume immersed boundary method for simulating 3D viscoelastic flows in complex geometries

S. De, S. Das, J.A.M. Kuipers, E.A.J.F. Peters, J.T. Padding

Research output: Contribution to journalArticleAcademicpeer-review

15 Citations (Scopus)
140 Downloads (Pure)

Abstract

We report on simulations of an unsteady three dimensional viscoelastic fluid flow through a model porous medium, employing a finite volume methodology (FVM) with a staggered grid. Boundary conditions at the walls of the porous structures are imposed using a second order immersed boundary method (IBM), allowing for accurate simulations using a relatively coarse grid. We compare the viscoelastic stresses obtained using this new IBM technique with those published in literature and find good correspondence. Next, we applied this methodology to model viscoelastic fluids with a FENE-P constitutive model flowing through closely spaced cylinders. Using periodic boundary conditions, we modeled the flow behavior for Newtonian and viscoelastic fluids for successive contractions and expansions. We observe the presence of counter-rotating vortices in between the closely spaced cylinders. The viscoelastic flow structure is symmetric for lower Deborah (De) number, but onset of an asymmetry occurs after a critical De for an infinite array of cylinders. In the presence of side walls, we observe that the onset of flow asymmetry happens at a much lower De, which can be related to higher viscoelastic stresses normal to the flow direction and larger extensional viscosities which affect the curved streamlines. The three-dimensional flow characteristics for viscoelastic flow at higher De number is quite different in comparison with Newtonian flow behavior.
Original languageEnglish
Pages (from-to)67-76
Number of pages10
JournalJournal of Non-Newtonian Fluid Mechanics
Volume232
DOIs
Publication statusPublished - 13 Apr 2016

Fingerprint

Immersed Boundary Method
Viscoelastic Flow
viscoelasticity
Complex Geometry
Finite Volume Method
Viscoelastic Fluid
Boundary conditions
Newtonian flow
Fluids
Geometry
Flow structure
geometry
Constitutive models
Asymmetry
Porous materials
Flow of fluids
Vortex flow
grids
asymmetry
methodology

Keywords

  • Computational fluid dynamics
  • viscoelastic flow
  • immersed boundary method

Cite this

@article{4ea4be844f254418956dc14b1f84ebc2,
title = "A coupled finite volume immersed boundary method for simulating 3D viscoelastic flows in complex geometries",
abstract = "We report on simulations of an unsteady three dimensional viscoelastic fluid flow through a model porous medium, employing a finite volume methodology (FVM) with a staggered grid. Boundary conditions at the walls of the porous structures are imposed using a second order immersed boundary method (IBM), allowing for accurate simulations using a relatively coarse grid. We compare the viscoelastic stresses obtained using this new IBM technique with those published in literature and find good correspondence. Next, we applied this methodology to model viscoelastic fluids with a FENE-P constitutive model flowing through closely spaced cylinders. Using periodic boundary conditions, we modeled the flow behavior for Newtonian and viscoelastic fluids for successive contractions and expansions. We observe the presence of counter-rotating vortices in between the closely spaced cylinders. The viscoelastic flow structure is symmetric for lower Deborah (De) number, but onset of an asymmetry occurs after a critical De for an infinite array of cylinders. In the presence of side walls, we observe that the onset of flow asymmetry happens at a much lower De, which can be related to higher viscoelastic stresses normal to the flow direction and larger extensional viscosities which affect the curved streamlines. The three-dimensional flow characteristics for viscoelastic flow at higher De number is quite different in comparison with Newtonian flow behavior.",
keywords = "Computational fluid dynamics, viscoelastic flow, immersed boundary method",
author = "S. De and S. Das and J.A.M. Kuipers and E.A.J.F. Peters and J.T. Padding",
year = "2016",
month = "4",
day = "13",
doi = "10.1016/j.jnnfm.2016.04.002",
language = "English",
volume = "232",
pages = "67--76",
journal = "Journal of Non-Newtonian Fluid Mechanics",
issn = "0377-0257",
publisher = "Elsevier",

}

A coupled finite volume immersed boundary method for simulating 3D viscoelastic flows in complex geometries. / De, S.; Das, S.; Kuipers, J.A.M.; Peters, E.A.J.F.; Padding, J.T.

In: Journal of Non-Newtonian Fluid Mechanics, Vol. 232, 13.04.2016, p. 67-76.

Research output: Contribution to journalArticleAcademicpeer-review

TY - JOUR

T1 - A coupled finite volume immersed boundary method for simulating 3D viscoelastic flows in complex geometries

AU - De, S.

AU - Das, S.

AU - Kuipers, J.A.M.

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

AU - Padding, J.T.

PY - 2016/4/13

Y1 - 2016/4/13

N2 - We report on simulations of an unsteady three dimensional viscoelastic fluid flow through a model porous medium, employing a finite volume methodology (FVM) with a staggered grid. Boundary conditions at the walls of the porous structures are imposed using a second order immersed boundary method (IBM), allowing for accurate simulations using a relatively coarse grid. We compare the viscoelastic stresses obtained using this new IBM technique with those published in literature and find good correspondence. Next, we applied this methodology to model viscoelastic fluids with a FENE-P constitutive model flowing through closely spaced cylinders. Using periodic boundary conditions, we modeled the flow behavior for Newtonian and viscoelastic fluids for successive contractions and expansions. We observe the presence of counter-rotating vortices in between the closely spaced cylinders. The viscoelastic flow structure is symmetric for lower Deborah (De) number, but onset of an asymmetry occurs after a critical De for an infinite array of cylinders. In the presence of side walls, we observe that the onset of flow asymmetry happens at a much lower De, which can be related to higher viscoelastic stresses normal to the flow direction and larger extensional viscosities which affect the curved streamlines. The three-dimensional flow characteristics for viscoelastic flow at higher De number is quite different in comparison with Newtonian flow behavior.

AB - We report on simulations of an unsteady three dimensional viscoelastic fluid flow through a model porous medium, employing a finite volume methodology (FVM) with a staggered grid. Boundary conditions at the walls of the porous structures are imposed using a second order immersed boundary method (IBM), allowing for accurate simulations using a relatively coarse grid. We compare the viscoelastic stresses obtained using this new IBM technique with those published in literature and find good correspondence. Next, we applied this methodology to model viscoelastic fluids with a FENE-P constitutive model flowing through closely spaced cylinders. Using periodic boundary conditions, we modeled the flow behavior for Newtonian and viscoelastic fluids for successive contractions and expansions. We observe the presence of counter-rotating vortices in between the closely spaced cylinders. The viscoelastic flow structure is symmetric for lower Deborah (De) number, but onset of an asymmetry occurs after a critical De for an infinite array of cylinders. In the presence of side walls, we observe that the onset of flow asymmetry happens at a much lower De, which can be related to higher viscoelastic stresses normal to the flow direction and larger extensional viscosities which affect the curved streamlines. The three-dimensional flow characteristics for viscoelastic flow at higher De number is quite different in comparison with Newtonian flow behavior.

KW - Computational fluid dynamics

KW - viscoelastic flow

KW - immersed boundary method

U2 - 10.1016/j.jnnfm.2016.04.002

DO - 10.1016/j.jnnfm.2016.04.002

M3 - Article

VL - 232

SP - 67

EP - 76

JO - Journal of Non-Newtonian Fluid Mechanics

JF - Journal of Non-Newtonian Fluid Mechanics

SN - 0377-0257

ER -