Dynamics of capillary transport in semi-solid channels

Capillary action has been described by Lucas and Washburn and extensively studied experimentally in hard materials, but few studies have examined capillary action in soft materials such as hydrogels. In tissue engineering, cells or dispersions must be often distributed within a hydrogel via microporous paths. Capillary action is one way to disperse such substances. Here, we examine the dynamics of capillary action in a model system of straight capillaries in two hydrogels. The channels had a circular cross-section in the micrometer size range (180–630 m m). The distance travelled over time was recorded and compared with the predictions of Lucas and Washburn. Besides water, we used a sucrose solution and a hydroxyethyl cellulose solution, both with viscosities slightly higher than that of water. The results showed that the distance travelled is proportional to the square root of time, x / ﬃﬃ t p , and that larger capillaries and lower viscosities result, as expected, in faster speeds. However, the absolute experimental values display large discrepancies from the predictions. We demonstrate that several possible reasons for these discrepancies can be ruled out and we describe a novel hypothesis for the cause of the retarded meniscus movement.


Introduction
Tissue engineering, prosthetic implantation, and controlled drug delivery are expanding fields where the control of artificial biomaterials is crucial for system success.2][3] The field has seen significant advances, such as the regulation of mesenchymal stem cell activity and proliferation by tuning the hydrogel's stress relaxation 4 or the injection of hyaluronic acid as a space filler. 5,6However, we lack a reliable method to disperse a concentration of soluble species, metabolites, therapeutics, or anti-fouling agents within a threedimensional network in the bulk. 7To address this deficiency, promising attempts have been made to transport substances within such three-dimensional networks.][8][9] Choi et al. and Cabodi et al. fabricated a microfluidic system made of a hydrogel for dispersing desired substances or for cell seeding and growth, with the flow driven by an externally applied pressure. 7,9Garg et al. used capillary action in a dried alginate gel as a matrix for cell growth, as they found capillary flow to be better in terms of cell seeding efficiency, cell dispersion, and cell viability than the previously described seeding approaches of injection, centrifugal, and orbital cultures. 8s we demonstrate here, the drying step is unnecessary, as capillary action in a templated hydrogel can be used to eliminate the need for the drying step.Capillary action has the advantages of being considerably faster and easier to use than external pressure-driven flow and not requiring extra equipment to distribute the liquid; in addition, less wetting liquid or cell solution is needed in the absence of a pump or the like.Besides its promising use in the above applications, capillary action is a phenomenon occurring in various other applications, such as wound care and hygiene products, printing on paper, oil recovery, and absorbent textiles.In all of these applications, it is important to be able to predict and tailor the mass transport in porous materials, resulting in desired properties.
The driving force for capillary action is the pressure difference created by the curved meniscus of a liquid evolving in a capillary, i.e. the Laplace pressure or capillary pressure: 10,11 Dp c ¼ 2g cos y r where g is the surface tension between liquid and gas, r the radius of the capillary, and y the contact angle.If y o 901 then cos y 4 0, so the liquid will penetrate the capillary.Or if the free surface energies are compared, penetration of the capillary occurs when

