An immersed boundary method for computing heat and fluid flow in porous media

David J. Lopez Penha; Lilya Ghazaryan; Bernard J. Geurts; Steffen Stolz; Markus Nordlund

An immersed boundary method for computing heat and fluid flow in porous media

David J. Lopez Penha, Lilya Ghazaryan, Bernard J. Geurts, Steffen Stolz, Markus Nordlund

Research output: Chapter in Book/Report/Conference proceeding › Conference contribution › Academic › peer-review

94 Downloads (Pure)

Abstract

A volume-penalizing immersed boundary (IB) method is presented that facilitates the computation of fluid flow in complex porous media. The computational domain is composed of a uniform Cartesian grid, and solid bodies are approximated on this grid using a series of grid cells (i.e., a ''staircase'' approximation). Solid bodies are distinguished from fluid regions using a binary phase-indicator function: Taking the value of ''1'' in the solid parts of the domain and ''0'' in the fluid parts. The effect of solid bodies on the flow is modeled using a source term in the momentum equations. The source term is active only within solid parts of the domain, and enforces the no-slip boundary condition. Fluid regions are governed by the incompressible Navier-Stokes equations. An extension of the IB method is proposed to tackle coupled fluid-solid heat transfer. The extended IB method is validated for Poiseuille flow, which allows for a direct comparison of the numerical results against a closed analytical solution. We subsequently apply the extended IB method to flow in a structured porous medium and focus on bulk properties such as the gradient of the average pressure and the Nusselt number. Reliable qualitative results were obtained with 16-32 grid points per singly-connected fluid region.

Original language	English
Title of host publication	Proceedings of the V European Conference on Computational Fluid Dynamics, ECCOMAS CFD 2010
Editors	J.C.F Pereira, A. Sequeira, J.M.C. Pereira
Place of Publication	Lisbon
Publisher	Technical University of Lisbon
Number of pages	16
ISBN (Print)	978-989-96778-1-4
Publication status	Published - Jun 2010
Event	5th European Conference on Computational Fluid Dynamics, ECCOMAS ECFD 2010 - Lisbon, Portugal Duration: 14 Jun 2010 → 17 Jun 2010 Conference number: 5

Conference

Conference	5th European Conference on Computational Fluid Dynamics, ECCOMAS ECFD 2010
Abbreviated title	ECCOMAS ECFD 2010
Country/Territory	Portugal
City	Lisbon
Period	14/06/10 → 17/06/10

Keywords

Incompressible Navier-Stokes equations
Immersed boundary method
Numerical simulation
Porous media
Volume penalization
Conjugate heat transfer

Access to Document

ArticleFinal published version, 685 KB

Cite this

@inproceedings{3cb50cc4481e4d7ebf0136cc40ea14d2,

title = "An immersed boundary method for computing heat and fluid flow in porous media",

abstract = "A volume-penalizing immersed boundary (IB) method is presented that facilitates the computation of fluid flow in complex porous media. The computational domain is composed of a uniform Cartesian grid, and solid bodies are approximated on this grid using a series of grid cells (i.e., a ''staircase'' approximation). Solid bodies are distinguished from fluid regions using a binary phase-indicator function: Taking the value of ''1'' in the solid parts of the domain and ''0'' in the fluid parts. The effect of solid bodies on the flow is modeled using a source term in the momentum equations. The source term is active only within solid parts of the domain, and enforces the no-slip boundary condition. Fluid regions are governed by the incompressible Navier-Stokes equations. An extension of the IB method is proposed to tackle coupled fluid-solid heat transfer. The extended IB method is validated for Poiseuille flow, which allows for a direct comparison of the numerical results against a closed analytical solution. We subsequently apply the extended IB method to flow in a structured porous medium and focus on bulk properties such as the gradient of the average pressure and the Nusselt number. Reliable qualitative results were obtained with 16-32 grid points per singly-connected fluid region.",

keywords = "Incompressible Navier-Stokes equations, Immersed boundary method, Numerical simulation, Porous media, Volume penalization, Conjugate heat transfer",

author = "{Lopez Penha}, {David J.} and Lilya Ghazaryan and Geurts, {Bernard J.} and Steffen Stolz and Markus Nordlund",

note = "The authors wish to thank Philip Morris International R&D for their financial support, and the Dutch Supercomputing Facility SARA.; 5th European Conference on Computational Fluid Dynamics, ECCOMAS ECFD 2010, ECCOMAS ECFD 2010 ; Conference date: 14-06-2010 Through 17-06-2010",

year = "2010",

month = jun,

language = "English",

isbn = "978-989-96778-1-4",

editor = "J.C.F Pereira and A. Sequeira and J.M.C. Pereira",

booktitle = "Proceedings of the V European Conference on Computational Fluid Dynamics, ECCOMAS CFD 2010",

publisher = "Technical University of Lisbon",

}

