Research ArticleENGINEERING

Poro-elasto-capillary wicking of cellulose sponges

See allHide authors and affiliations

Science Advances  30 Mar 2018:
Vol. 4, no. 3, eaao7051
DOI: 10.1126/sciadv.aao7051


We mundanely observe cellulose (kitchen) sponges swell while absorbing water. Fluid flows in deformable porous media, such as soils and hydrogels, are classically described on the basis of the theories of Darcy and poroelasticity, where the expansion of media arises due to increased pore pressure. However, the situation is qualitatively different in cellulosic porous materials like sponges because the pore expansion is driven by wetting of the surrounding cellulose walls rather than by increase of the internal pore pressure. We address a seemingly so simple but hitherto unanswered question of how fast water wicks into the swelling sponge. Our experiments uncover a power law of the wicking height versus time distinct from that for nonswelling materials. The observation using environmental scanning electron microscopy reveals the coalescence of microscale wall pores with wetting, which allows us to build a mathematical model for pore size evolution and the consequent wicking dynamics. Our study sheds light on the physics of water absorption in hygroscopically responsive multiscale porous materials, which have far more implications than everyday activities (for example, cleaning, writing, and painting) carried out with cellulosic materials (paper and sponge), including absorbent hygiene products, biomedical cell cultures, building safety, and cooking.


As a major constituent of plants, cellulose has been used as a source of energy (1), food (2, 3), building materials (4, 5), clothing (6), and hygiene products (7) throughout human history. In particular, transfer and preservation of information has relied on wetting of cellulosic materials for millennia, beginning with papyrus in 3000 BCE. Porous materials made of cellulose still abound around us as paper and sponges, among many others (8, 9). When we bring a dry paper or sponge into contact with water or ink, it absorbs the liquid while swelling simultaneously. Cellulose is a polymer whose chains are linked via a hydrogen bonding, and water molecules participate in the binding sites, causing the polymer volume to increase. Such physicochemical interaction of water and porous structure is called hygroscopic expansion, a mundane process observed in cleaning (10), painting (11), and writing (12).

Capillary imbibition of liquid in porous media (1215) is described by Darcy’s law, which gives the flow rate, q, as a function of the permeability k and gradient of driving pressure ∇p: q = −(k/μ)∇p, where μ is the liquid viscosity. Because the permeability and the driving capillary pressure are respectively scaled as the cross-sectional area of the fluid conduit and the meniscus curvature at the wetting front, both of them are determined by the pore size. When the porous structure deforms due to liquid infiltration, the poroelastic theory gives the pore size and the flow rate in general. The theory is built upon the basic assumption that the pore pressure increases with the water content, the amount of liquid inside the void, which in turn causes the pore to swell (16). It has successfully described the behavior of liquids in many porous media, such as soil (17), sandstone (18), and hydrogel (19). However, the present problem of wicking and hygroexpansion defies such classical theoretical understanding for the following reason. As water progressively wets the cellulose materials, the macroscale pores at the spreading front expand due to swelling of the surrounding walls or scaffolds, not by increased pore pressure. The pore pressure should rather decrease because of volume expansion. This disobeys the fundamental framework of poroelastic theory, and thus, a completely different approach should be devised to understand the hygroscopic expansion of cellulosic porous materials containing macro voids.


Characteristics of cellulose sponges

As a model system to study this problem, we bring a commercial cellulose sponge (VWR) into contact with water or various liquids (table S1) and observe the liquid front rise against gravity (Fig. 1A). When using aqueous liquids, the sponge swells while being wetted. Plotting the rise height h versus time t (Fig. 1B) reveals that the power laws of h versus t differ in the early (filled symbols) and late (empty symbols) stages. Here, the heights are scaled by hJ, Jurin’s height (20), a characteristic rise height at which the gravitational and capillary forces for a macropore are balanced: hJ = γ/(ρgR), with γ, ρ, and g being the liquid-air surface tension coefficient, the liquid density, and the gravitational acceleration, respectively. Figure 1C shows that the transition height at which the power law changes from h ~ t1/2 (filled symbols) to h ~ t1/5 (empty symbols) corresponds to Jurin’s height of macro voids. Rationalizing these power laws allows us to understand the fundamental wicking dynamics of the hygroexpansive, heterogeneous porous materials.