View Article Online View Journal | View Issue
This journal is © The Royal Society of Chemistry 2017 Soft Matter, 2017, 13, 2562--2570 | the surface energy of the dry wall, g SG , is greater than that of the wet wall, g SL (g SG 4 g SL ). 12 The opposing force is the viscous force known as the Poiseuille flow: where x is the travelled distance, Z the viscosity of the liquid, and : V the volume flow, which can be replaced with pr 2 Dx/Dt for a circular-cross-section channel of radius r and time t.By setting Dp c = Dp v while neglecting inertia and gravity, the Lucas-Washburn (LW) equation 10,11 is obtained: describing capillary action in a straight horizontal capillary of uniform diameter.As Fisher and Lark have demonstrated, the capillary must be sufficiently small, typically in the submillimetre range, for eqn (3) to be valid. 13The capillary length, l c ¼ ffiffiffiffiffiffiffiffiffiffiffiffiffi ffi g=ðrgÞ p , i.e. the length below which capillary forces are predominant over gravitational forces, is 2.7 mm for water, which explains the rough magnitude of the previous statement.2][23] The deviation of the experimental results from those predicted using the LW equation is often attributed to the exclusion of inertial forces, which, in particular, leads to such deviations in three-dimensional porous structures comprising channels of a range of diameters, such as foams, causing the flowing liquid to accelerate and decelerate throughout the system. 24Bosanquet, included inertial forces in eqn (3). 25 By balancing the acting forces and including the momentum of the liquid column, pr 2 rx : x, the equation of motion becomes: where r is the density, P an externally applied pressure, and : x the mass flow.On the left side of the equation are inertial forces and forces from the Poiseuille flow, and on the right side are the forces due to the externally applied pressure and capillary forces.Eqn (4), as modified by Bosanquet, becomes where the position x 2 4 x 1 and letting ][28][29][30] While capillary action has been studied extensively in many hard materials such as glass and silicon wafers, 13,17,18,21 to the authors' knowledge, no study describes the dynamics of capillary action in hydrophilic soft materials.However, it is observed that a drop spreads significantly more slowly on a soft substrate than on harder substrates 31 and forced wetting results in pinning effects on the contact line. 32These effects are caused by the local deformation of the substrate, which in turn leads to contact line deformation, termed the wetting ridge, which originates from the vertical component of surface tension in the wetting liquid. 32Swelling is another factor to be considered if the wetting liquid is soluble in the substrate, 33 which is enhanced at the three-phase contact line due to increased local flux in turn caused by increased evaporation in these areas. 348][9] It is a negatively charged polysaccharide derived from seaweed and consisting of repeating units of linearly (1-4)-linked b-D-mannuronic acid (M units) and a-L-guluronic acid (G units). 35In the presence of multivalent cations, such as Ca 2+ , alginate forms a gel linking only the G units and forming a network in the form of an ''egg box''model. 36,37The composition of M and G units depends, among other factors, on the plant species, season, and growth conditions. 35Another hydrogel commonly used in biomedical applications is agarose.It is the neutral fraction of the two major fractions in agar, which originates from red algae.Agarose powder is soluble in water at elevated temperatures and forms a gel upon cooling at 35-48 1C, with the methoxyl content influencing the gelling temperature. 37When the wateragarose mixture cools, agarose chains form a double helix.The properties of the gel depend on the cooling rate, since the system spatially optimizes its structure with increasing time.
The aim of this study was to determine the dynamics of capillary action in hydrogels in which a straight capillary was created.Capillaries with cross-sections of 180-630 mm were investigated and the distance travelled over time by a wetting liquid was recorded.Furthermore, we wanted to study the penetration behaviour of a sucrose solution having a high concentration of solutes and higher viscosity than water, generating osmotic pressure.The sucrose solution was also compared with a hydroxyethyl cellulose (HEC) solution having a low concentration and therefore a low osmotic pressure at a viscosity similar to that of the sucrose solution.Experiments were conducted using alginate and agarose gels.Note that in the following text, the terms wicking, imbibition, capillary action, and infiltration are used synonymously to refer to the surface-tension-driven flow.

Materials
Hydrogels were prepared using Protanal RF 6650 alginate containing 70% guluronate (kindly provided by FMC Biopolymer, Girvan, UK), Mikhart 2 calcium carbonate CaCO 3 (Provencale SA, Brignoles, France) with an average particle size of 10 mm, This journal is © The Royal Society of Chemistry 2017 Sigma-Aldrich, USA).The wetting liquids were Milli-Q water alone, with added hydroxyethyl cellulose (HEC; Hercules, Wilmington, USA), or with added sucrose (Nordic Sugar, Arlo ¨v, Sweden).All wetting liquids were stained with methylene blue (Fluka Chemie, Buchs, Switzerland) to increase the optical contrast to the material and ruthenium red (Sigma-Aldrich, Darmstadt, Germany) was used to stain alginate gel for analyzing microscopic cross-sections.

