Numerical Modeling of Capillary-Driven Flow in Open Microchannels: an Implication of Optimized Wicking Fabric Design

Total Page:16

File Type:pdf, Size:1020Kb

Numerical Modeling of Capillary-Driven Flow in Open Microchannels: an Implication of Optimized Wicking Fabric Design Scholars' Mine Masters Theses Student Theses and Dissertations Summer 2018 Numerical modeling of capillary-driven flow in open microchannels: An implication of optimized wicking fabric design Mehrad Gholizadeh Ansari Follow this and additional works at: https://scholarsmine.mst.edu/masters_theses Part of the Environmental Engineering Commons, Mathematics Commons, and the Mechanical Engineering Commons Department: Recommended Citation Ansari, Mehrad Gholizadeh, "Numerical modeling of capillary-driven flow in open microchannels: An implication of optimized wicking fabric design" (2018). Masters Theses. 7792. https://scholarsmine.mst.edu/masters_theses/7792 This thesis is brought to you by Scholars' Mine, a service of the Missouri S&T Library and Learning Resources. This work is protected by U. S. Copyright Law. Unauthorized use including reproduction for redistribution requires the permission of the copyright holder. For more information, please contact [email protected]. NUMERICAL MODELING OF CAPILLARY-DRIVEN FLOW IN OPEN MICROCHANNELS: AN IMPLICATION OF OPTIMIZED WICKING FABRIC DESIGN by MEHRAD GHOLIZADEH ANSARI A THESIS Presented to the Faculty of the Graduate School of the MISSOURI UNIVERSITY OF SCIENCE AND TECHNOLOGY In Partial Fulfillment of the Requirements for the Degree MASTER OF SCIENCE IN ENVIRONMENTAL ENGINEERING 2018 Approved by Dr. Wen Deng, Advisor Dr. Joseph Smith Dr. Jianmin Wang Dr. Xiong Zhang 2018 Mehrad Gholizadeh Ansari All Rights Reserved iii PUBLICATION THESIS OPTION This thesis has been formatted using the publication option: Paper I, pages 16-52, are intended for submission to the Journal of Computational Physics. iv ABSTRACT The use of microfluidics to transfer fluids without applying any exterior energy source is a promising technology in different fields of science and engineering due to their compactness, simplicity and cost-effective design. In geotechnical engineering, to increase the soil’s strength, hydrophilic wicking fibers as type of microfluidics have been employed to transport and drain water out of soil spontaneously by taking advantage of natural capillary force without using any pumps or other auxiliary devices. The objective of this study is to understand the scientific mechanisms of the capability for wicking fiber to drain both gravity and capillary water out of soil fills. In this work, wicking fibers were numerically modeled as open microchannels. The effect of the geometry and wettability on the spontaneous flow were analyzed. A 3D computational fluid dynamics-based model was developed to predict the transient movement of the meniscus in microchannels of different geometries. The Volume-of-Fluid (VOF) method was employed along with an in-house developed algorithm to refine the grid only at the air-water interface to reduce the computational effort. Analytical and numerical criteria were derived to determine the critical contact angle for spontaneous flow in different geometry open microchannels. It was found that the aspect ratio and channel geometry have a significant effect on the water filling velocity and mass flowrate. The study was extended to water drainage from unsaturated soils and the minimum negative inlet pressure which represents the maximum soil suction was determined numerically for spontaneous flow for different microchannel geometry. Results of this study can provide an implication in the design of wicking fibers and other microfluidic systems. v ACKNOWLEDGMENTS I would like to express my sincere gratitude to my advisor, Dr. Wen Deng, for his guidance, patience, and support. My special word of thanks goes to my committee members, Dr. Joseph Smith, Dr. Jianmin Wang and Dr. Xiong Zhang for their support of this research. I express my heart-felt gratitude to Dr. Joseph Smith for granting me access to the STAR-CCM+ license servers. His constant guidance, cooperation, motivation and support have always kept me going ahead. I owe my deepest gratitude towards my family, Farimah, Hamid and Mahta for their eternal support and understanding of my goals and aspirations. Their infallible love and support has always been my strength. A special thanks to my friends, Kiyavash, Kiarash, Behrouz and Elieh and all the other Amoos for their love, motivation and affection. vi TABLE OF CONTENTS Page PUBLICATION THESIS OPTION ............................................................................... iii ABSTRACT ................................................................................................................... iv ACKNOWLEDGMENTS ............................................................................................... v LIST OF FIGURES ...................................................................................................... viii LIST OF TABLES ........................................................................................................... x NOMENCLATURE ....................................................................................................... xi SECTION 1. INTRODUCTION .................................................................................................... 1 1.1 APPLICATIONS ............................................................................................. 1 1.2 SURFACE TENSION ..................................................................................... 5 1.3 YOUNG-LAPLACE LAW ............................................................................. 6 1.4 DIMENSIONLESS NUMBERS ..................................................................... 7 1.5 OBJECTIVES ................................................................................................. 8 2. METHODOLOGY ................................................................................................. 10 2.1 ADAPTIVE GRID REFINEMENT .............................................................. 11 2.2 ADAPTIVE TIME STEP .............................................................................. 14 PAPER I I. NUMERICAL MODELING OF CAPILLARY-DRIVEN FLOW IN OPEN MICROCHANNELS ........................................................ 16 ABSTRACT ........................................................................................................... 16 1. INTRODUCTION .......................................................................................... 17 vii 1.1 MOTIVATION ...................................................................................... 17 1.2 BACKGROUND ................................................................................... 17 1.3 SPONTANEOUS CAPILLARY-DRIVEN FLOW .............................. 22 2. NUMERICAL ANALYSIS OF CAPILLARY-DRIVEN FLOW .................. 26 2.1 GOVERNING EQUATIONS ................................................................ 26 2.2 DIFFERENCING SCHEMES ............................................................... 29 2.3 MODEL VALIDATION ....................................................................... 32 3. RESULTS AND DISCUSSION ..................................................................... 34 3.1 SCF CRITERIA ..................................................................................... 34 3.2 MOBILITY PARAMETER ................................................................... 36 3.3 NEGATIVE INLET PRESSURE .......................................................... 43 4. CONCLUSIONS ........................................................................................... 48 REFERENCES ....................................................................................................... 49 SECTION 3. CONCLUSIONS .................................................................................................... 53 4. FUTURE WORK ................................................................................................... 54 APPENDIX .................................................................................................................. 55 REFERENCES .............................................................................................................. 71 VITA ............................................................................................................................. 73 viii LIST OF FIGURES SECTION Page Figure 1.1. Blood collection device for Point Of Care Testing (POCT) [3]. ...................... 1 Figure 1.2. Schematic design of a self-pumping fuel cell[4]. ............................................. 2 Figure 1.3. Mechanism of energy transfer in a heat pipe[7]. .............................................. 3 Figure 1.4. Hydrophilic wicking fibers removing water from soil[10]. ............................. 4 Figure 1.5. Wicking fabric application to drain water from a road embankment. .............. 5 Figure 1.6. Shape equilibrium of a liquid droplet (L) on a solid surface (S), surrounded by the gas (G)................................................................................. 5 Figure 1.7. Schematic of the minimum and maximum principal osculating circle radius of the curvature at the interface [13]. ........................................... 7 Figure 2.1. Computational grid in an open circular microchannel. .................................. 12 Figure 2.2. Adaptive grid refinement algorithm. .............................................................. 13 Figure 2.3. Adaptive time step algorithm based on the maximum Co number in the computational domain. .......................................................................... 15 PAPER I Figure 1. Contact angle
Recommended publications
  • Wettability As Related to Capillary Action in Porous Media
    Wettability as Related to Capillary Action in Porous l\1edia SOCONY MOBIL OIL CO. JAMES C. MELROSE DALLAS, TEX. ABSTRACT interpreted with the aid of a model employing the concept of a cylindrical capillary tube. This The contact angle is one of the boundary approach has enjoyed a certain degree of success Downloaded from http://onepetro.org/spejournal/article-pdf/5/03/259/2153830/spe-1085-pa.pdf by guest on 26 September 2021 conditions for the differential equation specifying in correlating experimental results. 13 The general­ the configuration of fluid-fluid interfaces. Hence, ization of this model, however, to situations which applying knowledge concerning the wettability of a involve varying wettability, has not been established solid surface to problems of fluid distribution in and, in fact, is likely to be unsuccessful. porous solids, it is important to consider the In this paper another approach to this problem complexity of the geometrical shapes of the indi­ will be discussed. A considerable literature relating vidual, interconnected pores. to this approach exists in the field of soil science, As an approach to this problem, the ideal soil where it is referred to as the ideal soil model. model introduced by soil physicists is discussed Certain features of this model have also been in detail. This model predicts that the pore structure discussed by Purcell14 in relation to variable of typical porous solids will lead to hysteresis wettability. The application of this model, however, effects in capillary pressure, even if a zero value to studying the role of wettability in capillary of the contact angle is maintained.
    [Show full text]
  • Capillary Action in Celery Carpet Again and See If You Can the Leaf Is in the Sun
    Walking Water: Set aside the supplies you will need which are the 5 plastic cups (or 5 glasses or Go deeper and get scientific and test the differences between jars you may have), the three vials of food coloring (red, yellow & blue) and 4 the different types of paper towels you have in your kit (3), sheets of the same type of paper towel. You will also need time for this demo to the different food colorings and different levels of water in develop. your glass. Follow the directions here to wrap up your Fold each of the 4 paper towels lengthwise about 4 times so that you have nice understanding of capillary action: long strips. Line your five cups, glasses or jars up in a line or an arc or a circle. Fill the 1st, 3rd, and 5th cup about 1/3 to ½ full with tap water. Leave the in-between https://www.whsd.k12.pa.us/userfiles/1587/Classes/74249/c cups empty. Add just a couple or a few drops of food color to the cups with water appilary%20action.pdf – no mixing colors, yet. Put your strips of paper towel into the cups so that one end goes from a wet cup Leaf Transpiration to a dry cup and that all cups are linked by paper towel bridges. Watch the water Bending Water Does water really leave the plant start to wick up the paper. This is capillary action which involves adhesion and Find the balloons in the bags? through the leaves? cohesion properties of water.
    [Show full text]
  • Surface Tension Bernoulli Principle Fluid Flow Pressure
    Lecture 9. Fluid flow Pressure Bernoulli Principle Surface Tension Fluid flow Speed of a fluid in a pipe is not the same as the flow rate Depends on the radius of the pipe. example: Low speed Same low speed Large flow rate Small flow rate Relating: Fluid flow rate to Average speed L v is the average speed = L/t A v Volume V =AL A is the area Flow rate Q is the volume flowing per unit time (V/t) Q = (V/t) Q = AL/t = A v Q = A v Flow rate Q is the area times the average speed Fluid flow -- Pressure Pressure in a moving fluid with low viscosity and laminar flow Bernoulli Principle Relates the speed of the fluid to the pressure Speed of a fluid is high—pressure is low Speed of a fluid is low—pressure is high Daniel Bernoulli (Swiss Scientist 1700-1782) Bernoulli Equation 1 Prr v2 gh constant 2 Fluid flow -- Pressure Bernoulli Equation 11 Pr v22 r gh P r v r gh 122 1 1 2 2 2 Fluid flow -- Pressure Bernoulli Principle P1 P2 Fluid v2 v1 if h12 h 11 Prr v22 P v 122 1 2 2 1122 P11r v22 r gh P r v r gh 1rrv 1 v 1 P 2 P 2 2 22222 1 1 2 if v22 is higher then P is lower Fluid flow Venturi Effect –constricted tube enhances the Bernoulli effect P1 P3 P2 A A v1 1 v2 3 v3 A2 If fluid is incompressible flow rate Q is the same everywhere along tube Q = A v A v 1 therefore A1v 1 = A2 2 v2 = v1 A2 Continuity of flow Since A2 < A1 v2 > v1 Thus from Bernoulli’s principle P1 > P2 Fluid flow Bernoulli’s principle: Explanation P1 P2 Fluid Speed increases in smaller tube Therefore kinetic energy increases.
    [Show full text]
  • On the Numerical Study of Capillary-Driven Flow in a 3-D Microchannel Model
    Copyright © 2015 Tech Science Press CMES, vol.104, no.4, pp.375-403, 2015 On the Numerical Study of Capillary-driven Flow in a 3-D Microchannel Model C.T. Lee1 and C.C. Lee2 Abstract: In this article, we demonstrate a numerical 3-D chip, and studied the capillary dynamics inside the microchannel. We applied the level set method on the Navier-Stokes equation which incorporates the surface tension and two-phase flow characteristics. We analyzed the capillary dynamics near the junction of two microchannels. Such a highlighting point is important that it not only can provide the information of interface behavior when fluids are made into a head-on collision, but also emphasize the idea for the design of the chip. In addition, we study the pressure distribution of the fluids at the junction. It is shown that the model can produce nearly 2000 Pa pressure difference to help push the water through the microchannel against the air. The nonlinear interaction between capillary flows is recorded. Such a nonlinear phenomenon, to our knowledge, occurs due to the surface tension takes action with the wetted wall boundaries in the channel and the nonlinear governing equations for capillary flow. Keywords: Microchannel, Navier-Stokes equation, Capillary flow, Finite ele- ment analysis, two-phase flow. 1 Introduction Microfluidics system saw the importance of development of the bio-chip, which essentially is a miniaturized laboratory that can perform hundreds or thousands of simultaneous biochemical reactions. However, it is almost impossible to find any instrument to measure or detect the model, in particular, to comprehend the internal flow dynamics for the microchannel model, due to the model has been downsizing to a micro/nano- meter scale.
    [Show full text]
  • Capillary Action
    UNIT ONE • LESSON EIGHT Pert. 2 Group Discussion • Did you get similar results for the celery and Whet. IP? carnations? • Were there different results for different lengths Ask students: If you could change one thing about Maving Cln Up: of carnations or celery? this investigation to learn something new, what •Where do you think the water goes once it gets would you try? When we change one part of an r===Capillary J\cliian I to the top of the plant? experiment to see how it affects our results, this • What did you learn about water from this ex­ change is known as a variable. Use the chart in periment? be drawn through the your journal to record your ideas about what might Have you ever wondered • Did everyone have the same results? how water gets from narrow capillary tubes happen if you change some of the variables. Some •What did you like about this investigation? 1.eerning the roots of a tree to inside the plant. For possibilities are: • What variables did you try? capillary action to occur, • Will you get the same results if you use different Clbjecliive& its leaves, or why pa per •What surprised you? towels are able to soak up the attraction between the quantities of water or different amounts of food a soggy spill? It has to do water molecules and the coloring? Students will: The "Why" and The "How" tree molecules (adhesion) • What if you use soapy water? Salty water? Cold with the property of water This investigation illustrates the property of water known as capillary adion.
    [Show full text]
  • A Similarity Parameter for Capillary Flows
    INSTITUTE OF PHYSICS PUBLISHING JOURNAL OF PHYSICS D: APPLIED PHYSICS J. Phys. D: Appl. Phys. 36 (2003) 3156–3167 PII: S0022-3727(03)65940-1 A similarity parameter for capillary flows K A Polzin and E Y Choueiri Electric Propulsion and Plasma Dynamics Laboratory (EPPDyL), Mechanical and Aerospace Engineering Department, Princeton University, Princeton, New Jersey 08544, USA E-mail: [email protected] Received 11 July 2003 Published 25 November 2003 Online at stacks.iop.org/JPhysD/36/3156 Abstract A similarity parameter for quasi-steady fluid flows advancing into horizontal capillary channels is presented. This parameter can be interpreted as the ratio of the average fluid velocity in the capillary channel to a characteristic velocity of quasi-steady capillary flows. It allows collapsing a large data set of previously published and recent measurements spanning five orders of magnitude in the fluid velocity, 14 different fluids, and four different geometries onto a single curve and indicates the existence of a universal prescription for such flows. On timescales longer than the characteristic time it takes for the flow to become quasi-steady, the one-dimensional momentum equation leads to a non-dimensional relationship between the similarity parameter and the penetration depth that agrees well with most measurements. Departures from that prescription can be attributed to effects that are not accounted for in the one-dimensional theory. This work is dedicated to Professor Harvey Lam on the occasion of his recovery. 1. Introduction At the beginning of the 20th century, there were several works concerned with the dynamics of fluid penetration into 1.1.
    [Show full text]
  • Capillary Action Experiment
    Discover How Capillary Action Works with the Walking Water Experiment Grades: K - 6 Time: 20 mins set-up; view results the next day Lesson summary: Students will learn about capillary action. What’s the big idea? Outcomes or purpose: Water can move through tubes or Capillary action is when a liquid, like water, moves porous materials through a process up something solid, like a tube or into a material called capillary action. with a lot of small holes. Capillary action works against gravity because of three special forces called cohesion, adhesion and surface tension. Teacher background: Capillary action is at work in LGT gardens every day. Can you name three examples? If you said: Growboxes, peat pellets and GrowEase capillary mat trays then you’re right! You get bonus points if you said plants too! So what is capillary action anyway? Capillary action is when a liquid, like water, moves up something solid, like a tube or into a material with a lot of small holes. Capillary action works against gravity because of three special forces called cohesion, adhesion and surface tension. Tiny particles of water - also called water molecules - are like a group of good friends. They like to stick together! So when one molecule meets other molecules they tend to stick together in a group. This is called cohesion. Water molecules also like to stick to solid things. Sticking to solid things is called adhesion. Some examples of solids are: paper towels, the sides of a hollow tube or growing medium. All of these things have one thing in common: they are porous.
    [Show full text]
  • EXPERIMENT 17: CAPILLARY ACTION Objective: Can You Make Water Run Uphill?
    EXPERIMENT 17: CAPILLARY ACTION Objective: Can you make water run uphill? WHAT IS CAPILLARY ACTION? • Capillary action is defined as the movement of water within the spaces of a porous material due to the forces of adhesion, cohesion and surface tension. WHAT YOU NEED: STEP-BY-STEP: • Bowl of water 1. Add two or three • Food coloring drops of food • Two types of paper coloring to water. towels (choose ones with a different thickness or pattern) • Scissors, ruler 2. Cut paper towels 3. Place a mark on 4. Carefully dip into 2-inch each type of each paper wide strips. paper towel 1-inch from towel straight the edge. down into the bowl of water just to the 1-inch mark. 5. Slowly lift the towel straight up 6. Measure how far out of the water and watch the the water rises water “climb” onto the paper towel. after 30 seconds. QUESTIONS: • Which type of paper towel did the water rise higher on? Why? • What would happen if you put the paper towel into the water up to 2-inches? Do you think the water would rise more, less or the same amount? ? • Would the same thing happen with another liquid? Soda? Milk? • What other times might capillary action be used? For more fun and conservation-related experiments, please visit www.plumbingexcellence.org EXPERIMENT 17: CAPILLARY ACTION Instructor’s Guide ALIGNMENT WITH ILLINOIS STATE BOARD OF EDUCATION GOALS State Goal 11: State Goal 12: State Goal 13: Section A: 2a, 2b, 2c, Section D: 2b Section A: 2b and 2c 2d, 2e and 2f WHAT’S HAPPENING? Everyone knows that moving water tends to flow downhill and standing water will seek its own level, but sometimes water can run uphill.
    [Show full text]
  • Understanding Moisture Penetration of the Drainage Plane: a Building Science Principle Use Font Size 40
    Ralph E. Moon, Ph.D. CHMM, CIAQP, MRSA GHD, Inc. [email protected] 813-971-3882 (office) 813-257-0838 (desk) Understanding Moisture Penetration of the Drainage Plane: A Building Science Principle Use Font size 40 2019 IAQA Annual Meeting Why is this relevant to the IAQ professional? • The fundamental role of any structure is to keep water out • The drainage plane constitutes the elements that protect the structure and prevent water penetration • Understanding how water penetrates a structure is necessary to identify why indoor air quality may be impaired 2019 IAQA Annual Meeting The Four Basic Rules of Water Damage • Porous materials “suck” via capillary action • Water vapor diffuses from a higher to a lower concentration gradient • Air carrying water vapor moves from high pressure to low pressure • Water runs downhill by gravity 2019 IAQA Annual Meeting Basic Rule of Water Damage No. 4 • Water runs downhill by gravity (i.e., roof and flashing leaks, attic condensation) 2019 IAQA Annual Meeting Five Forces of Rainwater Penetration • Gravity • Surface Tension • Capillary Action • Momentum (Kinetic Energy) • Air Pressure Differences 2019 IAQA Annual Meeting 2019 IAQA Annual Meeting Which forces of drainage are shown? 1. Drip 2. Stucco Cracks 3. Failed vent cover Edge Gravity Surface Tension QUIZ Capillary Action Question Momentum (Kinetic Energy) Air Pressure Differences 2019 IAQA Annual Meeting Four Strategies to Prevent Rainwater Penetration: 4 D’s • Deflection • Drainage • Drying • Durable Materials Stave Church at Urnes, Norway – Norway’s oldest stave church that dates back to the early 12th century in its present form. Wooden components of an even older church were used to build it.
    [Show full text]
  • The Lost Chromatography Lesson
    The Lost Chromatography Lesson Background Science Summary: The wonder of the chromatography experiment, whether on earth or in the zero-G heavens, is the capillary action of water movement in the filter paper strip. Defying both the force of gravity and the expected inertia of the water drop, the liquid advances toward the ink spot. Were it not for an understanding of capillary action, the gradual soaking toward the ink spot would be a mystery. Indeed, what is the process which separates various colors of ink from the ink spot? How it works: Though ink appears to be composed of one color, it is not. Actually, several colored pigments blend together to make the color. This is known as a mixture. Artists mix different paint colors, like red, blue and green, to achieve a desired color which is a combination of the three basic colors. Because of chromatography, ink becomes soaked in the water drop so that the different pigments begin to separate or “bleed” apart. This results in the true colors being discovered. For Christa, one of the true colors in the black ink spot separated out by chromatography was blue. But more discussion of the science is needed. Why does one ink separate apart from another, and why does the water move toward the ink spot in the filter paper strip? The first question’s answer is: When a substance dissolves in a liquid like water, it is said to be soluble. Since this quality of being soluble differs in substances, each color will dissolve at a different rate so ink color separation occurs as the water soaks into the black ink spot.
    [Show full text]
  • Measurement of Dynamic Capillary Pressure and Viscosity Via the Multi-Sample Micro-Slit Rheometer Doyoung Moon, Kalman B
    ARTICLE IN PRESS Chemical Engineering Science ( ) -- Contents lists available at ScienceDirect Chemical Engineering Science journal homepage: www.elsevier.com/locate/ces Measurement of dynamic capillary pressure and viscosity via the multi-sample micro-slit rheometer Doyoung Moon, Kalman B. Migler ∗ Polymers Division, 100 Bureau Drive, NIST, Gaithersburg, Maryland 20899, USA ARTICLE INFO ABSTRACT Article history: We develop two direct methods to simultaneously measure the dynamic capillary pressure and the vis- Received 14 October 2008 cosity of fluids by application of differential forces during flow into micro-channels. In the first method, Received in revised form 30 January 2009 a series of external pressures is applied in conjunction with the dynamic capillary pressure and a “Bagley Accepted 18 February 2009 analysis” is applied to the flow front velocity, and in the second, we utilize differential gravitational forces. By explicitly measuring the dynamic capillary pressure, the measurement window of the recently Keywords: developed multi-sample micro-slit rheometer is extended to the regime where capillary forces are signif- Imbibition Dynamic capillary pressure icant. These measurement methods will be useful in understanding filling flows encountered in diverse Contact angle areas such as microfluidics, oil recovery and biological transport. Microfluidics Published by Elsevier Ltd. Viscosity Rheometer 1. Introduction (for the case of a rectangular capillary) where is the surface tension, is the contact angle, and h and w are
    [Show full text]
  • Investigation of Fluid Wicking Behavior in Micro- Channels and Porous Media by Direct Numerical Simulation
    Investigation of Fluid Wicking Behavior in Micro- Channels and Porous Media by Direct Numerical Simulation Ph.D. Dissertation By An Fu Master of Science, Mechanical Engineering, University of Cincinnati, 2016 April. 12th, 2019 Committee Chair: Dr. Milind Jog 1 A B S T R A C T Capillary transport phenomenon in porous media can be found in numerous natural processes and industrial applications such as inkjet printing, filtration, and enhanced oil recovery. In many of these applications, fluid penetration can be categorized into two orientations, viz., unidirectional (linear) and radial. In this thesis, the wicking process for both orientations has been investigated. Direct numerical simulation with Volume-of-Fluid (VOF) method has been implemented at pore scale (micro-scale simulation) to solve the transient governing equations, and global properties such as meniscus displacement and capillary pressure have been calculated to investigate macroscopic spontaneous imbibition (wicking) process. Linear fluid penetration kinetics within the viscous stage is first considered to investigate if the pore size and distribution are contributing factors which shift the overall kinetics, and how interface mobility is influenced by geometry layout. Various cases are simulated with increasing geometric complexity; from simple pore space connection with uniform pore size to complicated connection with random pore sizes. Average meniscus location, capillary pressure and other global variables are examined. The results indicate that within the viscous regime, the Lucas-Washburn behavior is valid in homogeneous porous media, even with random pore size and connection. However, in non-homogeneous structure, the wicking kinetic does depend on the distribution of pore spaces, and Lucas-Washburn equation fails to describe the linear macroscopic capillarity.
    [Show full text]