Fig. 1 Capillary rise in cellulose sponges.

(A) Optical images for wicking of water (main panel) and turpentine (inset) in the initially dry sponge. The sponge swells when contacting water or aqueous liquids (movie S1). From left to right, t = 0, 1, 10, 100, 1000 s. Scale bar, 10 mm. (B) Experimentally measured rise height of water (black symbols) and turpentine (red symbols) versus time. In the early stages (filled symbols), the rise height grows like t1/2 (gray line) for both the liquids. In the late stages (open symbols), the rise height of water follows the t1/5 rule (black line), whereas that of turpentine behaves like t1/4 (red line). (C) The power law of the height changes when the rise height h reaches Jurin’s height of macropores: hJ = γ/(ρgR), so the transition occurs at h/hJ ~ 1. The symbols for different liquids are listed in table S1. All the experimental data for the rise height are the average of three measurements, with the error bars smaller than the size of symbols.

We begin with characterizing the pore structure of the cellulose sponge. As shown in the scanning electron microscopy (SEM) images (Fig. 2, A to C), it consists of numerous cellulose sheets with two-dimensional microscale pores surrounding macro voids (13). The sheets approximately 10 nm in thickness are randomly stacked with nanometric spacings. Measuring the size distribution of the pores, we find the average radii of macro and micro voids to be R = 0.73 mm and r0 = 4 μm, respectively (fig. S1). Pores finer than the micrometric voids are hardly found in the sheets.

Fig. 2 Microscopic images of the cellulose sponge.

SEM images of macropores (A), micropores (B), and cross section of the sheets (C). Scale bars, 300 μm (A) and 10 μm (B and C). (D) Merging of micropores due to hygroscopic expansion of the cellulose sheet, as imaged by ESEM. The pores start to grow when the relative humidity (RH) exceeds 90%, and they coalesce with their neighbors. Scale bar, 10 μm. (E) Schematic illustration of micropore expansion and coalescence. The micropores grow from r0 to r1 in radius (pore expansion) and then merge to form large micropores of radius r (coalescence).

Capillary flows and volumetric expansion in porous media

Darcy’s law gives the wicking velocity u in a porous medium, which we now write as u = −(k/μ)dp/dz. The driving pressure arises as a consequence of capillary action so that Δp ~ γ/λ, where ~ signifies “is scaled as” and λ is the radius of curvature of the front meniscus. The permeability k is scaled as the cross-sectional area of fluid conduit, over which the viscous stress resisting the fluid flow develops. Noting that h measures the distance from the free surface of the liquid reservoir to the wet front, u and ḣ = dh/dt can differ when the media volume changes. For the isotropically expanding materials like cellulose sponges as shown in Fig. 3D, the bottom of the sponge descends by hεs, with εs being the hygroscopic strain of the saturated sponge, so that the total wet distance Hh(1 + εs). Not all the liquid flowing into the sponge contributes to the rise of H because of the transverse expansion of the sponge, leading us to write Ḣu/(1 + εs)2 (section S1). The hygroscopic strain of the used sponge is at most 0.23, allowing us to neglect higher-order terms of εs. Because hH/(1 + εs), we get ḣu/ζ, with ζ ≈ 1 + 3εs being the coefficient of volumetric expansion.

Fig. 3 The difference of wicking behavior between early and late stages.

(A) Optical image of a sponge wetted by water. At small h, or in the early stages, the macro voids are completely filled with liquid. However, at large h, or in the late stages, they are only partially filled with liquid. A white curved line indicates the deformation of sponge due to the constraint of the dry upper part (section S1). The blue and red boxes show macro voids completely and partially filled with liquid, respectively. Scale bar, 10 mm. (B) Schematic of liquid-filling behavior in the late stages. Macro voids are not completely filled with liquid due to gravitational effects, whereas micropores are fully occupied with liquid. The radius of curvature δ of the meniscus in the macro void is determined by the balance between gravitational and capillary forces: δ ~ γ/(ρgh). (C) Schematic of a macro void of radius R in a microporous sheet in the early stages. Both macro and micro voids are completely filled with liquid. (D) Simplified model of the sponge whose wetting behavior is mathematically analyzed. Black box shows the liquid path near the wetting front. The liquid permeates into micropores from the wet corner of macropores.

