in fluid mechanics , the non - slip boundary condition ( nsbc ) states that fluid at the interface between a fluid and a solid surface is stationary with respect to that surface .
surface slip arises when there is non - zero relative motion between the interfacial fluid and the solid surface , in which case the nsbc does not apply . in the vast majority of experiments described historically
, the nsbc has been an adequate boundary condition @xcite .
however , finite surface slip has been the subject of much recent interest and promises to play a significant role in the current and future development of micro- and nanofluidic devices @xcite .
navier @xcite first described surface slip as early as the 1820s and his mathematical definition of the slip length has been widely used , for example in recent reviews @xcite and analyses of slip - dependent flow @xcite .
slip is a well - studied phenomenon in non - newtonian fluids such as polymer solutions @xcite . for slip of newtonian fluids ( newtonian slip ) ,
theory is well developed @xcite , but interpretation of experiments has been much more controversial .
newtonian slip is important because technologies can usefully be applied to newtonian fluids , but also because understanding of newtonian slip will provide a more profound understanding of fundamental slip mechanics .
the physical mechanisms governing slip are affected by variables that include surface chemistry , shear rate , surface structure , newtonian viscosity , non - newtonian behaviour at the molecular scale and combinations of each of these variables @xcite . nanoengineered superhydrophobic surfaces trap gas between the fluid and the surface , producing large slip lengths of the order @xmath410 @xmath5 m @xcite .
experimentally , newtonian slip has been investigated using a variety of methods @xcite .
the most widely used and precise technique involves measurement of the drainage forces between closely - spaced surfaces in an atomic force microscope ( afm ) or surface force apparatus ( sfa ) and applying an analysis based on reynolds theory of lubrication @xcite .
robust , non - zero values of @xmath0 for newtonian fluids at relatively smooth surfaces are typically of the order 10 - 100  nm @xcite .
experiments have not established firm quantitative correlations between slip and major experimental variables , such as contact angle @xcite . recently
, it has emerged with some clarity that in afm or sfa experiments , the nsbc can be expected to hold on smooth hydrophilic surfaces @xcite , while finite slip lengths of order 10 - 20  nm have been measured on some hydrophobic surfaces @xcite .
additionally , the long - held idea that surface roughness should inhibit slip @xcite is now supported by reasonable experimental evidence @xcite .
shear - mode oscillations can be used to drive fluid flow , and such systems pose scientific questions that are distinct from those relating to the study of steady fluid flows .
oscillating piezoelectric components are likely to prove very useful in future applications .
ultrasonic oscillators , well established as important tools for investigating thin films and bulk fluids , can also be used to probe the interfacial forces that are relevant in the study of slip .
kanazawa and gordon @xcite analytically coupled a half - space of newtonian fluid to a piezoelectric quartz crystal microbalance ( qcm ) , in which a solid interface oscillates in shear mode , using the nsbc .
further reports have presented experiments @xcite and analysis @xcite describing slip of a newtonian fluid adjacent to an oscillating surface .
a relatively large set of slip lengths measured using a qcm was obtained by ellis and hayward @xcite .
torsional quartz crystal oscillators operating at khz frequencies were initially developed and used for viscoelasticity measurements @xcite and have been used to measure the complex and viscoelastic behaviour of polymer solutions . more recent work has investigated dispersions of colloidal @xcite and micellar @xcite spheres .
the shear rates accessible using torsional oscillators are considered advantageously high for such rheological measurements @xcite .
a particular configuration , the composite torsional ultrasonic oscillator , was developed and analysed by robinson et al .
@xcite and robinson and smedley @xcite , who also derived an equivalent circuit .
recently , the first experimental results investigating newtonian slip at a smooth surface using such a device have been reported @xcite , for oscillations at 40  khz
.    there has been some confusion relating to fundamental aspects of slip , perhaps best characterised by neto et al .
@xcite , who noted that dialogue between the qcm community and other researchers interested in slip has been scarce . in order to address this issue ,
navier s definition of the slip length is used consistently in the current work .
an intuitive , consistent and general methodology for slip measurement at an oscillating surface is presented .
analytic inconsistencies in the literature are identified and clarified , especially those concerning an oscillating interface .
the composite torsional oscillator is introduced along with a full analysis for measurement of newtonian slip using this device .
preliminary experiments are presented to demonstrate how this device is used to probe interfacial forces in oscillatory flow .
we discuss the interpretation of measurements and the analytical assumptions on which they are based .
finally , methods of slip measurement are compared and the advantages of using the composite torsional oscillator are explained .
by definition , viscosity @xmath6 is constant and independent of shear rate throughout a newtonian fluid .
the velocity field @xmath7 of an incompressible newtonian fluid of density @xmath8 in the presence of pressure gradient @xmath9 and body force @xmath10 can be determined by solving the navier - stokes equation ,    @xmath11    microfluidic flows have characteristically low reynolds numbers , so the convective derivative @xmath12 is insignificant when compared with the viscous term @xmath13 .
solutions to the navier - stokes equation are dependent on flow geometry and boundary conditions .    the relevant solution for the current paper concerns a smooth surface oscillating in shear mode in contact with a fluid half - space . in the absence of significant body forces or pressure differentials , the equation of the shear wave induced in the fluid
is    @xmath14    where @xmath15 is the fluid velocity in the direction of the surface oscillations , @xmath16 is the cartesian direction perpendicular to the solid surface and @xmath17 .
the value of @xmath15 at the interface is @xmath18 , which depends on the fluid - surface boundary condition .
the time dependence of @xmath15 and @xmath18 is @xmath19 , where the frequency of oscillation is @xmath20 and @xmath21 is time .
this solution is relevant for ultrasonic oscillators .
boundary slip occurs when the nsbc is contravened , so that there is finite relative motion between a solid surface and an adjacent fluid .
navier @xcite characterised the slip boundary condition by defining the slip length @xmath0 . for flow adjacent to a planar , impermeable solid surface ,    @xmath22    here
@xmath23 is the shear velocity relative to and parallel with the surface and @xmath24 is the value of @xmath23 at @xmath25 .
navier s definition of slip has been widely used over a long period and has a strong physical basis relating to viscous force at the interface . the fluid shear rate at the interface that appears in eq .
[ eq : gen slip length ] also determines the viscous force ,    @xmath26    slip therefore determines the shear stress at the interface , which is analogous to friction in character .
the coefficient @xmath27 , defined as    @xmath28    is finite only when there is slip . using eqs .
[ eq : gen slip length ] , [ eq : gen viscous stress ] and [ eq : friction define ] , we find    @xmath29    it is noted that the nsbc and boundary slip are , respectively , analogous to the regimes of static and sliding solid - solid friction . extending the analogy , we might expect a wide variety of material and interfacial properties to determine the onset and magnitude of solid - liquid friction .
if slip length depends on shear rate , eq
.  [ eq : gen slip length ] becomes non - linear and rate - dependent formulations of slip ( e.g. @xcite ) become relevant .
currently , there is conflicting experimental evidence relating to possible shear rate dependence of @xmath0 @xcite .    a different physical interpretation of slip length relates @xmath0 to the distance into the solid over which the tangent to the flow profile at the surface must be extrapolated for the velocity to reach zero .
in this context , navier s slip length has correctly been described as a fictional distance @xcite and is labelled in fig .
[ fig : slip definition ] .
the definition of slip length used in some recent slip studies @xcite , and formalised by ellis and hayward @xcite as ( their notation ) @xmath30 , or    @xmath31    is distinct from navier s definition in eq .
[ eq : gen slip length ] .
ellis and hayward s definition corresponds to extrapolation of the flow profile a distance @xmath32 ( fig .  [
fig : slip definition ] ) relative to the surface when slip occurs .
the value of @xmath32 is equal to @xmath0 only in the case of a linear velocity flow profile .
an often - used schematic diagram @xcite , which uses a linear flow profile to depict slip length , appears to have been an important factor in the emergence of eq .
[ eq : b_1 ] .
we wish to point out that when the flow of a newtonian fluid is driven by an oscillating surface ( as introduced above ) , the shear rate introduces non - zero phase into the value of @xmath0 , which is defined using @xmath23 ( eq .  [ eq : gen slip length ] ) .
the existence of this phase highlights the limitations of direct quantitative comparison of measurements of @xmath0 using different experimental techniques .
when measuring slip at an oscillating surface , it is both intuitive and useful to use the dimensionless complex parameter @xmath33 .
this parameter , which simply quantifies the velocity discontinuity between the solid surface and the first liquid layer , simplifies the algebra describing the effect of damping on the equivalent circuit of a piezoelectric oscillator .
it is defined by interfacial properties and does not depend on the flow profile ,    @xmath34    the velocity of the solid surface is @xmath35 and the fluid oscillation is assumed to be dominated by a sinusoidal oscillation of angular frequency equal to that of the solid surface . under this assumption , @xmath36 is equivalent to parameters introduced by ferrante et al .
@xcite and used elsewhere @xcite for modelling of slip at an oscillating surface ; to our knowledge , it has not previously been applied to experimental measurements with full consideration of phase dependence .
more generally , @xmath36 can be decomposed into fourier components , an approach that has been used when modelling nanobubble - induced slip @xcite .
the nsbc is consistent with the values @xmath37 and @xmath38 .    for an oscillating surface
, flow is relative to surface motion , so    @xmath39    the magnitude and phase of @xmath40 relative to @xmath35 are plotted in fig .
[ fig : slip phasor ] , for both slip and non - slip cases . figure
[ fig : c ] shows that the relative amplitude of oscillation is uniformly non - zero at the surface when slip is present , but varies significantly near the interface ( @xmath41 ) depending on the values of @xmath42 and @xmath43 .
further from the surface ( @xmath44 ) , the relative amplitude tends towards a stable value close to 1 but dependent on @xmath42 .
the relative phase follows a similar trend in fig .
[ fig : d ] , showing strong dependence on the combination of @xmath42 and @xmath43 near the interface .
there is a phase discontinuity at @xmath45 for the cases in which @xmath46 .
the plane of discontinuous phase is consistent with a negative slip length ( fig .
[ fig : alpha vs b ] ) and could , for example , be interpreted as the edge of a stagnant boundary layer .
further from the surface ( @xmath44 ) , the relative phase converges towards the nsbc value , which tends towards zero a long way from the interface .    using this flow profile ,
slip length can be directly related to @xmath36 ( eqs .
[ eq : gen fluid wave ] , [ eq : gen slip length ] , [ eq : gen alpha define ] and [ eq : gen osc flow profile ] ) :    @xmath47    equation [ eq : gen alpha equals ] is plotted in fig .
[ fig : alpha vs b ] using the notation @xmath48 .
figure [ fig : alpha_b_amp ] shows that @xmath42 decreases when @xmath49 is within @xmath50 of @xmath51 , so that the slip length @xmath0 is positive .
the phase plot ( fig . [ fig : alpha_b_phase ] ) demonstrates that @xmath52 when the fluid oscillates in phase with the surface ( @xmath53 ) .
for the nsbc , @xmath54 , so @xmath49 can take any value .
the maximum size of the phase shift @xmath43 is limited by slip length magnitude .
the complicated nature of the flow profile near the interface supports the use of the parameter @xmath36 , which does not depend on the flow other than at the interface itself , in contrast to @xmath0 .
however , calculation of the shear rate using the flow profile is still necessary for derivation of @xmath36 from surface force experiments .
the damping force on an oscillating surface is equal and opposite to the viscous force in the adjacent newtonian fluid .
when slip occurs , the damping stress is ( eqs .
[ eq : gen fluid wave ] , [ eq : gen viscous stress ] and [ eq : gen osc flow profile ] ) ,    @xmath55    analysis of ultrasonic oscillators is greatly simplified by using the equivalent circuit analogy , in which mechanical components are replaced by their electronic equivalents .
we now incorporate the damping stress into the analysis of a mechanical circuit .
the stress at the interface ( eq .  [ eq : gen interface stress ] ) , which has components both in phase and out of phase with surface velocity , can be written as    @xmath56    the impedance of the system is @xmath57 , as previously applied to a torsional oscillator by bergenholtz et al .
@xcite , who used the nsbc .
this impedance is dependent on the liquid , the surface , and the interaction between them .
the real part of the impedance ( @xmath58 ) determines the stress opposing , and in phase with , the surface velocity ,    @xmath59    this impedance is viscous in nature and dissipative of energy . in an equivalent circuit ,
velocity is analogous to current and @xmath58 is like a resistance .
energy dissipation is determined by the quality factor of the oscillation .
the inertial term ( @xmath60 ) gives the stress that is in phase with the acceleration of the surface ,    @xmath61    an increase in the inertial term is equivalent to the addition of mass to the surface of the oscillator . in the equivalent circuit , the time derivative of
current is analogous to acceleration ( modified by a factor of @xmath20 ) , so @xmath60 is an inductance .
a change in @xmath60 is measured by considering the change in the period of oscillation . under the nsbc , @xmath58 and @xmath60
are equal for newtonian fluid damping .
slip - induced change in the relative sizes of these terms indicates a finite value of @xmath43 .
the theory that has been presented implicitly addresses two conceptual pitfalls that have caused confusion and misinterpretation in a range of previous work .
the first pitfall concerns the variant slip length @xmath32 , defined in eq .
[ eq : b_1 ] , which is only equal to @xmath0 in the limit of a linear flow profile .
this limit is generally not applicable for slip measurements .
the difference between @xmath0 and @xmath32 is greater than 10% for typical slip measurements using afm , sfa or ultrasonic oscillation methods .
we further note that ellis and hayward s definition has been used inconsistently ; both with and without considering a taylor expansion @xcite .
navier s slip length @xmath0 should be used at the present time due to its strong physical basis and to promote consistency across experimental studies .
the second conceptual pitfall relates to those experimental studies of slip at oscillating surfaces which have imposed a restriction on the phase of slip @xcite , especially the assumption that the slip length @xmath0 is real - valued .
it is equivalent to assume that @xmath27 or its reciprocal @xmath62 is real - valued , or more generally , that the damping force is in phase with the surface velocity ( eq .  [ eq : b friction ] ) .
several lines of reasoning lead us to emphasize that @xmath0 must be allowed to have non - zero phase .
we have shown that , because the slip length is defined using the shear rate ( eq .  [ eq : gen slip length ] ) , the @xmath16-dependence of the fluid velocity ( a complex exponential , eq .  [ eq : gen fluid wave ] ) introduces phase to the value of @xmath0 .
the usual ( nsbc ) form of solid - fluid friction force is viscous ( eq .  [ eq : gen interface stress ] ) and therefore @xmath63 out of phase with the velocity above an oscillating surface .
figure  [ fig : alpha vs b ] shows that the phase of @xmath0 in the limit approaching the nsbc is non - zero and that both magnitude and phase must be determined for either @xmath0 or @xmath36 to be fully described . restricting @xmath0 to real values
is equivalent to restricting @xmath36 to one degree of freedom in eq .
[ eq : gen alpha equals ] .
furthermore , experiments @xcite have indicated that a fluid can be expected to oscillate in phase with , or slightly lag , the surface oscillation .
neither interpretation is consistent with a real value of @xmath0 .
having established that @xmath0 must be considered complex , we address specific arguments made in previous studies in support of a real - valued @xmath0 .
firstly , there has been a direct assumption @xcite that the value of @xmath27 ( or @xmath62 ) should not be phase dependent .
this assumption has been based on rodahl and kasemo s argument @xcite , which links @xmath27 to the ratio of force to velocity for low reynolds number drag on small bodies ( e.g. stokes drag ) , in which case @xmath27 is a real constant .
this argument ignores the fact that important factors in the determination of stokes drag include the contribution of pressure surrounding the particle , the requirement for steady non - oscillatory flow and use of the nsbc .
we note that friction and drag parameters are usually real because they are not used to describe oscillatory flow . in conjunction with this point
, we note that the term coefficient of friction , which has frequently been used to describe @xmath27 in slip studies , is distinct from the usual parameter of that name , which relates the shear force to the normal force in the case of solid - solid sliding friction .
secondly , an argument for phase restriction @xcite has drawn on the correlation between bond strength and bond length , which should determine magnitude and phase of slip respectively .
such an argument does not explain why the phase should take one particular value , or why such a correlation between magnitude and phase should not be experimentally investigated .
non - zero phase is not significant in previous work relating to lower - frequency oscillations in afm or sfa devices @xcite , in which @xmath64100  @xmath5 m for 39  hz oscillations in water , as opposed to @xmath4500  nm for a 1  mhz oscillator .
similarly , periodically patterned microchannels @xcite are adequately described by a real - valued slip length when the ratio of the flow velocity to the pattern wavelength does not exceed @xmath41  khz , corresponding to patterns of extremely short spatial wavelength .
the previous discussion of slip is now applied to the composite torsional ultrasonic oscillator ( fig .
[ fig : to_schematic ] ) .
the damping torque when the specimen rod is submerged in a newtonian fluid was derived by robinson et al .
@xcite under the nsbc , using an equivalent circuit analysis . in that derivation ,
the shear damping stress on a flat surface is integrated over the surface of the specimen rod .
appendix  a , in this paper , gives the same derivation with slip included and appendix  b discusses an assumption that is used .
similar derivations have considered damping of other torsional oscillators by viscoelastic fluids under the nsbc , yielding similar functional forms @xcite . for the oscillator of specimen rod radius @xmath65 described by fig .
[ fig : to_schematic ] , the total torque for @xmath66 or @xmath67 immersed flat surfaces is    @xmath68    where the subscript max refers to the value at the antinode of the standing wave , which has wavelength @xmath69 .
when a specimen rod is immersed , the resistance and inductance added to the equivalent circuit are @xmath70 and @xmath71 respectively , so    @xmath72    equation [ eq : to t equivalent ] is directly analogous to eq .
[ eq : components of interface stress ] for this specific configuration if we substitute notations @xmath73 and @xmath74 , consistent with the resistive and inductive descriptions of these parameters .
comparing eqs .
[ eq : to nsbc t ] and [ eq : to t equivalent ] and introducing a constant @xmath75 , we find    @xmath76    the equivalent circuit reduces to a series lcr circuit @xcite .
@xmath71 is therefore related to the change in period upon immersion @xmath77 and the resonant period @xmath78 , using    @xmath79    the inverse quality factor of the oscillation @xmath80 is given by the total resistance of the composite oscillator @xmath81 and the total inductance @xmath82 :    @xmath83    @xmath70 is calculated by assuming that @xmath82 is dominated by the mechanical inductance ( moment of inertia ) of the rods @xcite , so that ( @xmath84 )    @xmath85    @xmath86 is the change in @xmath80 upon immersion .
when piezoelectric elements are incorporated , the equivalent circuit for the torsional oscillator gives the inverse quality factor ,    @xmath87    where @xmath88 and @xmath89 are voltages applied to the drive and gauge crystals respectively , and @xmath90 is a calibrated constant for a particular oscillator arrangement @xcite .
the quality factor is also given by    @xmath91    where @xmath78 is the resonant period and @xmath92 and @xmath93 are measurements of the period when @xmath88 is double the resonance value with @xmath89 held constant . by plotting the right hand side of eq .
[ eq : to q^-1c ] against the ratio @xmath88/@xmath89 , eq .
[ eq : to q^-1 ] can be checked and the value of @xmath90 determined .
@xmath90 is independent of specimen rod damping , but it is dependent on the rigid moment of inertia of the specimen rod , so changes if mass is added to the specimen rod .
we now consider comparison of a base case ( subscript 0 ) in which the nsbc holds ( @xmath94 , @xmath38 ) , with a case in which slip is present ( subscript 1 ) . in practice ,
it is difficult to determine whether there is slip in the base case , so @xmath36 can be used to measure the relative amplitude and phase of slip between two cases investigated .
@xmath75 remains constant between measurements as long as @xmath95 , @xmath69 , @xmath6 , @xmath8 and @xmath20 are constant , so using eqs .
[ eq : r_eta ] , [ eq : to nsbc tor ] and [ eq : to r_eta ] , we define @xmath96 and @xmath97 as    @xmath98    and    @xmath99    rearranging eqs .
[ eq : r_foralpha ] and [ eq : l_eta ] , @xmath36 is calculated directly from measurement :    @xmath100    and    @xmath101    equations  [ eq : alpha measure mag ] and [ eq : alpha measure phase ] are plotted in fig .
[ fig : alpha from measurement ] . figure  [ fig : alpha amp ] shows that @xmath42 increases with both @xmath96 and @xmath97 , indicating that if mass is decoupled from the surface , the amplitude of fluid oscillation will decrease unless more power is applied .
figure  [ fig : alpha phase ] shows that to retain the nsbc value of phase ( @xmath38 ) , @xmath96 must decrease if @xmath97 decreases . therefore , an oscillation generating relatively low power will tend to produce a value of @xmath102 ( fluid oscillation leading ) unless there is some decoupling of mass from the surface , and therefore reduced inertia .
similarly , if more power is generated , the fluid will lag unless balanced by greater inertia .
considered together , fig .
[ fig : alpha amp ] and fig .
[ fig : alpha phase ] show that a non - zero measurement of either @xmath96 or @xmath97 necessarily indicates a departure from @xmath37 , or @xmath38 or both .
the introduction of a slip length is mathematically equivalent to considering a viscoelastic fluid . by inspection of eq .
[ eq : gen interface stress ] , and remembering that @xmath103 , we find that instead of using @xmath36 , we can introduce a complex viscosity @xmath104 to replace the newtonian viscosity @xmath6 such that @xmath105 .
due to coupling to the decaying wave in the fluid , we have @xmath106 , where @xmath107 is the complex shear modulus .
we find that the storage modulus @xmath108 and the loss modulus @xmath109 are given by combinations of the inertial and dissipative terms in the oscillation @xcite ,    @xmath110    and    @xmath111    it is therefore misleading to use the viscoelastic description of storage and loss moduli for the current situation .
it is also misleading to describe a newtonian fluid as viscoelastic until there is experimental evidence for viscoelastic behaviour rather than some other cause of slip . some previous qcm work based on fully phase - dependent equivalent circuit analyses
@xcite has considered slip in conjunction with a viscoelastic fluid and layering , thereby confusing possible newtonian slip with other effects ; for further effects of layering see the discussion ( below ) .
such an approach may be motivated by layering of a qcm electrode substrate , that is not present when using a torsional oscillator .
one previous investigation @xcite has noted the equivalence of viscoelasticity and slip in interpretation of slip measurements .
to demonstrate investigation of slip using a torsional oscillator , preliminary experiments have been carried out using a 40  khz oscillator , as described by fig .
[ fig : to_schematic ] .
the fused quartz specimen rod was uniformly ground using al@xmath112o@xmath113 grits to remove surface features of magnitude @xmath1145  @xmath5 m , and polished using cerium oxide powder . polishing debris
was removed by sonication in acetone , isopropanol , 10% nitric acid and ethanol .
any remaining organic material was removed using piranha solution ( 70  vol% sulphuric acid , 30  vol% hydrogen peroxide ) .
the rod was then thoroughly washed in deionised water ( 0.2  @xmath5 m filtered , 18.2  m@xmath115  cm ) and lightly flame - polished to minimize surface roughness .
the slip case was configured by coating the specimen rod with a hydrophobic fluoropolymer surfactant ( rs components ) , evenly dispensed from an aerosol can .
the base case was retrieved by dissolving the surfactant in isopropanol prior to re - cleaning the rod using piranha solution and washing with deionised water and ethanol .
a pre - optically polished , flame - polished piece of fused quartz has rms roughness @xmath116  nm and peak to trough roughness @xmath117  nm , as measured over several widely - distributed areas of size 5  @xmath5m@xmath118 using afm .
previously @xcite , scanning electron microscopy ( sem ) revealed peak to trough roughness of 24  @xmath5 m on rods left unpolished following grinding using al@xmath112o@xmath113 , with no significant alteration of the texture following the addition of the surfactant .
however , a flame - polished surface coated with the surfactant has additional texture characterised by features of spacing @xmath119
@xmath5 m with rms roughness 24  nm and peak to trough roughness @xmath120  nm , observed using optical microscopy and afm .
the advanced static contact angle of deionised water on fused quartz is @xmath121 .
the advanced contact angle on the same surface coated with the surfactant is @xmath122 .
the texture of the coating causes some hysteresis , and the receded contact angle is @xmath123 .
contact angles were measured using the sessile drop method ( cam 200 , ksv instruments ) .    prior to each immersion
run , the specimen rod was washed in deionised water and dried in air . for the base case ,
the uncoated rod was also rinsed with ethanol prior to drying .
deionised water was used as the test fluid in all experiments .
the oscillator was contained in an insulated bell jar standing on a plate maintained at 303  @xmath124  1  k and situated on a large mechanically damped pad usually used for a transmission electron microscope .
measurements from recent work are presented alongside previously reported experiments @xcite in table  [ tab : results ] and explained in this section .
the data points for @xmath96 and @xmath97 are each the result of several immersion events in each of the base and slip cases .
the specimen rod was immersed by slowly raising a filled vessel to the rod .
a value of @xmath86 or @xmath77 for any particular immersion run is the difference between the dry value and the value when the water was in contact with the node of the standing wave , adjusted to account for any background drift .
the fluid meniscus was aligned with the vertical mid - point of a rubber o - ring that was fitted half way along the specimen rod .
the o - ring , which fixes the position of the node , was located at a point where surface motion is zero or minimised , and is considered rigid during calibration .
each immersion was carried out after the oscillator had been left for at least 90 minutes , so that instrumental drift per second was @xmath125 for @xmath88/@xmath89 and @xmath126  s for @xmath78 , as measured over 10 minute intervals .    the standard error of the mean measurement was combined with error ascribed to drift to give the total experimental error in @xmath86 and @xmath77 in each case .
uncertainty is dominated by variation between different immersion events , probably due to the extreme thermal and mechanical sensitivity of the instrument . the lesser contribution due to instrumental drift
is accounted for by assuming that the mean drift correction is a random error .
earlier experiments @xcite have established that uncertainty due to surface tension effects and location of the node is small relative to the standard error in the mean value , because very little damping occurs at or near the node .    in order to calculate @xmath86 from measured voltages ,
the half width of the oscillation was measured directly following each immersion .
equations  [ eq : to q^-1 ] and [ eq : to q^-1c ] were then used to calculate a precise value of @xmath90 .
it has been shown @xcite that a plot having gradient @xmath90 is typically linear ( @xmath127 for all fits ) and independent of temperature , the immersion level and the presence of surfactant .
the two former variables are held constant during immersions , but the slip case experiments were calibrated separately from the base case in the most recent experiments .
the value of @xmath90 is slightly lower for the slip case , as expected from the first principles analysis @xcite in which @xmath90 is inversely proportional to the oscillator s moment of inertia .
the maximum error in @xmath90 , conservatively estimated as the difference between the slip case and the base case , is not significant in comparison with the experimental error of @xmath36 . @xmath90
changes when the torsional oscillator is mechanically altered , but was constant for each individual shear rate presented in fig  [ fig : alpha_expt ] .    for the nsbc , it is expected that @xmath128/@xmath78 @xcite .
this relation was followed to within 5% for the majority of immersions .
even accounting for the possibility of slip , immersion events with @xmath11410% deviation clearly indicate a significant instability in the oscillation and such data were discarded .
experiments have been performed at shear rate amplitudes @xmath129 between 230 and 6840  s@xmath3 , corresponding to linear displacement amplitudes of 3.2 and 96  nm at the antinode .
the linear amplitude of the oscillation is equal to @xmath130/@xmath20 .
the shear rate at the surface of the rod is calculated as described by robinson et al .
@xcite . for this calculation ,
rod diameters and masses were measured using a micrometer screw gauge and a mettler ae 200 mass balance , respectively . for each strain rate , intrinsic gauge capacitance @xmath131 was measured to better than @xmath132 by finding the gradient of straight line plots ( @xmath133 ) of the fractional change in resonant period as a function of total gauge capacitance @xcite .
the transformer ratio linking the electronic and mechanical parts of the equivalent circuit was calculated using @xmath90 .
the alternative calculation ( from first principles ) incorporates large error relating to the electrode separation angle 2@xmath134 ( fig .
[ fig : to_schematic ] , @xcite ) .
the final uncertainty in shear rate amplitude is @xmath135 .
however , the uncertainty in each value relative to the others depends only on @xmath89 , which is precise to better than @xmath136 .
the results for the five experiments in table  [ tab : results ] are plotted as slip parameters in fig .
[ fig : alpha_expt ] . figures  [ fig : ex alpha amp ] and [ fig : ex alpha phase ] appear similar because , for all five data points , mean values of @xmath137 and @xmath138 have the same sign .
all data lie within experimental uncertainty of the nsbc .
for the three data points from ref .
@xcite , the consistency of mean values ( @xmath139 , @xmath140 ) previously led to the conclusion that finite slip had likely been observed . although the newer results are not directly comparable due to differences in surface preparation , the mean values are now scattered around @xmath37 , @xmath38 , suggesting that the nsbc is observed with a high degree of accuracy .
the most recent experiments have covered a greater range of strain rates and afforded greater precision due to improved stability of the oscillation and better thermal and mechanical isolation .
the corresponding mean values of @xmath141 are now less than 30  nm , and the upper limits bounded by experimental uncertainty are less than 60  nm . for each data point , the lower error bar of @xmath142 is at zero because the nsbc falls within experimental uncertainty . because of the phase discontinuity at the nsbc ( fig .  [
fig : alpha vs b ] ) , any value of @xmath49 is possible for near - nsbc combinations of the measured variables .    previously , qcm measurements @xcite have obtained values of @xmath0 of similar precision , although in that case finite slip lengths of magnitudes between 10 and 80  nm were observed on bare as well as a range of hydrophobic surfaces . as discussed previously ,
the analysis used to derive these values contains inconsistencies .
when considering other results in the literature , we emphasize that comparisons are drawn from data sets taken from measurements of very different flows ( see discussion ) and that the current data are observed over a well - defined , wide range of shear rates .
it is therefore remarkable that the data is consistent with emerging trends from high - precision afm studies .
namely , slip at a hydrophobic smooth surface occurs on scales @xmath143  nm , if at all ; and surface roughness such as that observed for the slip case in the current work may play a part in suppressing any slip at a hydrophobic surface .
in order to discuss the interpretation of experimental slip measurements , we first clarify the terminology used in the discussion .
lauga et al.s definitions @xcite of indirect measurement , apparent slip and effective slip are used , although we extend the definition of effective slip to include temporal as well as spatial averaging .
true slip is an oft - used term not defined by lauga et al .
we define it here as the case when the mean tangential velocity component of fluid is different from that of a solid surface immediately in contact with it .
experimental methods are predominantly indirect at present .
even the most direct optical methods of slip measurement have been labelled indirect @xcite , in the sense that the velocity of tracers rather than the fluid itself is measured . in the case of the oscillator , the measured damping force is not necessarily consistent with the idealised analysis presented in the theory section . in any surface force measurement ,
the flow profile is inferred rather than explicitly determined . under these conditions ,
slip parameters @xmath0 and @xmath36 only describe the slip dynamics in full when rate - independent , true slip occurs and the fluid is entirely newtonian .
caution is required when analytically linking postulated mechanisms to apparent or effective slip measurements .
for example , it often appears convenient to assume the presence of layering at or near the slip interface @xcite .
typically , layers of unknown thickness and interfacial properties only introduce more modelling parameters .
as discussed previously , introduced viscoelasticity is indistinguishable from slip . in two studies employing @xmath36 @xcite , multiple viscoelastic layers with multiple boundary conditions
were modelled and several parameters were fitted .
the resulting magnitudes of @xmath36 ( ranging from @xmath144 to @xmath145 ) were larger than reasonable for a newtonian fluid .
at present , mechanisms of apparent slip are poorly characterized , so empirical measurement of slip parameters remains the most important step .
newtonian fluids should be considered newtonian right up to the interface to avoid introducing error .
sfa and afm techniques are currently the most well developed and widely - employed methods of measuring surface slip , and can reach resolution of better than @xmath146  nm for @xmath0 @xcite . optical methods such as the fluorescent recovery after photobleaching ( frap ) @xcite and particle image velocimetry ( piv ) @xcite techniques allow resolution of @xmath124100  nm , limited by diffraction and tracer particle size . recently ,
a nano - piv system has reportedly resolved true slip lengths of less than 100  nm with uncertainty less than @xmath12420  nm @xcite .
several groups have measured slip using capillary techniques @xcite , but accuracy has been limited by capillary fabrication issues , the difficulty of controlling surface chemistry within a capillary and ( in some cases ) the effects of electro - osmotic flow . a recently - developed technique derives slip from measurement of thermal motion of confined colloidal particles @xcite .
a slip length of 18  @xmath124  5  nm has been reported at near - zero shear rate .
torsional oscillator experiments presented in the current work confirmed the nsbc to within @xmath147 for @xmath42 , @xmath148 for @xmath43 and to within 60  nm of @xmath2 .
we stress that there are no significant limitations on the precision of the oscillator other than control of atmospheric conditions , vibration and electrical noise .
further , the technique offers advantages and points of difference when compared with other techniques
.    a particular advantage of the torsional oscillator is characterization of experimental shear rate .
the shear rate amplitude is easily varied by adjusting the applied gauge voltage and can at least cover the range @xmath149 to @xmath150  s@xmath3 . a major reason for use of torsional oscillators in rheological experiments is the large and difficult - to - reproduce range of accessible shear rates @xcite .
in contrast , the variable velocity distribution across a qcm causes analytic complications @xcite . in the afm and sfa techniques @xcite ,
the shear rate varies with respect to surface separation and distance from an axis passing through the centres of the approaching curved surfaces .
the shear rate is maximized away from this axis and also depends on the approach or withdrawal rate of the two spheres @xcite .
the second major advantage of the torsional oscillator is versatility .
surface treatments can easily be applied to the specimen rod , either while remote from the oscillator or attached to it
. if the specimen rod is detached , the response of the oscillator is easily recalibrated .
electronic voltage and period measurements are captured near - instantaneously , so time - dependent trends can be measured .
the oscillator is perfectly suited to measurement of in - situ changes of surface chemistry or other parameters , and experimental precision would be significantly enhanced for such experiments .
the oscillator can readily be used at variable temperatures and even variable pressures .
although it is best suited to differential measurements ( similarly to optical techniques ) , the zero - slip response can be determined by careful calibration of all the parameters described in the original equivalent circuit analysis @xcite .
in contrast , use of qcms and other torsional oscillator designs @xcite for slip measurement presents difficulties with regard to surface modification and incorporating piezoelectric equations into an equivalent circuit analysis .
an important distinction between damping - oscillation methods and other techniques is that the shear damping force is measured .
sfa and afm methods measure drainage forces caused by pressure in the flow around curved surfaces ; viscous forces are not significant .
another difference is that , when using the torsional oscillator , slip occurs over a relatively large fluid - solid interface , so measurement of molecular - scale effects is limited by surface homogeneity .
however , measurement of effective slip will give more accuracy and precision than a value derived from a small system .
it is not true that measurements of effective slip must take place on a scale @xmath151 in order to have a measurable impact @xcite , because any fundamental slip mechanism should be consistent over the whole interface .
regardless , we note that the linear displacement amplitude for the torsional oscillator is of comparable magnitude to the slip length , so molecular - scale interactions are certainly probed using this device .
in this paper , we have drawn particular attention to the occurrence of slip at an oscillating interface .
the ultrasonic composite torsional oscillator is a versatile instrument for measurement of viscous damping forces , and therefore slip . in order to fully characterize this device for slip measurements
, we have revisited fundamental aspects of the slip analysis .
the analysis has then been applied at a generalised oscillating surface , followed by the surface of the oscillator s specimen rod in particular .
preliminary experiments at shear rates between 230 and 6800  s@xmath3 have probed slip at a hydrophobic surface of rms roughness 24  nm , obtaining resolution better than @xmath147 for @xmath42 , @xmath148 for @xmath43 and within 60  nm of @xmath2 .
navier s slip parameter , while retaining the units of length , is best suited to interpretation as a surface force effect .
measurement of surface forces has played an important role in developing the understanding of slip to date , while applications are likely to utilise the lubricating effect at a slipping surface .
the actual velocity discontinuity at an oscillating interface , as encapsulated by the parameter @xmath36 for oscillating surfaces , is being directly observed with increasing precision using optical methods .
the link between a velocity discontinuity and a surface force necessarily utilizes the spatial gradient of an analytical flow profile , unless both types of experiment can be performed on the same system .
ongoing experiments will aim to significantly develop the stock of published slip measurements at oscillating surfaces .
particular experimental topics will be the influence of shear rate , hydrophobicity and surface texture , while the development and in - situ characterisation of surfaces or fluids with switchable properties will also be a priority .
thanks are due to dr .
shen chong for assistance with setting up the apparatus described here and used previously @xcite .
in this section , the damping torque on a composite torsional ultrasonic oscillator rod , as illustrated in fig . [
fig : to_schematic ] , is calculated .
we follow robinson and smedley @xcite and incorporate the case in which there is finite slip , as introduced in the theory section . using the cylindrical polar co - ordinate system defined in fig .  [
fig : to_schematic ] , the velocity in the azimuthal ( @xmath152 ) direction is given by    @xmath153    where the subscript max refers to the value at @xmath154 and at the antinode of the standing wave , which has wavelength @xmath69 .
there is a velocity gradient in the fluid in both the @xmath155- and @xmath16-directions .
viscous forces in the fluid exert azimuthal shear stresses on the rod , which dampen the torsional oscillation .
forces due to the @xmath16-gradient are neglected ( appendix  b ) .
the total damping torque is derived by integrating the tangential stress ( eq .
[ eq : gen interface stress ] ) along the immersed , curved length of the rod @xmath156 .
it is assumed that the fluid surface coincides with a node and the free end is an antinode , so that @xmath157 .
the torque is    @xmath158    the torque on the oscillating flat end of a rod is derived similarly .
the surface velocity at radius @xmath155 from the cylindrical axis is given by    @xmath159    and the torque is then integrated over the surface of the free end , so that    @xmath160    neglecting edge effects between the flat end and the curved surfaces , the total torque for @xmath161 immersed flat surfaces is    @xmath162
the velocity flow field around the torsional rod has gradients in both the @xmath155- ( radial ) and @xmath16- ( axial ) directions .
the azimuthal stress due to these gradients are given by                                                                                and @xmath32 defined in the text , with the shaded area representing solid material .
navier s slip length @xmath0 ( eq .  [ eq : gen slip length ] ) is only equal to the variable @xmath32 ( eq .  [ eq : b_1 ] ) when the velocity profile is linear .
the labels @xmath166 and @xmath16 denote cartesian coordinate axes.,width=321 ]