Lopez Penha, DJ, Ghazaryan, L, Geurts, BJ, Stolz, S & Nordlund, M 2010, An immersed boundary method for computing heat and fluid flow in porous media. in JCF Pereira, A Sequeira & JMC Pereira (eds), Proceedings of the V European Conference on Computational Fluid Dynamics, ECCOMAS CFD 2010., 01199, Technical University of Lisbon, Lisbon, 5th European Conference on Computational Fluid Dynamics, ECCOMAS ECFD 2010, Lisbon, Portugal, 14/06/10.

An immersed boundary method for computing heat and fluid flow in porous media. / Lopez Penha, David J.; Ghazaryan, Lilya; Geurts, Bernard J. et al.
Proceedings of the V European Conference on Computational Fluid Dynamics, ECCOMAS CFD 2010. ed. / J.C.F Pereira; A. Sequeira; J.M.C. Pereira. Lisbon: Technical University of Lisbon, 2010. 01199.

Research output: Chapter in Book/Report/Conference proceeding › Conference contribution › Academic › peer-review

TY - GEN

T1 - An immersed boundary method for computing heat and fluid flow in porous media

AU - Lopez Penha, David J.

AU - Ghazaryan, Lilya

AU - Geurts, Bernard J.

AU - Stolz, Steffen

AU - Nordlund, Markus

N1 - Conference code: 5

PY - 2010/6

Y1 - 2010/6

N2 - A volume-penalizing immersed boundary (IB) method is presented that facilitates the computation of fluid flow in complex porous media. The computational domain is composed of a uniform Cartesian grid, and solid bodies are approximated on this grid using a series of grid cells (i.e., a ''staircase'' approximation). Solid bodies are distinguished from fluid regions using a binary phase-indicator function: Taking the value of ''1'' in the solid parts of the domain and ''0'' in the fluid parts. The effect of solid bodies on the flow is modeled using a source term in the momentum equations. The source term is active only within solid parts of the domain, and enforces the no-slip boundary condition. Fluid regions are governed by the incompressible Navier-Stokes equations. An extension of the IB method is proposed to tackle coupled fluid-solid heat transfer. The extended IB method is validated for Poiseuille flow, which allows for a direct comparison of the numerical results against a closed analytical solution. We subsequently apply the extended IB method to flow in a structured porous medium and focus on bulk properties such as the gradient of the average pressure and the Nusselt number. Reliable qualitative results were obtained with 16-32 grid points per singly-connected fluid region.

AB - A volume-penalizing immersed boundary (IB) method is presented that facilitates the computation of fluid flow in complex porous media. The computational domain is composed of a uniform Cartesian grid, and solid bodies are approximated on this grid using a series of grid cells (i.e., a ''staircase'' approximation). Solid bodies are distinguished from fluid regions using a binary phase-indicator function: Taking the value of ''1'' in the solid parts of the domain and ''0'' in the fluid parts. The effect of solid bodies on the flow is modeled using a source term in the momentum equations. The source term is active only within solid parts of the domain, and enforces the no-slip boundary condition. Fluid regions are governed by the incompressible Navier-Stokes equations. An extension of the IB method is proposed to tackle coupled fluid-solid heat transfer. The extended IB method is validated for Poiseuille flow, which allows for a direct comparison of the numerical results against a closed analytical solution. We subsequently apply the extended IB method to flow in a structured porous medium and focus on bulk properties such as the gradient of the average pressure and the Nusselt number. Reliable qualitative results were obtained with 16-32 grid points per singly-connected fluid region.

KW - Incompressible Navier-Stokes equations

KW - Immersed boundary method

KW - Numerical simulation

KW - Porous media

KW - Volume penalization

KW - Conjugate heat transfer

UR - http://www.eccomas-cfd2010.org/proceedings.php

M3 - Conference contribution

SN - 978-989-96778-1-4

BT - Proceedings of the V European Conference on Computational Fluid Dynamics, ECCOMAS CFD 2010

A2 - Pereira, J.C.F

A2 - Sequeira, A.

A2 - Pereira, J.M.C.

PB - Technical University of Lisbon

CY - Lisbon

T2 - 5th European Conference on Computational Fluid Dynamics, ECCOMAS ECFD 2010

Y2 - 14 June 2010 through 17 June 2010

ER -

An immersed boundary method for computing heat and fluid flow in porous media

Abstract

Conference

Keywords

Access to Document

Other files and links

Fingerprint

Cite this