When the macro voids are completely filled with infiltrating liquid, as observed for the early stages of capillary rise (Fig. 3, A and C), λ ~ R and k ~ R2. Then, we get ḣ ~ γR/(ζμh) so that h ~ [γRt/(ζμ)]1/2, which is consistent with Lucas-Washburn’s rule (21, 22) except the fact that the effect of swelling (ζ) is included. We plot the rise heights of Fig. 4A according to our scaling law to find the scattered data to collapse onto a single straight line in Fig. 4B together with the data of nonaqueous liquids (εs = 0). When plotting h versus γRt/μ without ζ in Fig. 4C, we find two distinct lines depending on whether the sponge swells or not.

Fig. 4 Wicking dynamics in different stages.

(A) Experimentally measured rise height versus time in the early stages. (B) The early-stage data are collapsed onto a single line when plotted according to our scaling law h ~ [γRt/(ζμ)]1/2. The law holds for both the aqueous and nonaqueous liquids. (C) Two collapsed lines for aqueous (black line) and nonaqueous liquids (red line). (D) Experimentally measured rise height of aqueous liquids versus time in the late stages. (E) The late-stage data are collapsed onto a single line when plotted against our scaling law h ~ (Bt)1/5. (F) The experimental results appear to collapse onto a single line but disobey the t1/4 rule (Eq. 1). The symbols for different liquids are listed in table S1.

Darcy’s law for late stages

Beyond Jurin’s height, macro voids cannot be completely filled because the capillary pressure based on the void radius cannot withstand the hydrostatic pressure. Then, the foregoing model fails (13), which is consistent with the change of the slopes in Fig. 1B. In this regime, the rise is rather driven by the capillary pressure provided by the micropores of characteristic radius r so that λ ~ r. However, the liquid does not flow only through microporous sheets, but it can also wet the corners of macro voids in such a way that the radius of corner meniscus δ balances capillary and hydrostatic pressure: γ/δ ~ ρgh. We display the image of the macro void partially filled with liquid and its schematic in Fig. 3 (A and B). Because the liquid that advances the wetting front can be supplied from the wet corners rather than the network of micropores owing to reduced viscous stress (for δ >> r) as shown in the box of Fig. 3D, we should take the permeability k as the cross-sectional area of the wet corner. A simple geometric consideration allows us to write the area as δ2 so that k ~ δ2. Then, Darcy’s law, u ~ δ2γ/(μrh), givesuγ3μ(ρg)2rh3(1)

Pore growth due to hygroscopic swelling

When the micropore size is invariant as r = r0, we easily get h ~ t1/4, which was shown to hold for the rise of nonaqueous liquids without causing hygroscopic swelling (13). However, in case of water wicking, we have discovered drastic shape changes of micropores in an environmental SEM (ESEM) chamber, where the RH around a sponge specimen has increased from 10 to 100%. See Fig. 2D and movie S2 for experimental images and Materials and Methods for experimental procedure. The expansion process of the porous sheet can be decomposed into two steps (Fig. 2E). First, the pore size grows from r0 to r1, allowing us to write r1 ~ r0(1 + ε). As the expanded pores get closer, they merge to form large pores of radius r accompanied by the decrease of the number of pores from N0 to N. Because of the insignificant change of the total area of pores upon coalescence, we get N0r12 ~ Nr2. Similar growth and coalescence of pores can be easily observed by stretching a macroscopic perforated polymer film as well as in microscopic ductile fracture (23) and polymer crack (24). The ratio N0/N > 1 should increase with ε, which we simply estimate as N0/N ~ 1 + βε (section S2 and figs. S1 to S4). The prefactor β is a function of a distribution of interpore distances, and a detailed discussion of the estimate is given in section S2. With the measured values of β ≈ 50 and ε ~ εs ≈ 0.23 in the cellulose sponges, we obtain N0/N ~ βε, and thus, r ~ (βε)1/2r0. When the wet sponges are dried, a similar size distribution of micropores to the originally dry state is recovered (section S3 and fig. S5).

