University of Birmingham Measuring the impact of channel length on liquid flow through an ideal Plateau border and node system

The phenomenon of foam drainage is a complex multi-scale process that unites molecular level interactions with bulk foam characteristics. Foam drainage is primarily governed by the flow of liquid in the channels and junctions that form between bubbles, which are known as Plateau borders (PBs) and nodes respectively. Existing theoretical work predicts the surface rheology of the PB and node air–liquid interface to influence liquid flow rates; however, direct experimental observations of this phenomenon remain scarce. This study recognises the clear need for a reproducible, accurate and standardised approach to directly studying liquid flow at the scale of a theoretically ‘ideal’ PB–node architecture. Measurements of PB geometric profiles and their apparent surface shear viscosities, m s , were made for an aqueous solution of Sodium Dodecyl Sulphate (SDS) at varying PB lengths, l 1 , and liquid flow rates in the range 10 m l min (cid:2) 1 r Q r 200 m l min (cid:2) 1 . Geometric profiles displayed previously unobserved transitions between PB relaxation and expansion towards the node, with expansion dominating under conditions approaching conventional foam drainage. Average values of m s in the PB relaxation regions showed virtually inviscid behaviour, with magnitudes of 10 (cid:2) 8 g s (cid:2) 1 o m s o 10 (cid:2) 4 g s (cid:2) 1 for l 1 in the range 27.5 mm \ l 1 \ 8.0 mm. Decreasing magnitudes of m s and degrees of shear thinning were observed with increasing l 1 . This was attributed to a compressibility of the interface that was limited by an SDS concentration dependence on l 1 . Numerical evaluation predicted the appearance of Marangoni forces that scaled strongly with liquid shear rates, and could therefore have been responsible for the apparent shear thinning behaviour.


