{"id":206,"date":"2020-12-09T12:54:59","date_gmt":"2020-12-09T17:54:59","guid":{"rendered":"https:\/\/mae.ncsu.edu\/pmmf\/?page_id=206"},"modified":"2024-12-20T14:59:40","modified_gmt":"2024-12-20T19:59:40","slug":"fluid-absorption-in-porous-media","status":"publish","type":"page","link":"https:\/\/mae.ncsu.edu\/pmmf\/fluid-absorption-in-porous-media\/","title":{"rendered":"Two-Phase Flow in Porous Media"},"content":{"rendered":"\r\n
\r\n

 <\/p>\r\n\r\n\r\n\r\n\r\n\r\n\r\n\r\n\r\n\r\n\r\n\r\n\r\n\r\n\r\n\r\n\r\n\r\n\r\n\r\n\r\n\r\n\r\n\r\n\r\n\r\n\r\n\r\n\r\n\r\n\r\n
Capillarity is the driving force that plays an important role in the displacement of one fluid by another at the microscale on a solid wall. Understanding and quantifying capillary forces is crucial in designing products and applications that deal with immiscible fluids. In context of fluid transport in fibrous materials, the rate of fluid absorption is often predicted using the Lucas\u2013Washburn (LW) equation. Typical LW models are developed on the basis of treating a porous medium as bundles of 1-D parallel capillary tubes with a given single-value contact angle. A more accurate (but more rigorous), model for fluid transport is the Richards equation<\/strong>, which requires mathematical relationships for the capillary pressure pc<\/sub>=pc<\/sub>(S) and relative permeability K(S) of the porous media as a function of its instantaneous fluid saturation.<\/span><\/td>\r\n<\/tr>\r\n
\"\" Richards equation can be solved numerically to predict fluid saturation as a function of time and space S=S(x,y,z,t) in a fluid absorbing (or fluid releasing) porous material. Instantaneous saturation is needed to predict the rate of fluid absorption into (or release from) the material [Jaganathan et al., 2009].<\/span><\/em><\/span><\/td>\r\n<\/tr>\r\n
Capillary pressure<\/strong> is the pressure required to force a nonwetting fluid into a porous media (e.g., forcing air into a hydrophilic media initially filled with water) as can be seen in the figure below.\u00a0<\/span><\/td>\r\n<\/tr>\r\n
\"\"2-D numerical simulation of water drainage from a fibrous media comprised of fibers with different degrees of hydrophilicity. The red and black fibers have a Young-Laplace contact angle (YLCA) of 65 and 35 degrees, respectively. It can be seen that fluid saturation percentage decreases (from left to right) by increasing the intrusion (capillary) pressure [Bucher and Tafreshi 2014].<\/span><\/span><\/em><\/td>\r\n<\/tr>\r\n
Such simulations can be used to produce a mathematical relationship between capillary pressure and fluid saturation in single- or multi-component (comprised of fibers with different fiber diameter of contact angles) as can be seen in the figure below.<\/span><\/td>\r\n<\/tr>\r\n
\"\"Effects of reducing the mass fraction of hydrophilic fibers (from A to E) on capillary pressure of bicomponent media comprised of hydrophilic (black) and relatively hydrophobic (red) fibers. All media have the same porosity and fiber diameters and are under the same intrusion pressure of 7 kPa [Bucher and Tafreshi 2014].<\/span><\/span><\/em><\/td>\r\n<\/tr>\r\n
\r\n

A fast and flexible method to estimate the capillary pressure of a 3-D fibrous media is the Pore Morphology Method (PMM). PMM employs a voxel-based approach to present a quasi-static description of how a non-wetting phase penetrates into a porous medium filled with a wetting phase. PMM simulations are faster than the convention methods (e.g., the Volume-of-Fluid or Lattice Boltzmann method)\u00a0 by orders of magnitude as the PMM simulations are merely geometry-based calculations (i.e., do not require solving the Navier\u2013Stokes equations). The PMM approach can make 3-D two-phase flow simulations computationally affordable for industrial applications.\u00a0<\/span>PMM relies on two primary morphological operations, referred to as dilation and erosion, to simulate the spatial distribution of the phases in response to an increase in the non-wetting phase pressure. The pictures below show example results obtained from\u00a0 our in-house PMM computer code [Gautam et al. 2025].<\/span><\/p>\r\n<\/td>\r\n<\/tr>\r\n

\"\"
Our MATLAB implementation of the PMM algorithm is demonstrated a 2-D periodic domain comprised of nonspherical solid objects (a).\u00a0 The dilation process is shown in (b) and (c). The no-wetting phase voxels isolated from their reservoir are turned into wetted voxels in (d). The erosion process is shown in (e). The final result is shown in (f) for the given intrusion pressure [Gautam et al., 2025].<\/span><\/em><\/td>\r\n<\/tr>\r\n
The PMM algorithm was then used in larger geometries to conduct realistic simulations of fluid saturation\/desaturation in 2-D fibrous media as shown the figure below.\u00a0 \u00a0<\/td>\r\n<\/tr>\r\n
\r\n
\"\"<\/figure>\r\n

Simulation of pressurized air intruding into a fully-wetted fibrous media is shown. The fibers have two different diameters of 10 \u00b5m and 20 \u00b5m and two different Young-Laplace contact angles of 20 and 60 degrees. The bottom figure is from the semi-numerical simulations of Bucher and Tafreshi (2014) and the top figure is the PMM simulation results\u00a0 generated using MATLAB [Gautam et al., 2025].\u00a0<\/span><\/em><\/p>\r\n<\/td>\r\n<\/tr>\r\n

Our ongoing effort is now focused on extending the above algorithm to the case of 3-D geometries. The figures below show our preliminary PMM simulations in 3-D.\u00a0<\/td>\r\n<\/tr>\r\n
\r\n

\"\"<\/p>\r\n

The figure above shows preliminary results for the MATLAB implementation of the PMM algorithm in a 3-D packed bed [Gautam et al. 2024].\u00a0<\/span><\/em><\/p>\r\n<\/td>\r\n<\/tr>\r\n

Preliminary PMM simulation results are shown below for desaturation of a 3-D fibrous medium comprised of fibers with a diameter of 20 \u00b5m and an YLCA of 10 degrees. The wetting phase here is water (shown in blue).<\/td>\r\n<\/tr>\r\n
\r\n

\"\"<\/p>\r\n

The figure above shows preliminary results for the MATLAB implementation of the PMM algorithm in a 3-D fibrous domain [Bhatta et al. 2024]<\/span><\/em><\/p>\r\n<\/td>\r\n<\/tr>\r\n

The PMM simulation results show good qualitative agreement with those obtained via VOF method as shown below.\u00a0\u00a0<\/td>\r\n<\/tr>\r\n
\r\n

\"\"<\/p>\r\n

The figure compares the VOF and PMM simulations of bubble point formation in a 3-D fibrous domain [Bhatta et al. 2024]<\/span><\/em><\/p>\r\n<\/td>\r\n<\/tr>\r\n

The PMM simulations can also handle fluid penetration in media with heterogenous YLCAs and fiber diameter. The picture below is an example of such calculations.\u00a0<\/span><\/td>\r\n<\/tr>\r\n
\"\" The YLCA for the red and green fibers is 10 and 80 degrees, respectively.\u00a0 Fiber diameter is 20\u00b5m [Bhatta et al. 2024].<\/span><\/em><\/span><\/td>\r\n<\/tr>\r\n
Despite the advances in computational methods, capillary pressure-saturation relationship for media comprised of fluid-storing swelling-contracting porous fibers (e.g., Rayon or cotton fibers) can only be found experimentally. Figure below shows an example of such an experiment [Jaganathan et al., 2009].<\/span><\/td>\r\n<\/tr>\r\n
\"\"Fluid height rise experiment is conducted to obtain a capillary pressure-saturation relationship for nonwoven fabrics made of Rayon (swelling) fibers [Ashari et al., 2010].<\/span><\/em><\/td>\r\n<\/tr>\r\n
Relative permeability<\/strong> is the inverse of the resistance that a partially-wetted media shows against fluid flow. The relative permeability tensor can be obtained by solving the Stokes equations numerically at each saturation levels and post processing the resulting data using Darcy’s law. It is important to note that relative permeability (unlike capillary pressure) is a tensorial property. This means that fluid absorption in the x,y and z directions can be quite different in a fibrous media dependent on the x,y, and z orientation of the fibers. This also means that fibrous media with different overall fiber orientations (but identical solid volume fractions and fiber diameters) can absorb (or release) fluids at different rates [Jaganathan et al. 2008].\u00a0<\/span><\/td>\r\n<\/tr>\r\n
\"\"\u00a0 \u00a0\u00a0\"\"
Relative permeability tensor can be obtained by solving the Stokes equations for flow of a wetting fluid in a fibrous media under different saturation levels. Figure to the right is an example of relative permeability calculation performed for the 2-D macroscale simulations shown in the next figure below [Ashari et al. 2010].
<\/em><\/span><\/td>\r\n<\/tr>\r\n
With relative permeability and capillary pressure obtained from the above 3-D microscale (fiber-level) simulations in terms of wetting saturation, we solved the Richards equation on macroscale (scales comparable with dimesions of the media) using a finite element PDE solver [Jaganathan et al., 2009]. Figures below show an example of such dual-scale calculations conducted for fluid absorption in fibrous sheets with identical parameters but different in-plane fiber orientations.
<\/span><\/td>\r\n<\/tr>\r\n
\"\"Contour plots of saturation at t = 0.41 s planar media with different in-plane fiber orientations. Different colors from red to blue represent different saturation values from one to zero, respectively. Coordinates are normalized by the actual sample dimensions. The volume of the absorbed liquid in these media is plotted versus time. The figure to the left shows that the rate of fluid spread increases with increasing the in-plane alignment of the fibers.\u00a0 [Ashari et al. 2010].\u00a0<\/i><\/td>\r\n<\/tr>\r\n
We also used our dual-scale modeling approach to simulate the rate of fluid release<\/strong> from moving partially saturated nonwoven sheets in contact with a <\/span>solid surface. This is a challenging task as the release rate depends on many parameters, some of which are difficult to quantify. In this concern, we developed a diffusion-controlled boundary treatment to simulate fluid release from partially saturated porous materials onto surfaces with different hydrophilicity. Our numerical simulations were accompanied with experimental data obtained for lotion-saturated Rayon and PEY nonwoven wipes using a custom-made test rig.\u00a0<\/span><\/td>\r\n<\/tr>\r\n
\"\"The figure on the left is a schematic illustration of our experimental rig built to test nonwoven wet wipes. The figure on the right compares wipes\u2019 average saturation obtained from our simulations with those obtained experimentally for Rayon wipes at at different speeds [Ashari et al., 2011].<\/span><\/em><\/td>\r\n<\/tr>\r\n
Representative Publications:<\/strong>\r\n
    \r\n
  • S. Jaganathan, H.V. Tafreshi, and B. Pourdeyhimi, A Realistic Modeling of Fluid Infiltration in Thin Fibrous Sheets, Journal of Applied Physics<\/strong><\/em> 105, 113522 (2009) <\/span><\/li>\r\n
  • T.M. Bucher, H.V. Tafreshi, Modeling Air\u2013Water Interface in Disordered Fibrous Media with Heterogeneous Wettabilities, Colloids and Surfaces A<\/em><\/strong> 461, 323 (2014)<\/span><\/li>\r\n
  • S. Jaganathan, H.V. Tafreshi, and B. Pourdeyhimi, A Case Study of Realistic Two-Scale Modeling of Water Permeability in Fibrous Media, Separation Science & Technology<\/em> <\/strong>43, 1901 (2008)<\/span><\/li>\r\n
  • A. Ashari, T.M. Bucher, H.V. Tafreshi, M. A. Tahir, and M.S.A. Rahman, Modeling Fluid Spread in Thin Fibrous Sheets: Effects of Fiber Orientation, International Journal of Heat & Mass Transfer<\/strong><\/em> 53, 1750 (2010)<\/span><\/li>\r\n
  • A. Ashari, M. Bucher, and H.V. Tafreshi, A Semi-Analytical Model for Simulating Fluid Transport in Multi-Layered Fibrous Sheets Made up of Solid and Porous Fibers, Computational Materials Science<\/em> <\/strong>50, 378-390 (2010) <\/span><\/li>\r\n
  • A. Ashari and H.V. Tafreshi, General Capillary Pressure and Relative Permeability Expressions for Through-Plane Fluid Transport in Thin Fibrous Sheets, Colloids and Surfaces A<\/em><\/strong> 346, 114 (2009)<\/span><\/li>\r\n
  • A. Ashari, T.M. Bucher, and H.V. Tafreshi, Modeling Motion-Induced Fluid Release from Partially-Saturated Fibrous Media onto Surfaces with Different Hydrophilicity, International Journal of Heat and Fluid Flow<\/em><\/strong> 32, 1076 (2011)<\/span><\/li>\r\n
  • N. Bhatta, S. Gautam, H.V. Tafreshi, and B. Pourdeyhimi,\u00a0 Accuracy of Pore Morphology Method in Modeling Fluid Saturation in 3-D Fibrous Media, Ind. Eng. Chem. Res.<\/strong><\/em>, 63, 18147 (2024)
    <\/span><\/li>\r\n
  • S. Gautam, N. Bhatta, H.V. Tafreshi, and B. Pourdeyhimi, MATLAB Implementation of Pore Morphology Method for Modeling Liquid Residue in Porous Media with Heterogenous Wettabilities, Powder Technology<\/strong><\/em>, 452, 120509 <\/span><\/span>(2025)<\/span>\u00a0<\/li>\r\n<\/ul>\r\n

    Click here for a more complete list of publications<\/a><\/span><\/p>\r\n<\/td>\r\n<\/tr>\r\n

Acknowledgement:<\/strong>\r\n
    \r\n
  • Ahlstrom Company
    <\/span><\/li>\r\n
  • The Nonwovens Institute<\/span><\/li>\r\n
  • Johnson & Johnson Company<\/span><\/li>\r\n
  • Kenan Institute<\/span><\/li>\r\n<\/ul>\r\n<\/td>\r\n<\/tr>\r\n<\/tbody>\r\n<\/table>\r\n<\/figure>\r\n

    <\/p>\r\n

     <\/p>\r\n

     <\/p>","protected":false},"excerpt":{"rendered":"

      Capillarity is the driving force that plays an important role in the displacement of one fluid by another at the microscale on a solid…<\/p>\n","protected":false},"author":394,"featured_media":0,"parent":0,"menu_order":0,"comment_status":"closed","ping_status":"closed","template":"","meta":{"_acf_changed":false,"footnotes":""},"class_list":["post-206","page","type-page","status-publish","hentry"],"acf":[],"_links":{"self":[{"href":"https:\/\/mae.ncsu.edu\/pmmf\/wp-json\/wp\/v2\/pages\/206","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/mae.ncsu.edu\/pmmf\/wp-json\/wp\/v2\/pages"}],"about":[{"href":"https:\/\/mae.ncsu.edu\/pmmf\/wp-json\/wp\/v2\/types\/page"}],"author":[{"embeddable":true,"href":"https:\/\/mae.ncsu.edu\/pmmf\/wp-json\/wp\/v2\/users\/394"}],"replies":[{"embeddable":true,"href":"https:\/\/mae.ncsu.edu\/pmmf\/wp-json\/wp\/v2\/comments?post=206"}],"version-history":[{"count":10,"href":"https:\/\/mae.ncsu.edu\/pmmf\/wp-json\/wp\/v2\/pages\/206\/revisions"}],"predecessor-version":[{"id":1269,"href":"https:\/\/mae.ncsu.edu\/pmmf\/wp-json\/wp\/v2\/pages\/206\/revisions\/1269"}],"wp:attachment":[{"href":"https:\/\/mae.ncsu.edu\/pmmf\/wp-json\/wp\/v2\/media?parent=206"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}