Hygroscopic swelling in cellulose sheets

The degree of hygroscopic expansion of a porous sheet as the wetting front propagates is related to the amount of water absorbed in the sheet. See fig. S6 (section S4) for the schematic of the porous cellulose sheet being wetted. The hygroscopic strain ε = αη, where α is the hygroscopic swelling coefficient measured to be α ≈ 0.33 (section S5 and fig. S7). The volume fraction of aqueous liquid, η, in the front sheet is given by η = Vl/Vc, with Vl and Vc being the volume of liquid and of cellulose, respectively. Because a liquid is absorbed into a sheet of thickness s by the distance ld, Vl ~ N0r0lds. The diffusion length ld is scaled as ld ~ (Dτ)1/2, with τ being the characteristic time taken for the rising liquid to pass the sheet: τ ~ s/u. The diffusivity of a dense hygroscopic medium (19), D, is given by D = 2Gkc(1 − ν)/(1 − 2ν)/μ, where the shear modulus G = 1.62 MPa (section S6 and fig. S8), the permeability kc ~ rc2/32, with rc being the typical pore radius (order of 1 nm) of a cellulose sheet (25), and Poisson’s ratio ν = 0.3 (26). The cellulose volume Vc ~ N0r02s(1 − φ)/φ, where the porosity φ is the ratio of the micropore volume to the total sheet volume (section S4 and fig. S6).

Although the moisture diffusion at the wetting front determines the size of newly wetted pores, the diffusion length itself is insignificant as compared with the overall rise height. Namely, we find that the typical increment of the diffusion length Δld is much smaller than that of the rise height Δh, to give Δldh ~ 0.01 for a given duration even in the late stages.

Dynamics of wicking and swelling in late stages

The foregoing considerations allow us to write ε as ε ~ α(Ds/u)1/2φ/(1−φ)/r0, which leads tor(αβφr01φ)1/2(Ds)1/4u1/4(2)

The relation indicates that the micropore radius at the wet front increases as the liquid rising velocity (u) decreases, which results in the decrease in the capillary pressure (~γ/r). Combining Eqs. 1 and 2 and recalling u/ζ, we obtain a power law for the rise height in the late stages, h ~ (Bt)1/5, with B given byB=1ζ[γ3(1φ)1/2μ(ρg)2(αβφr0)1/2(Ds)1/4]4/3(3)

We plot the experimental results in the late stages of Fig. 4D according to scaling law (Eq. 3) in Fig. 4E. We see that the experimental data for various liquids are collapsed onto a single master curve despite the variations of γ, μ, ρ, and D, consistent with our theory. Figure 4F plots h based on the scaling law suggested for nonswelling sponges (13): h ~ {γ3t/[μ(ρg)2r0]}1/4. Although the data appear to collapse onto a single line, h does not follow h ~ t1/4, invalidating the nonswelling theory for the current situation. We show in section S7 (fig. S9) that even with aqueous liquids, the rise height follows the t1/4 law in the late stages for the sponges that have been pre-wetted and fully swollen in advance.


We have constructed the theoretical framework to analyze the liquid flow in hygroexpansive porous medium, where multiscale pore dynamics cannot be adequately accounted for by conventional poroelasticity theory. In particular, our analysis correctly captures the t1/5 behavior of the rise height in the late stages by considering the expansion of micropores as a function of the liquid rise speed, which determines the diffusion length within a cellulose sheet. Although the evaporation of water from the wet sponge occurs in the course of capillary rise, it has been ignored in our analysis because its rate (q ≈ 6 × 10−9 m/s as experimentally measured) is negligibly small compared with the typical wicking velocity u ~ 10−4 m/s in the late stages. We note that although cellulose is a major constituent of wood, the porous structure of artificially manufactured cellulose sponges is distinguished from wood that has grown in nature. One of those remarkable differences is that wood has submicrometric pores in a hierarchical order (25), which are invisible in the cellulose sponges.