Gel preparation
An aqueous alginate solution containing 2% w/w alginate was prepared by adding alginate at room temperature (RT) under vigorous stirring.The dispersion was then heated to 80 1C for 30 min, after which it was left to cool to RT and the pH was adjusted to 7. Immediately before gel preparation, the alginate solution was deaerated under vacuum for approximately 10 min.Controlled gel formation was achieved by mixing CaCO 3 and GDL in water just before adding the alginate solution to yield a final alginate concentration of 1.7% w/w.The stoichiometric relationship of calcium ions to G units of alginate was 0.8.To form straight capillaries with circular cross-sections, immediately after mixing, the mixture was suctioned into 1 mL disposable syringes, each with a nylon string (of various diameters; nylon, monofilament) inside to create a channel with a circular crosssection.After standing for 24 h at RT, the gel was used for wicking experiments and subsequently discarded.The strings were gently removed before the use of the gel to leave a straight empty capillary.The actual diameter was measured by slicing the gel, staining it with ruthenium red and measuring the cross-section of the channel under a microscope (Axio Scope.A1; Zeiss, Oberkochen, Germany).The following diameters were determined: 180 (AE10), 250 (AE7), 410 (AE2), and 630 (AE29) mm.Agarose gels were prepared by boiling 2% (w/w) agarose on a hot plate for 5 min and subsequently replacing the evaporated water.The procedure continued in the same way as described above, resulting in gels upon cooling to RT. Channels of 250 mm were produced in agarose gels.

Wicking measurements and analysis
Wicking measurements were carried out by connecting the channel in a horizontal position to a drop of ultra-pure water stained with methylene blue (0.03% w/w) and recording the spontaneous infiltration using a camera (EX-ZR300; Casio, Tokyo, Japan) acquiring up to 120 frames per second.The video was separated into single frames using FFmpeg software (2013-07-09 git-00b1401) and the liquid was tracked using ImageJ (1.50b).

Surface tension
The surface tension of the water after adding methylene blue, sucrose, and HEC solutions was measured using an optical tensiometer (Attension Theta Lite; Biolin Scientific, Stockholm, Sweden) using the pendant drop technique.Due to the short durations of the wicking measurements, the surface tension was determined only within the first 3 s after a droplet was created.

Contact angle measurements
The material surfaces were characterized using a Fibro DAT1100 dynamic contact angle analyser (Fibro System AB, Ha ¨gersten, Sweden) after depositing a drop of pure water on the flat surface of the gel.Images were automatically acquired at 20 ms intervals directly after deposition of the drop.The mean contact angle of water on the gel surface was determined for at least five droplets.Staining the water with methylene blue was found not to influence the results.In addition, the dynamic contact angle of the moving meniscus was measured directly in a capillary.A highspeed camera (DCU223m; Thorlabs, Newton, NJ, USA) was connected to an objective (10Â Plan Achromat Objective, 0.25 NA; Olympus, Tokyo, Japan), serving as a microscope, which was pointed towards the channel at a specified distance from its entrance.The camera, operating at approximately 140 fps, captured four to eight images per experiment, depending on the speed of the meniscus.The speed at the point of measurement could be determined from the image series using the same method as described above for the wicking measurements.The contact angles were determined using ImageJ.

Rheology
Oscillatory rheology measurements were carried out on the gels using an MCR300 rheometer (Anton Paar Physica, Ostfildern-Scharnhausen, Germany) with a cone-plate geometry.The cone had a diameter of 50 mm and an angle of 11, and the gap between the cone and the plate was 52 mm.The temperature was controlled using a Peltier system and set to 20 1C.The strain was measured at a constant frequency of 1 Hz at shear stresses of 0.1-12 Pa.
The viscosity of the wetting liquids was measured using the same cone and plate and by applying rotational motion between shear rates of 1-600 s À1 and 600-1 s À1 .

Results and discussion
When a sufficiently small capillary is connected to a wetting liquid, the liquid will automatically penetrate the channel at a decreasing speed over time.Below, we will explain how our results are congruent with the LW predictions for y o 901, consider what is not congruent, and rule out a number of reasons for the deviations.Finally, we will demonstrate how our results can be interpreted in a novel way.

