Application of an immersed boundary method to flow in cerebral aneurysms and porous media


Mikhal, Julia and Lopez Penha, David J. and Stolz, Steffen and Geurts, Bernard J. (2010) Application of an immersed boundary method to flow in cerebral aneurysms and porous media. In: 3rd Joint US-European ASME Fluids Engineering Summer Meeting and the 8th International Conference on Nanochannels, Micorchannels and Minichannels, FEDSM-ICNMM 2010, 1-5 Aug 2010, Montreal, Quebec, Canada (pp. p. 31042).

[img] PDF
Restricted to UT campus only
: Request a copy
Abstract:We present the development and application of an immersed boundary (IB) method for the simulation of incompressible flow inside and around complex geometrical shapes and cavities. The IB method is based on a volume-penalization method that is applied throughout the domain, rendering the velocity in stationary solid parts negligibly small, while the flow in the open parts of the domain is governed by the Navier-Stokes equations. The flow solver is based on a skew-symmetric finite-volume discretization in combination with explicit time-stepping for the convective and viscous fluxes, and implicit time-stepping for the IB forcing term. The complex geometry is characterized in terms of a so-called ‘masking function’ which equals unity in the solid parts and zero in the open parts of the domain. The focus is on the accuracy with which gradients of the solution close to solid walls can be approximated using the IB methodology. We investigate this for flow through a model of an aneurysm as may develop in the circle of Willis in a human brain, and to flow in a structured porous medium composed of a regular spatial arrangement of square rods. The shear stress acting on the vessel walls in case of flow through an aneurysm and the permeability of the porous material were analyzed. The computational method converges as a first order method for Poiseuille flow, with a considerable influence derived from the precise definition of the masking function near solid-fluid interfaces. We identify the best masking function strategy and show that for plane Poiseuille flow even second order convergence may be obtained. Qualitatively reliable results are obtained already at modest resolutions of 8-16 grid cells across a characteristic opening in the flow domain, e.g., the vessel diameter or the size of the gap between individual square rods.
Item Type:Conference or Workshop Item
Copyright:© 2010 ASME
Electrical Engineering, Mathematics and Computer Science (EEMCS)
Research Group:
Link to this item:
Conference URL:
Export this item as:BibTeX
HTML Citation
Reference Manager


Repository Staff Only: item control page