Background
In general terms, a foam can be described as a large number of gas bubbles that are closely packed together. 3 Owing to the competing action of capillary and gravitational forces, the liquid fraction gradually decreases over time as liquid travels through the network of channels between bubbles and back into the bulk solution beneath. This phenomenon of drainage is a complex physicochemical process that ultimately leads to the collapse of the foams due to eventual film rupture and subsequent bubble coalescence. 33,42 The complexity of foam systems have necessitated a bottom-up approach from those attempting to model macroscopic drainage. By simplifying the extensive liquid network into its component channels, known as Plateau borders (PBs), and junctions or 'Nodes', the microscale building-blocks of the foam system can be better characterised and subsequently reassembled. A large body of theoretical work exists to describe PB flow; 20,24,28 however, corresponding experimental verification is limited due to the difficulties in observing and measuring such PB-node architectures.
Perhaps the most detailed direct observation of liquid flow through an individual PB within a bulk foam comes from Koehler, et al., 21 who used confocal microscopy to measure flow profiles in individual channels and a node. Since then, a greater focus has been given to forced drainage experiments through isolated PB systems where experimental and formulation parameters are easier to control and define. 9,14,[29][30][31] Much of the current theoretical work on isolated PB systems concerns the measurement and characterisation of interfacial properties, whose molecular origins have significant macroscopic effects on liquid flow and thus foam drainage. 1,36 Numerous experimental techniques are available to quantify the surface rheology of aqueous solutions, 11,13,35,46 however none of these account for the complex geometries and stresses of a PB-node system. 16,34,38 Development of such in-situ techniques is therefore paramount for both the verification and improvement of current microscale drainage theory. 7,22,23,29,43 The surface parameter that has seen the most in situ experimental investigation in foam systems to date is an apparent surface shear viscosity, m s , that ultimately relates the shear force per unit length of interface to an applied shear rate, _ g s . 44 While the true surface shear viscosity is a material function of the surfactant solution, separation of this value from dilational viscous effects and Marangoni stresses has been shown to be highly troublesome even for analysis of simple planar interfaces. The combination of these effects into an apparent surface viscosity is however, sufficient for analysing the combined effects of the interface on liquid flow tangential to the PB vertical axis. Other publications have featured this parameter in its dimensionless form, the Boussinesq Number, B 0 , which is scaled by the bulk liquid viscosity, m, and the PB's radius of curvature, R, according to eqn (1): PB systems with a low m s impart minimal shear on liquid flow, producing a more uniform flow profile resembling 'plug flow'. A high m s would produce profiles closer to that of the Poiseuille flow observed in rigid pipes. 37 The contrast between low and high m s is thought to be responsible for the two macroscopic regimes of node-dominated and channel-dominated drainage respectively, which describe whether the bulk of flow dissipation is thought to occur at the nodes or within the PBs themselves. 29,36 The 'Plateau Border Apparatus' of Pitois, et al. 29 (Fig. 1A) has provided a wealth of information on this topic, with identification of 'mobile' and 'rigid' interfacial behaviour depending on the choice of surfactant. By measuring pressure variations within the PB, the authors were subsequently able to determine the apparent PB surface viscosities and resultant hydrodynamic resistances that inform models describing macroscopic foam permeability. 28,30 However, despite average values of these resistances being in agreement with theory, the observed variation with liquid flow rate differed significantly. 2,25,28,30 Despite its successes, there are still shortfalls of the Plateau Border Apparatus that need to be addressed. Firstly, for channels laden with high mobility surfactants, the node is known to dominate flow dissipation, and therefore must be included in investigations. Pitois, et al. 31 made one such study in which the Plateau Border Apparatus was adapted to include a bubble at the base of the PB (Fig. 1B); however, the limited amount of data published was insufficient to account for large discrepancies in the current theory predicting hydrodynamic resistance of the node. 4,19,27,37 In addition, the adaptation of the setup meant that the apparent surface viscosity needed to calculate hydrodynamic resistance was taken from measurements of the original PB-only system. This neglected both the potential impact of PB length variation and of variations in PB surface rheology that may have been introduced by the very presence of a node. These factors cannot be assumed and must therefore be experimentally verified.
Further challenges of such a system include the complex process of pressure measurement, where the significant timescales to achieve equilibrium 29 and induced pressure fluctuations often compromise PB and node stability. Overall, this can make the measurement procedure extremely labour intensive and therefore severely limit the amount of data that can be collected.
In order to begin to address the challenges above, we propose the development of a purpose-built experimental setup to measure the apparent surface viscosity of 'ideal' PB-node systems for a wide range of flow rates and PB lengths. This reproducible and accurate technique will add to the limited body of experimental data that exists to describe such systems, while helping to standardise the approach to studying them. Overall, it is hoped that this will improve the comparability of findings between both previous and future studies, thus enabling further progress to be made in understanding these systems.