Infiltration speed in hydrogels follows the x 2 p t relationship
A droplet of pure water connected to a hydrogel penetrates a 6 cm-long micro-channel in the hydrogel in less than 3 s.Fig. 1 represents these dynamics for water in an alginate gel.The squared distance, x 2 , that the water travelled is plotted over time, t, through capillaries of 180, 250, 410, and 630 mm inner diameters.The data points represent the means of at least five experiments, while the shaded areas represent the standard deviation.The smaller the capillary diameter, the slower the wicking speed.This trend corresponds to the LW predictions given by eqn (3).In Fig. 1, one can see that the squared travelled distance is proportional to the time, x 2 p t, as also predicted by eqn (3).15][16][17][18][19][20]25 The gels were investigated for damages and the real diameter of the capillaries was determined by slicing the gels into cross-sections, as shown in the inset in Fig. 1.Using a 10Â magnification lens revealed no roughness left by the strings and in line with that very little force was necessary to remove the strings from the gel.Wicking experiments were also conducted using liquids with differing concentrations and thus osmotic pressures, but similar viscosities, namely, a sucrose solution (21.4% w/w) and a HEC solution (0.1% w/w); the properties of these liquids are presented in Table 1.Fig. 2 shows the slopes of the squared distance over time, k = Dx 2 /Dt, plotted over the viscosity, Z, of the wetting liquids and the capillary diameter, D, in alginate gels.The ball sizes represent the capillary diameters, which are also indicated in mm next to the balls; the ball colours show the viscosity, as indicated.It can be seen that low viscosities and large capillary diameters result in the steepest slopes, i.e. the capillary action is the fastest, while at high viscosities and small diameters the capillary action is the slowest, as the example of HEC illustrates.Because of the low concentration of HEC, the osmotic effect is assumed to be low, although the viscosity is similar to that of sucrose.The tendencies for higher viscosities to result in slower capillary action and larger diameters to result in faster capillary action are accordingly predicted by eqn (3).Our results also indicate that capillary action in empty capillaries interspersed within a hydrogel follows trends similar to those of the more commonly studied systems, such as capillaries in glass.Unlike the other systems, however, the mechanisms of capillary action considered here are based mainly on water-water interactions, because hydrogels consist predominantly of water (about 98%).The reason why capillary action occurs at all is that the water in the hydrogel has a tendency to minimize its surface energy, g SL o g SG , where g SL is the surface tension at the solid-liquid interface and g SG is the surface tension at the solid-air interface.By filling the capillary with water, the surface area exposed to air is minimized and so is the surface energy.

Absolute values do not indicate LW behaviour
The previous section describes the similarities between our experimental results and the LW predictions.If we insert the absolute values from our experiments into eqn (3), we obtain values deviating from the LW predictions.To quantify this deviation, we calculated the ratio of the experimental to the predicted slopes, k exp /k LW .The result is the relative difference between the experimental and predicted results, or how much slower our results are compared with LW predictions, which are plotted over the capillary diameters in Fig. 3.The closer the points in the graph to the LW results (the dashed line at k exp /k LW = 1), the   This journal is © The Royal Society of Chemistry 2017 more similar the experimental results to the predictions.The parameters used in eqn (3) are listed in Table 1.
The experiments with water and HEC conducted in alginate resulted in large discrepancies between the predicted and experimental results.The speed of liquid, v, and with it also the slope, k, were only 0.3-0.4times as high as predicted, though they were 0.6 times as high for the sucrose solution, specifically.Considering the high concentration of the sucrose solution (21.4%), we would expect an osmotic effect in which water from the gel flowed into the highly concentrated sucrose solution.The sucrose solution would become diluted, which in turn would reduce the viscosity with increasing x of the solution in the channel.Inserting a lower viscosity into eqn (3) would result in a higher k LW , and as k exp is determined directly from the experiment, it does not change.The outcome would be a larger discrepancy between predicted and experimental values even for sucrose.The results for sucrose will be discussed further in Section 3.5.
Table 2 presents the slopes for all tested capillary diameters and hydrogels with water as the wetting liquid.The results are derived by calculating the mean of the individual slopes for each tested capillary diameter.The standard deviation and number of replicates, n, are given in parentheses.We fitted the lines from 50% of t to ensure a linear relationship and exclude inertial effects.In a capillary of diameter 180 mm, the slope, k, in alginate was 1.2 Â 10 À3 compared with the LW prediction of 3.5 Â 10 À3 m 2 s À1 .In the tested capillary of the largest diameter, i.e. 630 mm, the difference was even greater, with an experimentally measured slope of 3.4 Â 10 À3 in alginate versus a predicted slope of 11.1 Â 10 À3 .To rule out the possibility that this effect was related to something specific to the calcium alginate gel, we also investigated capillary action in agarose gels.Agarose gels containing a capillary of diameter 250 mm were prepared as described in Section 2, ''Experimental''.As in alginate, in agarose with water as the wetting liquid, the experimental value of k was considerably lower than that predicted by the LW equation.If, on the other hand, the experimental value of k is used in eqn (3) and the contact angle is calculated, the result is a contact angle of 741 for water in 630 mm capillaries, decreasing as the capillary diameter decreases; for example, in 180 mm capillaries, the contact angle is 661.HEC results in the highest value, 861 for all tested diameters.These values are rather high, and the experimentally measured dynamic contact angles will be discussed below.