Besides cellulose sponges, a variety of mundane and industrial materials of heterogeneous porosity can benefit from our theory, including biomedical devices such as the cytosponge (27), cell cultures (28), shape-morphing microneedles (29), soft actuators (30, 31), and plant seeds (32). Paper itself involves only microscale pores (33, 34), but crumpled or folded paper forms macro voids, wetting of which should be in line with the current problem. In addition to cellulosic porous materials, hygroexpansive porous bread composed of starch (35) was found to exhibit micropore coalescence and to follow the t1/5 law in the late stages of vertical wicking (section S8 and fig. S10). This enables us to start thinking about applying our theoretical framework to a relevant field in the science of cooking. Although we have concentrated on the power-law behavior of liquid dynamics, hygroscopic deformation dynamics of the heterogeneous porous solid structure would help us to fully understand and control the soft materials of ever-growing importance.


Experimental procedures of ESEM

In the ESEM chamber (XL-30 FEG, Philips), we placed a piece of cellulose sponge, 5 × 5 × 1 mm3 in volume, on a Peltier plate, which was kept at 2°C. Water vapor was supplied into the chamber to increase the RH. The sponge underwent drastic shape change as the environmental humidity reached approximately 90% by absorbing water molecules. The deformation of the sponge was insignificant when RH was below 90%, indicating that water molecules hardly infiltrate the sponge until the vapor pressure reaches a critical value. The imaging results are shown in Fig. 2D.


Supplementary material for this article is available at

section S1. Effects of isotropic volumetric expansion on liquid rise height

section S2. Correlation between hygroscopic strain and pore coalescence

section S3. Recovery of microporous structure of cellulose sponges

section S4. The volume fraction of aqueous liquid in a cellulose sheet

section S5. Effects of water concentration on hygroscopic strain

section S6. Mechanical properties of cellulose sponges

section S7. Capillary rise in pre-swollen sponges

section S8. Scaling laws of water rise within bread made from starch

fig. S1. The measurement data of the cellulose sponge structure.

fig. S2. Macroscopic experiments for pore coalescence.

fig. S3. Numerical analysis of porous sheet deformation.

fig. S4. Moisture flux into cellulose sheet in ESEM chamber.

fig. S5. Microporous structure of cellulose sponges after cycles of wetting and drying with water.

fig. S6. Analysis of the cellulose sheets.

fig. S7. Hygroscopic strain of saturated sponge for different water contents in aqueous glycerin and ethylene glycol.

fig. S8. Shear modulus of dry and wet cellulose sponge.

fig. S9. Capillary rise height of water versus time in an initially dry sponge (black) and pre-swollen sponge (red).

fig. S10. Capillary rise in porous bread.

table S1. List of liquid properties and symbols.

movie S1. Wicking and swelling in the cellulose sponge.

movie S2. The merging of micropores in the cellulose sheets.

References (3638)

This is an open-access article distributed under the terms of the Creative Commons Attribution-NonCommercial license, which permits use, distribution, and reproduction in any medium, so long as the resultant use is not for commercial advantage and provided the original work is properly cited.


Acknowledgments: We are grateful to L. Mahadevan for stimulating discussion. Funding: This work was supported by Samsung Research Funding and Incubation Center of Samsung Electronics (project no. SRFC-MA1301-05) and National Research Foundation of Korea (grant nos. 2016901290 and 2016913167) via SNU-IAMD. Author contributions: J.H. and J.K. carried out the experiments. J.H. and Y.J. analyzed the pore coalescence process. G.Y. and D.-N.K. performed numerical simulations. J.H. and H.-Y.K. developed the mathematical model and wrote the paper. H.-Y.K. conceived and supervised the project. Competing interests: The authors declare that they have no competing interests. Data and materials availability: All data needed to evaluate the conclusions in the paper are present in the paper and/or the Supplementary Materials. Additional data related to this paper may be requested from the authors.

Stay Connected to Science Advances

Navigate This Article