Theory
The bulk of existing theory relating PB surface viscosity to PB geometry, pertains to idealised 'infinite' PBs. Owing to the complexity of PB geometries, these simplified systems help to identify cross-sectional areas, S, such as that shown in Fig. 2, where S is a simple function of film thickness, w, and PB Radius, R, described by eqn (2): 20 A common approximation is applied to cases where films are thin (w { R). In these cases, S can be approximated as a function of R only (eqn (3)) using the geometric factor c E 0.161. 20 Liquid flow through a PB is defined within the regimes of viscous and inertial flow, with viscous flow generally dominating at low flow velocities and high liquid viscosity, and inertial flow becoming more prominent with increasing flow velocity and decreasing liquid viscosity (i.e. higher Reynolds numbers). As the viscous regime is known to dominate for liquid foam systems, this has seen the majority of theoretical and experimental study. 3,9,20,36 Therefore, in order to maximise the applicability of the current study, this work will focus only on the viscous component of the liquid flow, D.
The current benchmark relationship between the viscous flow parameter, D, and Boussinesq number, B 0 , is the numerical solution of Nguyen, 28 which describes the data given by solving the Navier-Stokes equation for flow in an ideal 'infinite' PB (eqn (4)). This has been shown to agree well with subsequent model extensions and experimental analysis in macroscopic foam systems by Koehler, et al.
In order to relate eqn (4) to measureable variables, further considerations of PB geometry are required. Elias, et al. 9 describe the ideal relaxation profile of a vertical PB with thin films according to: which describes the transition of the PB profile from some initial radius, R 0 , to an equilibrium radius, R e , along the vertical axis, Z, according to a relaxation length, L. The equilibrium radius is analogous to the constant radius of the theoretical 'infinite' PB, which allowed further links to be drawn. Eqn (6) gives the theoretical relationship between R e and D according to measurable parameters: 9 R e = (DmQ/crg) 1/4 (6) where r is the liquid density, Q is the liquid flow rate and g is the acceleration due to gravity. As such, combining eqn (4) and (6) into eqn (7) yields a relationship between PB geometry, bulk liquid variables and surface mobility of the equilibrium (or 'infinite') PB with a thin film approximation.
In order to maximise the comparability of the results presented here to existing work, values obtained for B 0 are converted into surface shear viscosities, m s , using eqn (1) (with R = R e ).

The PB-node setup
The current setup formed an ideal PB and node by submerging and withdrawing a closed 3-legged frame geometry into surfactant solution (Fig. 1C). Flow rates, Q, were accurately controlled via a Cole-Parmer Dual Syringe Pump, which could directly inject surfactant solution into the upper PB through brass 3D print nozzles with outlet diameters of 0.4 mm. Larger diameter nozzles were found to cause pinning of the PBs to one side of the nozzle at low flow rates, thereby distorting the ideal profiles.
The frame was designed using Tinkercad software (Autodesk inc., USA) and 3D printed using stereolithography on a FormLabs Form s 2 printer, as this approach offered fast and simple production of precision tailored geometries. The distance between the legs and central tripod axis was set at x = 6.40 mm and heights were varied to produce PB lengths of 27.5 mm, 15.0 mm and 8.0 mm with absolute errors of approximately 0.5 mm. These heights were both above and within the ranges used in existing setups. 9,[29][30][31] The error in PB length resulted from the uncertainty in the exact end of the PB and beginning of the node.
The frame was mounted and enclosed within an acrylic and glass fronted box to allow imaging, minimise liquid evaporation and to remove air current disturbances. Access from a rear panel allowed investigators to lift and withdraw containers of surfactant solution from the frame. The whole system was then mounted on an optical rail to aid in PB profile imaging. PB profiles were obtained by illumination from the rear of the setup with a diffuse LED panel light and imaging via CCD camera from the front. A micrometre translation stage with translation in both the Z and Y directions, allowed precision movement of the CCD camera relative to the PB and node (Fig. 3a).
The resolution of images varied based on the choice of lens. Image resolution was calculated based on calibration images of a syringe needle with measured diameter 1.250 AE 0.001 mm. The maximum resolution of 1.5 mm per image pixel was used for curve fitting. Due to the high magnification at this resolution, it was necessary to combine multiple images along the Z-axis in order to visualise the full length of the PB profile ( Fig. 3b). Images were obtained in triplicate at 1000 AE 5 mm increments along the Z-axis, then processed and combined in Matlab (MathWorks, USA).
PB profile widths, d, were measured on a pixel row-by-row basis, using a Canny edge detection algorithm to determine PB boundaries. This method has been shown to be robust in detecting the presence of strong edges, while removing false edges produced by image noise. 39 The sensitivity of the edge detection could be varied by adjusting the upper and lower intensity thresholds used by the algorithm to determine edge strength. As such, it was possible to identify only PB profile edges, neglecting any fringe patterns in the adjoining films due to localised thinning.
Profile widths, d, were converted into PB radii, R, using the geometric correction factor 9 2 ffiffi ffi 3 p such that: Values for the PB Radii were averaged for each of the 1 mm incremental images along the Z-axis and subsequently combined into a single matrix of Z values and corresponding R values in the range 0 r Z r l 1 .
As was observed by Elias, et al., 9 around Z = 0 mm the PB exhibited a significant distortion due to its attachment to the liquid injection nozzle. This distortion was thought to mark a transition from a circular PB cross-section to an ideal crosssection ( Fig. 2) just below the point of injection. A vertical offset, Z offset , was therefore defined at the inflexion point of the distortion (Supplementary S1, ESI †), below which the ideal PB cross-section could be assumed and curve fitting undertaken.
The PB equilibrium radius was determined by fitting eqn (5) to values of R and Z using a non-linear least square fit method. R e values were then combined with measurements of m and r to calculate the left hand side of eqn (7). A least squares method was used to solve eqn (7) for B 0 , with these values subsequently converted to m s using eqn (1) with R = R e .
A solution of 0.50 wt% SDS was selected for analysis owing to its extensive use in previous work and its ability to easily produce stable films. This concentration was approximately double the Critical Micelle Concentration (CMC) of B0.235 wt%, 6 ensuring that observations would be of a fully populated interface, where little variation in surface tension could be assumed.
SDS solution was loaded into the PB-node setup and PB profiles measured for flow rates in the range 10 ml min À1 r Q Z 200 ml min À1 , therefore extending beyond the range of conventional foam drainage experiments. 29 This was in order to maximise the available data with which to investigate flow dependent trends.

Preparation of SDS solutions
SDS (499.9%) from Fisher Scientific (Loughborough, UK) was weighed using a digital balance to an accuracy of three decimal places. Purified water (15.0 MO cm) was weighed into borosilicate glass beakers and magnetically stirred at room temperature while SDS was added. Beakers were covered with cling-film to reduce evaporative water loss or contamination. Stirring continued for a minimum of 30 minutes prior to use of the solution. Great care was taken to ensure that glassware and stirring bars were thoroughly cleaned and rinsed with purified water, then air dried prior to use. Nitrile gloves were worn throughout handling to avoid sample contamination.
Bulk properties r and m were measured in triplicate using a Krüss Processor s Tensiometer K100 (Krüss GmbH, Germany) with density hook attachment and silicon density standard, and a Malvern Kinexus s Pro rheometer (Malvern Panalytical, UK) with Double-Gap geometry respectively. The averaged results with accompanying measurement errors were then used in the calculation of m s .

SDS solution properties
Values for specific bulk properties of the 0.50 wt% SDS solution, that are subsequently used for the calculation of m s (according to eqn (1) and (7)), are presented in Table 1.

PB relaxation and expansion
Images of PB profiles revealed a previously unreported phenomenon of PB relaxation followed by rapid localised expansion at discrete distances, Z swell , from the liquid injection nozzle at Z = 0 mm (Fig. 4a). Following the initial rapid expansion, a more gradual expansion of the PB followed, transitioning smoothly into the node. Measurements of the PB profiles clearly highlighted the Q dependence of Z swell , which determined the extent to which the PB profile was dominated by relaxation or expansion. The Expansion Fraction, F e , of the total PB length was described by F e = 1 À Z swell /l 1 , and can be seen to vary accordingly with Q in Fig. 4b. At lower flow rates, expansion was seen to almost completely dominate the PB profile until a critical flow rate, Q crit , after which F e decreased. Q crit increased with decreasing l 1 as shown in the inset of Fig. 4b, resulting in the increasing dominance of expansion for shorter PBs. Values of Q crit were 50 ml min À1 , 100 ml min À1 and 160 ml min À1 for l 1 E 27.5 mm, 15.0 mm and 8.0 mm respectively. These findings are highly significant. As flow rates and PB lengths approached those more typically observed in foam drainage experiments, 29 PB relaxation all but disappeared in favour of an expansion limited by the spatial geometry of the adjoining node. This therefore represents an important consideration to be made when choosing suitable theory to describe PB profiles.
The differing nature of all existing experimental work in this area meant that this phenomenon was not observed. 18,[29][30][31] Firstly, the few isolated PB and PB-node experiments mainly focused on shorter PB lengths, l 1 o 15 mm and the range of Q was smaller, 0 ml min À1 r Q r 100 ml min À1 in the case of Pitois et al. (2005aPitois et al. ( , 2005b. [29][30][31] In the only case where PB length was longer (l 1 E 40 mm), flow rates were in significantly higher increments of 5 ml min À1 in the range 0 ml min À1 r Q r 40 ml min À1 , with considerably poorer image resolution. 9 This indicates that previous experiments would either have observed only expansion-dominated PB profiles (e.g. Pitois et al. [29][30][31] ), or relaxation-dominated PB profiles (e.g. Elias, et al. 9 ), where the swelling region may have been obscured by a combination of distortion at the node or the point of liquid injection and insufficient image resolution. The range of flow rates and PB lengths studied by the current experimental setup therefore make it unique in its ability to observe both expansion and relaxation states simultaneously.
As the current theory does not describe the swelling phenomenon observed here, both this and the PB expansion regions should not be considered in the subsequent analysis of PB profiles. Instead, curve fitting of eqn (5) to experimental data is limited to the relaxation region only (see Supplementary S2, ESI †).

Length and flow dependence of PB geometry
Due to the fundamental differences noted between the present experimental setup and others, only limited comparisons could be made to existing experimental data. Neither the influence of a node on the upper PB nor the specific impacts of height variations are well documented in the literature; however, the differences in data obtained here serve to highlight the potential significance of these factors.

View Article Online
Experimental data were fitted to eqn (5) in order to extract values for the equilibrium radius, R e , which describes the theoretically ideal infinite PB, therefore validating the use of the theory described in Section 1.2. The fit quality was high in all cases, with a minimum r-squared value of 0.87 and the majority over 0.95. The values of R e were converted to equilibrium cross-sectional areas, S e , using eqn (3) with R = R e . These were compared to the minimum cross-sectional areas, S min , measured by Pitois, et al. 29 for a 3 g l À1 SDS solution with l 1 o 15 mm (Fig. 5). It is important to note here that values of S min may not have corresponded to equilibrium cross-sectional areas, which were the basis for the theory used in their subsequent calculations. For PBs of insufficient length, it is likely that values for S min would have been higher than S e , with a potentially significant impact on the observed flow dependent trends and magnitudes of calculated values.
As it was predicted above that Pitois, et al. 29 may only have observed an expansion-dominated PB, measurements were also made of the minimum cross-sectional area in the PB expansion region, S Exp , for Q 4 Q crit (Fig. 5). This would serve to highlight any changes in PB geometry that were characteristic of expansion rather than relaxation, where S Exp would better represent Pitois' S min than S e . The measurement of S Exp was taken from the PB profiles at the approximate point where the swelling region had transitioned into a continuous expansion towards the node (see Supplementary S2, ESI †). The difficulty in precisely defining these points is reflected in the error margins.
From Fig. 5 it can be seen that PB length had a significant impact on the magnitude of S e before Q crit was reached. In this region, S e was invariant to Q and increased in magnitude with decreasing l 1 . At flow rates above Q crit , values of S e from the different PB lengths began to converge. This suggested that the impact of PB length on the equilibrium PB cross-section could become insignificant for flow rates at the top end of those measured here and above.
The projected trend of increasing S e with decreasing l 1 could have conceivably produced magnitudes of S e comparable to the S min values reported by Pitois, et al., 29 with a PB length of approximately 4 mm. This would fall within their range of l 1 o 15 mm, however it does not account for the positive correlation between S min and Q. It was only from the values of S Exp for Q 4 Q crit that a significant increase in area with flow rate was observed. This suggested that the increase in S min with flow rate observed by Pitois, et al. 29 was indeed a characteristic of the PB expansion region, albeit without the delay imposed by Q crit . It is hypothesised that the PB-only nature of Pitois' setup may have been responsible for the lack of a Q crit in their data. By terminating in the bulk liquid instead of the strict geometry of a node, the expansion profile of the PB would have been far less restricted, potentially explaining the immediate increase in S min from Q = 0 ml min À1 .
Based on the observations made here, it seems unlikely that S min accurately represented an equilibrium PB cross-section in Pitois et al.'s measurements. This would have introduced a degree of error in values subsequently calculated using idealised PB theory, the consequences of which are discussed further in Section 3.4. It is emphasised here that the proposed method for extracting the true values of S e from experimental PB profiles is vital for the accuracy of subsequent calculations using ideal infinite PB theory.

Length and flow dependence of apparent PB surface viscosity
The flow rate dependence of the apparent surface shear viscosity, m s , was calculated as described in Section 1.2, the results of which can be found in Fig. 6.
It was seen that that decreasing l 1 produced values of m s that approached those of Pitois, et al. 29 Their average m s in the range 10 ml min À1 r Q r 100 ml min À1 was 1.57 Â 10 À5 g s À1 for l 1 o 15 mm, in comparison to a value of 2.16 Â 10 À5 g s À1 for l 1 E 8.0 mm determined here. It should be stressed that values published by Pitois, et al. 29 were for a slightly lower concentration (3 g l À1 ) of TTAB solution rather than SDS, however the authors reported their SDS data to have been similar. Eqn (6) and (7) highlight the importance of accurate measurements of the equilibrium PB cross-section when using idealised 'infinite' PB theory to calculate m s . Here it can be seen that overestimations of S e would have resulted in overestimations of m s . For the case of S min 4 S e therefore, as described in Section 3.3, it can be seen how this may have obscured flow dependent behaviour and increased the magnitude of m s measured by Pitois, et al. 29 The range of values shown here for m s with varying l 1 , showed a good agreement with those measured in the high precision setup of Zell, et al. 46 who reported m s t 10 À5 g s À1 for the true surface shear viscosity, in comparison to values of 10 À8 o m s o 10 À4 g s À1 for the apparent surface viscosity presented here. Zell's results are argued to be the most accurate measurements of SDS surface shear viscosity to date, where 10 À5 g s À1 represented virtually inviscid behaviour at the limit of their technique's sensitivity. The ability of the new technique to extract values as low as 10 À8 g s À1 therefore represents a significant step forward in measurement sensitivity. Furthermore, the ability to measure the flow dependent variations in m s at this scale could provide a wealth of information regarding the dynamics of the interface.
In order to interpret the origin of the flow and length dependent trends in m s , it was most useful to consider the shear imparted by the bulk liquid on the interface, which takes into account the variations in the PB cross-sectional area. The average liquid shear rate within the equilibrium PB, _ g s , for a given geometry and flow rate was approximated from eqn (9) (Elias,et al. 9 ), with the results shown in Fig. 7a. (9) Due to the flow rate independence of S e prior to Q crit , the resulting shear rates in these regions approximately scaled with Q. Here the apparent surface viscosity increased with decreasing _ g s , initially being well described by power laws, but beginning to deviate towards finite values of m s for the lower values of _ g s . For flow rates larger than Q crit a distinctive change in the behaviour of m s with _ g s was observed, reflecting the transition to a Q dependent S e (Fig. 7a). From this point, the data for all PB lengths began to converge abruptly, following what appeared to be a single trend of decreasing m s with _ g s . While additional PB length data is required to confirm this, it would imply that this flow region marked a response of the PB interface that was similar for all PB lengths. Indeed, the appearance of maximum shear rates for l 1 E 27.5 mm and l 1 E 15.0 mm, where m s appeared to increase, implied a resistance of the interface to  further apparent shear thinning. These maximum shear rates also corresponded to the beginning of expansions of the PB cross-sections (see Fig. 5) implying a causal link between interfacial resistance and expansion. The apparent shear thinning of the PB interfaces could have resulted from two primary effects: the compression of the interface, and/or the non-trivial combination of surface dilational viscosity and Marangoni forces with the surface shear viscosity, 8,15 both of which are influenced by rates of adsorption/desorption. A full numerical analysis of these effects represents a complex task that is beyond the scope of this study; however, initial comparisons to theory and experimental observations by other authors were able to provide a physical basis for the observed trends.
In a recent study by Elfring et al. 8 it was shown that very low surface viscosity interfaces, such as those described here, are expected to experience increasing Marangoni forces, F 1 , upon applied forces to the interface according to: where a = O(1) and represents the ratio of surface dilational to shear viscosities, and B 0 { 1. The divergence of this term as B 0 -0 was justified as a result of the limits of their problem, where instead it would be expected that F 1constant, as B 0 -0. The results of eqn (10) when applied to the values of B 0 obtained here are shown in Fig. 7b, where the scaling of F 1 with B 0 exhibited a very similar trend to that between m s and _ g s . This showed that the contribution of Marangoni forces to the apparent surface viscosity would have been expected to scale with increasing shear rates (see Supplementary S3, ESI †) and could therefore have been responsible for the apparent shear thinning behaviour. Furthermore, the apparent tendency of m s towards a constant value with decreasing _ g s noted previously, reflected the behaviour that would be expected from Fig. 7b were it not for the limitations of eqn (10) discussed above.
The theory of Elfring et al. 8 relies upon a limited degree of compressibility of the interface in order to form the concentration gradients necessary to produce Marangoni forces. The increase in surfactant concentration, G, at the interface results in a decrease in the surface tension that is well described by theory and experiment alike. 10,12,17,26,40,45 Fig. 7c shows the theoretical relationship using the Frumkin isotherm adsorption model with model parameters for pure SDS solution as described by Kinoshita et al. 17 At the CMC, the formation of micelles in the bulk liquid and the subsequent reduction in the energy barrier limiting surfactant desorption tends to reduce further compression of the interface. 5 This would imply that around this point, any differences in the interfacial structures caused by PB length dependence should have become less pronounced. The convergence of the m s data in Fig. 7a at flow rates above Q crit could therefore represent the tendency towards a common limit of interfacial compressibility for the different PB lengths.
Based on the analysis so far, a logical prediction for the strong PB length dependence of m s observed below Q crit was that it was likely to have been the result of differences in bulk and/or surface concentrations of SDS. Such differences would have affected the rates of adsorption/desorption of SDS as well as the compressibility of the interface at given liquid shear rates. Indeed, the effect of varying bulk concentrations on the compressibility of SDS interfaces has been directly observed in experimental work by Vogel. 41 In order to obtain an approximation of how SDS concentrations may have varied for different PB lengths, the average moles of surfactant contained in the liquid volume of each full-length PB profile was compared to the average PB interfacial area over the same length. This was done using eqn (11), producing a value for the maximum average available moles of SDS to populate an average unit area of interface, G max . Here the initial bulk concentration of SDS, C wt% , was 0.5 for all solutions, and the molar mass of SDS was 32 M = 288.38 g mol À1 . It is important to note that interfacial area and volume calculations were based on ideal PB crosssections with negligible film thickness. 9 The results of eqn (11) are plotted in Fig. 8, and show a marked difference between values of G max for the different PB lengths. This difference was approximately independent of the liquid shear rates, giving average values of 3.56 mol cm À2 , 2.90 mol cm À2 and 2.37 mol cm À2 , for l 1 E 8.0 mm, 15.0 mm and 27.5 mm respectively.
The results of Fig. 8 clearly show that an average variation in either the bulk or surface concentration of SDS must have occurred as a result of varying PB length. The reduction in concentration with increasing l 1 fits well with the observed trends in m s and the equilibrium cross sectional area. A lower initial surfactant population at the interface would have reduced the interfacial area due to the increase in the surface free energy, therefore increasing the liquid shear rates at given liquid flow rates according to eqn (9). Furthermore, one would expect a reduced m s for a lower surfactant population at the interface, due to the increased freedom of individual surfactant molecules. During compression experiments of SDS interfaces, Vogel 41 observed that the reduction in bulk concentration of SDS solutions resulted in a reduction in the compressibility of the interface. This would have served to reduce the variation in m s at lower SDS concentrations, in agreement with the results shown here.

Conclusions
In this study, we have proposed a novel technique for studying liquid flow through an isolated Plateau border of foam terminating in a geometrically 'ideal' node. The high resolution of the imaging achieved, combined with the flexibility of the setup to vary both liquid flow rate and the length of the PB, have proven capable of exploring a broad range of experimental conditions not seen in previous work.
An unexpected outcome of this study was the appearance of a clear distortion in the PB profiles that marked a sudden transition from the anticipated PB relaxation to a gradual PB expansion that transitioned smoothly into the node. These 'swelling regions' were found to be both PB length dependent and flow rate dependent after a critical flow rate, Q crit . After Q crit , swelling regions propagated discrete distances, Z swell , along the PB, thus determining whether the PB profiles were dominated by relaxation or expansion. While previous work was found to have observed profiles dominated by either expansion or relaxation, this marks the first instance where this dominance could be precisely controlled. As conditions approached those of conventional foam drainage (i.e. PB lengths less than 10 mm and flow rates less than 100 ml min À1 ), PB profiles were increasingly dominated by expansion into the node. The variable of PB length is therefore deemed an important consideration when selecting relevant theory to describe the evolution of PB profiles. We propose a further study of this region, wherein theory describing both PB relaxation and PB expansion is applied to the profiles studied here. A successful theoretical description of both regions should be able to describe the geometric evolution of their profiles, while predicting the appearance of the swelling regions at varying PB lengths and liquid flow rates.
Due to the nature of the theory presented here, further analysis of PB profiles was restricted to the relaxation regions.
Theory describing ideal PB relaxation agreed well with measured PB profiles, producing equilibrium cross-sectional areas, S e , that could then be accurately interpreted using idealised 'infinite' PB theory. The values of S e obtained here were shown to be consistently lower than the minimum crosssectional areas for similar systems in previous studies, 29,30 suggesting that direct measurement of minimum PB crosssectional areas provides a poor indication of the equilibrium PB geometry. In turn, it was shown that inaccurate measurements of the equilibrium PB geometry would have resulted in artificially high measurements of the surface shear viscosity when applying idealised PB theory. It is also thought that the PB-only nature of previous setups may have influenced measurements, requiring a more direct investigation of the limits imposed on the PB by the node geometry.
The magnitudes of the apparent surface shear viscosity, m s , and its dimensionless equivalent, the apparent Boussinesq number, B 0 , agreed well with values from the most precise existing measurements of true surface shear viscosities of SDS interfaces in the literature, 46 exceeding the measurement sensitivity by approximately three orders of magnitude. Values of m s were in the range 10 À8 g s À1 o m s o 10 À4 g s À1 corresponding to PB lengths within the range 27.5 mm \ l 1 \ 8.0 mm respectively, while corresponding values of B 0 were in the range 10 À4 o B 0 o 10 À3 .
The high sensitivity of the measurements here were able to detect a PB length dependence and flow rate dependence of m s has not been measured in previous studies. While a full numerical analysis of these effects was beyond the scope of the present study, a partial numerical and qualitative analysis is provided that indicates likely causes for the observed dependencies.
The apparent shear thinning behaviour of the PB interface is likely to have been the effect of a limited compressibility of SDS around the CMC, resulting in the significant increase in the contribution of Marangoni forces predicted numerically 8 for virtually inviscid interfaces such as those described here. Marangoni forces were predicted to have increased with the liquid shear rate that was driving the compression of the interface.
PB length dependence of m s has been attributed to differences in the bulk and interfacial concentrations of SDS that occurred at the different PB lengths. The average available moles of SDS in the bulk liquid per unit area of interface for the different PB lengths was measured to have decreased with increasing l 1 . This decrease was roughly independent of liquid shear rate, scaling with factors of 1.00, 0.82 and 0.67 for l 1 = 8.0 mm, 15.0 mm and 27.5 mm respectively. The variations in the bulk and surface concentrations of SDS are hypothesised to have resulted in changes to the compressibility of the PB interface, as observed in previous SDS compression experiments. 41 These observations indicated a lesser degree of interfacial compressibility for lower surfactant concentrations, which at the same time would be expected to lower initial values of m s , as was seen from the length dependence of m s here. The convergence of the length dependent m s data with increasing shear rates above Q crit , suggested a limit in the compressibility of the interfaces, wherein a common interfacial structure was approached for all PB lengths.
The complex interactions of the forces and geometries described here represent a significant challenge for the application of a full numerical analysis. Indeed, further work is still required in order to untangle these effects at the level of simple planar interfaces. 8 In order to better understand the contributions of the factors discussed, we therefore propose further studies of soluble surfactant systems with known variations in chemical and physical characteristics. Variations in bulk concentrations and viscosities, for example, would be expected to elicit changes in the measured m s that could then be used to better understand the effects proposed above.
Overall, this study has shown that the present technique represents a potential milestone in measurement sensitivity for the interfacial responses of soluble surfactants to applied shear. Furthermore, it directly relates these values to flows through ideal foam channels, a key aspect of macroscopic foam drainage models where the significance of m s is often contested for soluble surfactants. 46 It is hoped that the further exploration of soluble surfactants in this manner will yield insights into the dynamic nature of these foam interfaces, ultimately improving the current understanding of these complex systems.

Conflicts of interest
There are no conflicts to declare.