Possible reasons for the deviation between experimental and LW-predicted values
Agreement between capillary inflow and x / ffiffi t p is commonly reported.Agreement between experimentally obtained values for the slope, k, and the predicted values obtained using the LW equation is, however, rarer.As mentioned above, the absence of inertial forces in the LW equation, i.e. eqn (3), and the use of a static contact angle in the LW equation, while the contact angle in practice changes with speed, 22,26-30 are often cited as reasons for the discrepancy between the predicted and experimentally obtained values of k.However, neither the presence of inertial forces nor the use of a static contact angle in the LW equation can explain the deviation between the predicted and experimental values obtained here (Fig. 3).
Inertial flow was described by Bosanquet 25 and further developed by Ridgway et al. 24 as follows: The equation is valid in the absence of external pressure, as is the case here, and if aÁt { 1; for term a, see eqn (6).The distance travelled, x, is therefore directly proportional to the time, t.If we set t { 1/a, we get Using water as the wetting liquid in a capillary of 630 mm in diameter, we should consider inertial flow at t { 0.012 s, equivalent to a distance of 1.6 mm from the capillary entrance.
As we have calculated the slopes, k, based on much later times (50% of the entire penetration time, which is 0.39-0.79s in a capillary of 630 mm diameter and using water, equivalent to a distance of 35-50 mm from the entrance), we can therefore neglect inertial flow.Reducing the diameter of the capillary or increasing the density of the wicking fluid (as in the case of a sucrose solution) still does not increase the time at which inertial flow becomes important for the distance and time of interest here.The dynamic contact angle of water over alginate and agarose gels was measured experimentally on a flat surface of each gel using a Fibro DAT1100 dynamic contact angle analyser (Fibro System AB).After initially decreasing, the contact angle stabilized at approximately 151.This value should be considered an overestimation due to the inherent difficulty of measuring smaller contact angles (o101) using the Fibro device, leading to the overestimation of the mean contact angle.We therefore regard the contact angle between water and alginate or agarose gels to be o151.Since the aspect ratio of liquid volume to gel surface is different in small capillaries, we also measured the dynamic contact angles directly in the capillary during an experiment (see Fig. 4a).Fig. 4a shows a series of three images acquired just over approximately 2 cm from the capillary entrance with water coming from the left-hand side.The speed is plotted against the contact angle in Fig. 4b.The contact angle is largest, approximately 331, at higher speeds and declines to approximately 101 at rest.An increasing contact angle with speed is in agreement with reported values for wicking in glass tubes. 22,26It should be noted that the contact angle at rest obtained using this method agrees reasonably well with the value of o151 obtained using the Fibro DAT1100 dynamic contact angle analyser (Fibro System AB), given the inherent overestimation of the contact angle.
As shown in Fig. 4, the measured contact angles are much lower than the calculated values (see Section 3.2); it therefore follows that the dynamic contact angle does not solely explain the deviation between our experimental results and the LW predictions.The two commonly cited reasons for such discrepancies, i.e. the presence of inertial flow and the use of a static instead of dynamic contact angle, fail to explain the difference observed between predicted and experimentally obtained values in the systems studied here.
To rule out the effect of the difference between the water's chemical potential in the channel and in the hydrogel pore wall on our experimental results, we used water enriched with Ca 2+ ions as the wetting liquid at the same concentration as present in the alginate gel.Contrary to what was discussed above for the high-concentration sucrose solution, in the case of pure water as the wetting liquid, we assumed an osmotic flow of water from the flow channel into the gel wall.However, the results indicated no difference in the slope, k, when using calcium-enriched versus pure water.Further, swelling of the gel is not expected on the time-scale used in this study. 40nterstitial flow occurs when a liquid flows over a porous solid medium 41 or over hydrogels such as agarose 42 and calcium alginate. 43In all cases, the boundary layer of the fluid flow extends into the porous medium, inducing increased movement of the fluid within the gel.It follows that interstitial flow has the effect of causing a non-zero velocity at the boundary of the channel, 42 increasing the flow speed and decreasing the friction.Still, as energy is needed to accelerate the interstitial fluid, a lack of energy could slow down the meniscus.To estimate the importance of this effect, we assume that the interstitial fluid is accelerated uniformly with acceleration a for time T over a distance d = aT 2 /2 until it reaches the speed v I = aT.By assuming a uniform pressure distribution across the channel and gel, and using Darcy's law, it follows that da = v I 2 /2, where the interstitial velocity can be estimated to be v I E 8k : x/(fr 2 ), where k is the permeability of the gel and f its porosity.The accelerating force is thus given by where R is the radius of the gel cylinder.Because of the low hydraulic permeability of alginate gels (k E 10 À14 m 2 ), 44 and thus the very slow interstitial flow, inserting relevant values into eqn ( 9) and comparing them with the viscous force term from eqn (4), we found that F acc /F visc E 10 À9 .Therefore, the acceleration effect is very minor and not strong enough to cause the deviation seen in Fig. 3.

Proposal for the existence of a wetting ridge
A deformable substrate undergoes local deformation upon deposition of a water droplet, causing pinning effects, 32 and it has been demonstrated that the local deformation of the substrate slows the spreading speed, significantly increasing the spreading time of a droplet. 31To increase the spreading time by more than one order of magnitude, a nanometre-scale vertical displacement of the solid is enough. 31The effect has been modelled by Shanahan and Carre ´. 45 We expect the dynamics on flat surfaces to be similar to those of a moving meniscus in a circular capillary.However, to our knowledge, experimental wicking studies investigating deformation or wetting ridges in soft materials have yet to be conducted.Shanahan and Carre ´45 gave the following expression for the total braking force due to local deformation: where n and U 0 are fitting parameters depending on the material, m is the shear modulus of the material, and d is a cutoff length.Hence, the equivalent of eqn (4), ignoring inertial effects but taking deformations into account, is By solving this equation numerically using MATLAB (The MathWorks Inc., Natick, MA, USA), the impact of the viscoelastic braking on the capillary flow was estimated.The slopes of the function x(t) 2 were computed in the same way as for the experimental data.The parameter n was set to 0.5 and the value of K el = mdU n 0 was fitted to match the model to the data for water in Fig. 3.The value obtained was K el = 0.0027, which is not unreasonable compared to the value 0.40 for rubber, 45 which has a larger shear modulus than the gels used here.The lower value for K el indicates a higher rate of dissipation and a higher braking force on the meniscus in the gel.In Fig. 5, the output of the model is shown for water and for sucrose solution in comparison with the experimental data obtained from Fig. 3.The same value of K el was used for both liquids; only the surface tension and viscosity were changed.Still, the fit is fairly good, with the behaviour being correct for different diameters, which makes it plausible that the local deformations are responsible for slowing down the meniscus.The fit to the data for water could be made better by fitting the value of n, but since the number of data points is limited, it is considered sufficient to use one free parameter to prevent overfitting.
Shanahan and Carre ´also gave an equation for the dynamic contact angle as a function of the wetting front velocity during a spreading drop experiment, when the local deformations were assumed to be the cause of the change in contact angle. 45This equation does not fit our data in Fig. 4b, most likely because the velocities are orders of magnitude larger compared to the data for the spreading drop.This is reasonable considering that an increase in dynamic contact angle with increasing speed is also observed for hard materials like glass and thus that other effects may be responsible for the increase. 22,26e have ruled out numerous reasons for the deviation between k determined in our experiments and k predicted using the LW equation.We therefore hypothesize that the advance of the meniscus will be slower through a capillary in a soft substrate than through a capillary with a hard wall, for example, of glass, due to local deformation, i.e. a wetting ridge, along the moving contact line.The next section describes and discusses an indication for this hypothesis.

Deformation of the gel upon stress exertion
Small deformation rheology was used to determine the deformation (strain) that each gel exhibits upon exposure to the shear stresses relevant to this study (see Fig. 6).Both gels displayed a linear increase in strain with increasing shear stress, also see the work by Schuster et al. 46 for a more detailed mechanical and microscopic characterization of the alginate gels.The alginate gel was more sensitive to shear stress in that it deformed more at a given stress than the agarose gel.At the highest tested shear stress of 12 Pa, for example, the alginate gel displayed a deformation with strain of almost 0.9%, while the agarose gel deformed considerably less at 0.3%.The height of a wetting ridge is related to the elastic modulus of the substrate, 47 scaling as h E g sin y/Y, where Y is the Young's modulus. 45For the alginate gel, this gives an estimated height of the wetting ridge in the order of 1 mm, which is fairly large but still not visible in the images in Fig. 4a.Since the alginate gel was more prone to deformation than the agarose gel at a given shear stress, we assume a more pronounced wetting ridge in the alginate gel during capillary infiltration.This assumption is the object of further investigation.
Based on the greater deformation at a given shear stress in alginate than in agarose gels, we compared the speed of infiltration in the two gels (see Fig. 7).The squared travelled distance, x 2 , through a capillary of 250 mm diameter is plotted over time.The shaded areas represent the standard deviation.The plot shows, first, that the results are similar for the two hydrogels, independent of whether they consist of agarose or calcium alginate.Generally, as described above, the speed of the meniscus advance in both hydrogels deviates substantially from the absolute values calculated using the LW equation.It can be seen in Fig. 7 that before 0.25 s the infiltration speed is  very similar in both gels; after that, it starts deviating with the speed being faster in agarose than in alginate, despite the capillaries being of the same diameter.Eqn (10) demonstrates mathematically that the higher the shear modulus, the lower the viscoelastic braking due to the deformation of the substrate. 45lthough further investigation is needed, this strengthens our hypothesis that the observed deviation of k between the experimental and LW-predicted values is due to local deformation of the substrate.

Conclusions
We have demonstrated that it is possible to make use of capillary action in hydrogels despite their consisting predominantly of water.In capillaries of 180-630 mm in diameter, we also investigated more viscous wetting liquids, all of which follow the dynamics of the LW relationship x 2 p t. Also congruent with LW predictions, smaller capillary diameters and higher wetting liquid viscosities result in slower penetration speeds.The absolute experimentally determined meniscus speed was substantially lower than the absolute predicted speed.We hypothesize that this delay is caused by deformation near the three-phase contact line, i.e. the wetting ridge, of the moving meniscus.The viscoelastic energy dissipation of the substrate slows the penetration speed, leading to the deviation from the LW-predicted behaviour.This study also suggests that capillary action can be used to disperse a concentration of soluble species, metabolites, therapeutics, or anti-fouling agents through a spatial network of micro-channels in a hydrogel consisting predominantly of water.

Fig. 1
Fig. 1 Mean values at each time point of the squared distance, x 2 , that water travelled in alginate gels over time, t, through capillaries of 180, 250, 410, and 630 mm inner diameters; the shaded areas represent the standard deviation (n Z 5) for each diameter.The inset shows a micrograph of a cross-section of the alginate channel displaying the measured diameter for one sample.

Fig. 2
Fig.2The slopes, k, of the penetration speed over the viscosity, Z, of the wetting liquids (i.e.water, sucrose, and HEC solution) and the capillary diameter, D, in alginate gels.The colours and sizes of the depicted balls indicate different liquid viscosities and capillary diameters, respectively; the diameter in mm is also indicated next to each point.

Fig. 3
Fig. 3 Slope ratio, k exp /k LW , over the diameter, D, for the indicated liquids in alginate gels.The shaded areas show the standard deviation (n Z 5) at each point.The dashed line shows the LW-predicted behaviour.

Fig. 4
Fig. 4 Dynamic contact angles measured directly in the capillary: (a) a series of images of an alginate capillary of 250 mm in diameter; (b) the dynamic contact angle over the speed.

Fig. 5
Fig.5Comparison of the output of model(11) with the experimental slope ratio, k exp /k LW , for water and sucrose solutions.The model was fitted to the water data and used to predict the sucrose data.

Fig. 6
Fig. 6 Deformation of agarose and alginate gels at shear stress values experienced by the gels during capillary action, (n = 1).

Table 1
Properties of the hydrogels and the wetting liquids (i.e.water, sucrose solution, and HEC solution)

Table 2
Slopes, k = Dx 2 /Dt, of water in alginate and agarose gels compared with the LW-predicted values a This journal is © The Royal Society of Chemistry 2017 Soft Matter, 2017, 13, 2562--2570 | 2567