Commun. Comput. Phys., 2 (2007), pp. 723-745.

Multiscale Finite Element Modelling of Flow Through Porous Media with Curved and Contracting Boundaries to Evaluate Different Types of Bubble Functions

V. Nassehi 1*, M. Parvazinia 2, A. Khan 1

1 Department of Chemical Engineering, Loughborough University, Loughborough, Leicestershire, LE11 3TU, UK.
2 Department of Mathematical Modelling and Simulation, Iran Polymer and Petrochemical Institute, Tehran, P.O. Box 14185-458, Iran.

Received 15 August 2006; Accepted (in revised version) 9 November 2006
Communicated by Pingwen Zhang
Available online 22 January 2007


The Brinkman equation is used to model the isothermal flow of the Newtonian fluids through highly permeable porous media. Due to the multiscale behaviour of this flow regime the standard Galerkin finite element schemes for the Brinkman equation require excessive mesh refinement at least in the vicinity of domain walls to yield stable and accurate results. To avoid this, a multiscale finite element method is developed using bubble functions. It is shown that by using bubble enriched shape functions the standard Galerkin method can generate stable solutions without excessive near wall mesh refinements. In this paper the performances of different types of bubble functions are evaluated. These functions are used in conjunction with bilinear Lagrangian elements to solve the Brinkman equation via a penalty finite element scheme.

Notice: Undefined variable: ams in /var/www/html/issue/abstract/readabs.php on line 163
PACS: 47.11.-j, 47.11.Fg, 47.11.St
Key words: Finite element, multiscale method, porous media, Newtonian fluid flow, bubble function, static condensation.

*Corresponding author.
Email: (V. Nassehi), (M. Parvazinia), (A. Khan)

The Global Science Journal