Darcy's law–based numerical simulation for modeling 3D liquid absorption into porous wicks

Liquid imbibition into polymer wicks, where a clear liquid front can be seen rising during the wicking process, is modeled using the concepts of flow in porous media. The flow of liquid behind the moving liquid front is modeled using the physics of single-phase flow in a porous medium where the Darcy’s law is combined with the continuity equation and a capillary suction pressure is imposed at the liquid front. A novel numerical simulation PORE-FLOW V C based on the finite element/control volume method is proposed to model such imbibitional flows in wicks of complex shapes. A validation of the simulation is obtained by achieving an excellent comparison of its predictions with an experimental result, an analytical solution, and the Washburn equation for the case of wicking against gravity in a cylindrical wick. The simulation is also used to predict a case of two-dimensional (2D) wicking in the altered cylindrical wicks with two different cross-sectional areas. Once again an excellent match is obtained with the experimental results, while analytical solutions for the single and double cross-section cases along with the Washburn equation fail to predict the 2D wicking. Later, some other types of altered wicks with sharp changes in their cross-sectional areas were analyzed numerically for their wicking behavior. It was observed that the height of liquid front in a vertical wick as a function of time, which is proportional to the history of liquid imbibed, is strongly dependent on the extent of reduction in the wick cross-sectional area as well as its location vis-a `-vis the wick entrance. V C 2010 American Institute of Chemical Engineers AIChE J, 57: 1132–1143, 2011

[1]  Reza Masoodi,et al.  Role of Hydraulic and Capillary Radii in Improving the Effectiveness of Capillary Model in Wicking , 2008 .

[2]  H. Tan,et al.  Finite element implementation of stress-jump and stress-continuity conditions at porous-medium, clear-fluid interface , 2009 .

[3]  J. Szekely,et al.  The rate of capillary penetration and the applicability of the washburn equation , 1971 .

[4]  S. Whitaker.  The method of volume averaging , 1998 .

[5]  Adrian E. Scheidegger,et al.  The physics of flow through porous media , 1957 .

[6]  E. Kissa,et al.  Wetting and Wicking , 1996 .

[7]  F. Dullien.  Porous Media: Fluid Transport and Pore Structure , 1979 .

[8]  Suresh G. Advani,et al.  Wicking across a Fiber-Bank , 1996 .

[9]  J. Parlange,et al.  Capillary absorption in porous sheets and surfaces subject to evaporation , 2007 .

[10]  Richard Lucas,et al.  Ueber das Zeitgesetz des kapillaren Aufstiegs von Flüssigkeiten , 1918 .

[11]  J. Bear.  Dynamics of Fluids in Porous Media , 1975 .

[12]  E. W. Washburn.  The Dynamics of Capillary Flow , 1921 .

[13]  Patrick R. Amestoy,et al.  Multifrontal parallel distributed symmetric and unsymmetric solvers , 2000 .

[14]  N. Pan,et al.  Thermal and moisture transport in fibrous materials , 2006 .

[15]  H. Tan,et al.  Effect of fiber‐mat anisotropy on 1D mold filling in LCM: A numerical investigation , 2008 .

[16]  R. Horton,et al.  An Analytical Solution for One-Dimensional Water Infiltration and Redistribution in Unsaturated Soil , 2009 .

[17]  J. Parlange,et al.  Wetting front instability in unsaturated porous media: A three-dimensional study in initially dry sand , 1990 .

[18]  K. Pillai,et al.  Darcy's law‐based models for liquid absorption in polymer wicks , 2007 .

[19]  R. Skalak,et al.  THE HISTORY OF POISEUILLE'S LAW , 1993 .