Finite Element Modelling of Flow Through a Porous Medium Between Two Parallel Plates Using The Brinkman Equation

Despite the widespread use of the Darcy equation to model porous flow, it is well known that this equation is inconsistent with commonly prescribed no slip conditions at flow domain walls or interfaces between different sections. Therefore, in cases where the wall effects on the flow regime are expected to be significant, the Darcy equation which is only consistent with perfect slip at solid boundaries, cannot predict velocity and pressure profiles properly and alternative models such as the Brinkman equation need to be considered. This paper is devoted to the study of the flow of a Newtonian fluid in a porous medium between two impermeable parallel walls at different Darcy parameters (Da). The flow regime is considered to be isothermal and steady. Three different flow regimes can be considered using the Brinkman equation: free flow (Da > 1), porous flow (high permeability, 1 > Da > 10−6) and porous flow (low permeability Da < 10−6). In the present work the described bench mark problem is used to study the effects of solid walls for a range of low to high Darcy parameters. Both no-slip and slip conditions are considered and the results of these two cases are compared. The range of the applicability of the Brinkman equation and simulated results for different cases are shown.

[1]  Fathi M. Allan,et al.  Fluid mechanics of the interface region between two porous layers , 2002, Appl. Math. Comput..

[2]  J. N. Reddy,et al.  Applied Functional Analysis and Variational Methods in Engineering , 1986 .

[3]  Vahid Nassehi,et al.  Practical Aspects of Finite Element Modelling of Polymer Processing , 2002 .

[4]  M. Kaviany Laminar flow through a porous channel bounded by isothermal parallel plates , 1985 .

[5]  Massoud Kaviany,et al.  Non-Darcian effects on natural convection in porous media confined between horizontal cylinders , 1986 .

[6]  V. Sugunamma,et al.  Finite Element Analysis of Convection Flow Through a Porous Medium in a Horizontal Channel , 1999 .

[7]  Curtis F. Gerald Applied numerical analysis , 1970 .

[8]  V. Nassehi,et al.  Modeling the transient flow of rubber compounds in the dispersive section of an internal mixer with slip‐stick boundary conditions , 1997 .

[9]  J. Z. Zhu,et al.  The finite element method , 1977 .

[10]  G. Lauriat,et al.  Analytical solution of non-Darcian forced convection in an annular duct partially filled with a porous medium , 1995 .

[11]  Bruce M. Irons,et al.  A frontal solution program for finite element analysis , 1970 .

[12]  Kankanhalli N. Seetharamu,et al.  Finite element analysis of heat transfer by natural convection in porous media in vertical enclosures: Investigations in Darcy and non‐Darcy regimes , 1997 .

[13]  C. L. Tien,et al.  Boundary and inertia effects on flow and heat transfer in porous media , 1981 .

[14]  C. Hsu,et al.  The Brinkman model for natural convection about a semi-infinite vertical flat plate in a porous medium , 1985 .

[15]  Kambiz Vafai,et al.  Analysis of flow and heat transfer at the interface region of a porous medium , 1987 .

[16]  A. Mohamad,et al.  Non-Darcy effects in buoyancy driven flows in an enclosure filled with vertically layered porous media , 2002 .

[17]  A. Bejan,et al.  Convection in Porous Media , 1992 .

[18]  M. Kaviany,et al.  Non-Darcian effects on vertical-plate natural convection in porous media with high porosities , 1985 .