We continue earlier studies on the indicated configuration, improving previous estimates, providing explicit expressions for the relevant forces and a formal algorithmic procedure for their calculation, and sharpening and extending the predictions for qualitative distinctions among varying types of behavior that can occur. We include graphical representations for some of the more significant relations, as an aid to interpretation and for eventual design of experiments to test the physical relevance of the new material.
1 Background remarks
The present work is a continuation of  and of , where the behavior of the solutions of the capillary equation for the surface height of liquid in an infinite tank is described, in terms of the contact angles of the liquid with two infinite parallel plates that are partially immersed into the liquid and held rigidly. We describe here an algorithmic procedure for explicitly calculating the forces of attraction or repulsion between the plates and , depending on the respective contact angles and with the fluid, on the sides of the plates that face each other. As pointed out in  and in , the net forces in question do not depend on the contact angles at the triple interfaces on the opposite (outer) sides of the plates; that is a consequence of the hypothesis that the fluid surface extends to infinity in the two directions exterior to the plate configuration and orthogonal to it. We may thus concentrate attention on the integral curves for the fluid height on a section of the channel joining the plates; these curves are determined as solutions of the ‘capillary equation’
where ψ is inclination of the curve with the horizontal, and u is the height above the level at infinity. This equation asserts geometrically that the planar curvature of the interface is proportional to the height above the (uniquely determined) level at infinite distance from the plates (see, e.g., Theorem 2.1 of ). We assume in this work that the proportionality factor , as occurs for a non-zero gravity acceleration g directed downward toward the fluid. Physically, , where ρ is density change across the interface and σ is the surface tension arising from the fluid/fluid interface.
We are interested in categorizing the ranges of qualitatively distinct behavior that can occur. In accord with engineering practice and in cognizance of relevant uniqueness properties, the distinctions are best displayed in terms of non-dimensional parameters: setting , , , where 2a is the distance between the plates, (1.1) takes the form
In the non-dimensional coordinates, the plates are always two units apart. The physical concept of plate separation is replaced by the magnitude .
We proved in  that for arbitrary , there is a unique solution of (1.1) in the interval between the plates, achieving prescribed inclinations , (equivalently, contact angles , ) on the respective plates. Correspondingly, there is a unique solution of (1.1+) for arbitrary and contact angles. We obtain an explicit representation in terms of ψ as parameter by rewriting (1.1+) in the form
The latter relation separates:
on any interval on which ψ is monotonic in ξ. Note that (1.5) relates three distinct points on the solution curve, any two of which may coincide. We shall have to take pains in each instance to ensure that the correct branches for the roots are used. One sees easily that aside from the trivial solution of (1.1+), inflections of a solution curve occur exactly at crossing points of that curve with the ξ-axis, that at most one such point occurs on any solution, and that the sense of monotonicity of ψ reverses at every such point. Thus when such a crossing occurs, the integral in (1.5) must be split into two parts with the senses of integration reversed. We observe that only a single inflection can occur on any solution curve, see the assertions (1) to (5) in Sec. II of .
Solutions in the infinite intervals exterior to the plates are uniquely determined by the contact angle on the plates facing the respective interval and by the requirement of being defined in an infinite interval; see Theorem 2.1 in . These solutions are asymptotic to the ξ-axis but do not contact it; they admit the representations
Here refers to the plate position, and is the inclination on that plate. Note that four distinct solutions appear, depending on the signs of the non-constant roots. The constant root is taken as positive.
In Figures 1, 2, 3, particular integral curves T, I, II, III, IV, IV0, and V are sketched for the subset of solutions of (1.1) in the interval between the plates and meeting in the prescribed angle , and for successively decreasing plate separations. We have chosen for convenience . In the figures, some of the curves are extended beyond the plates as solutions, at least to the extent to which they can be represented as graphs. The sketched curves serve as barriers for distinguishing the qualitative structures of general solutions. No two of the curves in can cross each other within the interval in which both are graphs, and the regions between adjacent barriers serve to distinguish specific global behaviors.
Figure 1. Large plate separation.
Figure 2. Intermediate plate separation.
Figure 3. Small plate separation.
In the ensuing context, we discuss in detail the specific roles of the indicated curves. For convenience, we provide here a preliminary outline of essential features.
• The curves I and V are determined by (1.6) with and choosing respectively the positive and negative signs for the non-constant roots. IV0 is the unique curve of meeting at its crossing point with the ξ-axis. These three curves are rigidly attached toand are independent both of the position of, and of the contact angle.
When the plates are sufficiently far apart (B large enough), neither IV0 nor V extends to , and IV lies above both these curves between the plates. As the plates come together, a critical separation is attained, for which IV passes through and can then be shown to coincide with IV0. Figure 1 prevails when the plates exceed this separation so that ; we refer to such configurations as large separations. Closing the gap further, IV moves below IV0, remaining at first above V, and the relevant picture becomes Figure 2. Further gap closure leads to a value at which IV and V coincide, again with the common curve extending exactly to . For , V will lie above IV in the interval connected to in which both curves are graphs, and Figure 2 must be replaced by Figure 3. From (1.6) we find that the critical separation for this change is
using positive roots. The corresponding value is determined by (3.1). The two crucial values for are illustrated in Figure 4.
Figure 4. The half plate separationsand.
We established in  and in  that solution curves joiningandare attracting when the extended curve has a positive minimum or negative maximum. Denoting that height by , and the (horizontal) attracting force by F, we are led to a non-dimensional attracting force
When the extended curve is asymptotic to the ξ-axis at infinity, then and there is no force between the plates. The final alternative is that the extended curve meets the axis at some point . When that occurs, the curve is either the trivial solutionor else it crosses the axis in an angleat a uniquely determined pointand there will be a repelling force
We consider the family of solutions in the interval between and . We have from (1.4) that if , then the solution will not cross the ξ-axis; if , the curve will then attain a minimum height at a point where ; thus at this point we find ; in further continuation the curve becomes vertical at a height for which . If in (1.5) we set ξ to be the position of the plate , then we find
For given ‘dimensionless plate separation’ , relation (2.1) uniquely determines the (positive) heightfor which the solution meetsin the contact angleandin the prescribed angle. The valuethus found is the highest initial positionfor which the solution inextends to meet the plate. We designate this solution with T. It forms an upper barrier for all solutions in that extend as graphs to meet .
The dependence of the plate heights , on is illustrated in Figure 5.
Figure 5. The intersection heights,with the plates, as functions of half plate separation;.
Referring to any of Figures 1, 2, 3, we see that the two plates together with T and I determine a non-null closed region , topologically a disk, which is simply covered by a subset of solutions in , all of which yield attracting forces (or zero force in the unique case of I). We obtain such a region for every choice of ‘plate separation’ .
If , then for the givenno further attracting solutions can be found; all further solutions inwill be repelling. But if , then IV and V coincide, providing a negative solution yielding zero force; if , then IV moves below V and there is a new region of negative solutions providing attracting forces (Figure 3). Thus when the plates are close enough to each other, two complementary regions appear, one of positive solutions and the other of negative solutions, both of which yield attracting forces between the plates.
The curve II is the unique element of whose height vanishes on . The region between I and II is again simply covered by solutions in . All these curves lie in the upper half-plane within but intercept the ξ-axis when extended, and thus provide repelling forces.
III is the symmetric solution, achieving on the contact angle . This curve has special properties as we shall see below. Within the region , all curves of are repelling; all curves cross the axis and the sense of monotonicity of ψ in x reverses on crossing, so that special precautions must be taken in the representations. Corresponding comments apply for when . In that event, there are no solutions in below IV, and IV then plays to some extent the role at the bottom that T plays at the top, the adjacent solutions being, however, repelling rather than attracting.
3 Barrier curves
3.1 The barriers T and IV
By definition, for the upper curve T, we have ; equivalently, . Equation (2.1) determines the height . The counterpart for negative solutions is the curve IV at the bottom, which needs a bit more discussion. We introduce the further barrier IV0 which meets in angle at the level . This curve becomes vertical at a critical ‘dimensionless separation’ , and we find
If , then IV contains both positive and negative heights, and account must be taken of the change in sense of its curvature at the crossing point with the ξ-axis, where . Denoting the inclination at that point by , we obtain using (1.3) separately on the negative and positive portions of the curve,
3.2 The barriers I and V
3.3 The barrier II
3.4 The barrier III
The barrierIIIhas the unique property, that ifis fixed and the separation, the configuration remains repelling, with III asymptotic to the symmetric linear segment inclined at to the axis and joining the plates.
4 Force calculations
We proceed to calculate the forces between the plates in varying configurations. In practice, the accessible parameters will generally be the contact angles and on the sides of the plates facing each other, and the dimensionless plate separation . Other parameters, such as the height of the contact points and with the plates or the height of a local extremum or position of the crossing point with the ξ-axis, can be substituted via relations (1.4)-(1.6). Our basic force relations are (1.8) and (1.9), corresponding respectively to the attracting and repelling cases. The results for varying configurations are illustrated in Figures 6-13.
Figure 6. The anglefrom (4.2).
Figure 7. The net attracting force, from (4.4), with. The force vanishes when .
Figure 8. Attracting forces, negative solutions: plots of (4.5) and (4.6) assuming. Solid lines indicate forces at the specified (right vertical axis). Dashed line indicates (left vertical axis).
Figure 9. Repelling forces, positive solutions: plots of (4.2), (4.7), (4.8),. Solid lines indicate forces at the specified (right vertical axis). Dashed lines indicate specified angles (left vertical axis).
Figure 10. Equations (4.7) and (4.9);.
Figure 11. Equation (4.10);.
Figure 12. Forces on the symmetric curve III.
Figure 13. Plots in rangebased on (4.18) and (4.19). Note that for , and hold.
4AP Attracting forces, positive solutions
These are encountered only in (see Section 2). As shown in , the dimensionless force is determined, in units of σ, from
In view of (1.8), we may rewrite this relation in the form
which determines ℱ uniquely in terms of the contact angles on the plates and the separation.
4AN Attracting forces, negative solutions
If (see (1.7)), there are no such solutions. When the plates are close enough so that , a new region appears (Figure 3) in which the (extended) solutions have negative maxima, and the net force will again be attracting. We may emulate the 4AP discussion. The crossing angle of V with is now determined by
using the positive root.
4RP Repelling forces, positive solutions
Repelling solutions all cross the ξ-axis and thus change sign; however, we may characterize those that are positive between the plates as those lying in the region . Denoting by the inclination of II at the crossing point with , we find by procedures analogous to those above
4RPN Repelling forces, changing sign
In the region , we have . Solutions continue to repel, however their heights change sign between the plates. We shall see below that this has significant effects on limiting behavior as . Since the orientations change at the crossing points, we must split the integration in (4.8) into two parts. Applying (1.5) separately over each of the two segments and into which the crossing point divides the interval 2a and adding, we find
4S The symmetric curve III
This curve has a special interest. Regardless of plate separation, it crosses the axis at the midpoint between the plates, and thus yields repelling force for every separation. Relation (4.9) simplifies to
From (4.11) follows
leading to bounds in both directions for the repelling force . These bounds are precise asymptotically as . For large separations, (4.13) loses in precision. We improve it by observing that the actual solution in the positive portion of the interval between the plates lies below the line segment joining the crossing point to the height on . Thus,
yielding a perfunctory but conceptually useful bound for ℱ that could be improved in detail by using again the left-hand side of (4.13). Actually, the force for large a vanishes exponentially in a as follows from the general estimates of Siegel .
If B is small, we find from (4.13) and the monotonicity of ψ in x on each side of the halfway point between the plates that asthe inclination ofIIIbetween the plates tends uniformly to. As a consequence, we find that the normalized repelling force of the solutionIIItends in the limit to the magnitude
The above force calculations apply for any choice of separation 2a. To continue with solutions joining the plates but situated below III, we distinguish cases according to the plate separation.
In this event IV lies above IV0, see Figure 1. The range of inclinations achieved in the corresponding region is . All solution curves cross the x-axis between the plates, and the force calculation proceeds as in (4.9), (4.10), with the extended range for . We shall see, however, in Section 5 that from the point of view of limiting behavior as the plates approach each other, it would not be appropriate to join this region with the preceding one. No solutions meeting both plates and achieving the contact angleonexist belowIV.
IV lies below IV0 but above V, see Figure 2. IV0 meets in an inclination , determined by
In the region , we have , the solution curves cross the x-axis between the plates, and we may again use the procedure indicated by (4.9), (4.10). In the remaining region , there is no crossing point between the plates, the fluid level is negative with ψ decreasing in ξ, and the force is obtained from the modified version of (4.8):
Now IV lies below V, see Figure 3. We obtain a region which falls in the range of 4.2 above, then a region yielding repelling solutions with no axis crossing between the plates. lies in the range , where is determined from
5 Limiting behavior for small separation
With given contact angles on the two plates (corresponding to prescribed materials), we investigate the consequences of varying the separation of the plates. We effect the change conveniently by holding fixed and displacing in either direction. It is crucial to observe that in such a displacement, the barriers I, IV0 and V are rigidly attached to (and hence remain fixed), and the barrier III continues to pass through the midpoint on the x-axis. II and IV move downward as a decreases. The set of solutions examined is rigidly attached to and does not change as is shifted; however, the choice of elements within must change to maintain prescribed conditions on . Note that the geometric locus of II - when considered as an element of determined by its contact angles with the plates - moves upward in the family with decreasing a, but when considered as defined by its property of passing through the intersection of with the x-axis, it moves downward.
We distinguish the initial ℛ-regions and examine what happens to a typical solution curve in each such region, with decreasing a.
5.1 Curves above T
In a given configuration, there are no solutions above T that meet in angle and extend to . If we allow a to decrease with and fixed, a new T+ appears, lying above T. There are no solutions above T+, but the original is extended with the new (attracting) solutions in the region , consisting of solutions that previously did not extend to . Every curve above T and meeting in angle eventually falls into this category, as a decreases toward zero.
5.2 Curves above I
We consider a particular such curve of , displace toward , and ask what must be done to preserve the original angle on . Since is convex upward, the angle with will increase when is displaced toward . To retain the original angle , one must move to a curve above the original one. As a consequence, every curve oflying aboveImoves upward and remains attracting following the displacement. Additionally, new attracting curves will appear above T, as noted in 5.1.
Each of the curves considered has a positive minimum , or else achieves (in the particular case of I) a minimum zero at infinity. We can determine the net attracting force by estimating and using (4.1). Adapting (1.5) to the configurations considered, we find
Here can be arbitrary in the interval , where is the inclination of the barrier I at its intersection with . We distinguish the cases in which the minimizing point occurs outside the interval between the plates () from those for which it occurs within that interval (). In the former case, we have , and by inserting the end values or into the root under the integral sign, we get
The normalized attracting force ℱ can now be estimated using (4.1). Relations (5.2a), (5.2b) provide an explicit version of Laplace’s discovery  that the attracting forces remain attracting and become unbounded as the inverse square of the distance between the plates, as the separation decreases to zero.
5.3 Negative attractors
In the event , a second interval of (negative) attracting solutions appears above IV and below V. For in the range we have chosen, the discussion for these solutions is analogous and somewhat simpler than the one just given, as in no case does the (negative) maximum appear between the plates. Again attracting solutions remain attracting as plate separation decreases; the estimate (5.2a) prevails, albeit with and interchanged.
This case is discussed in explicit detail in ; we include here in outline form some essential features, returning for explicit convenience to direct physical notation. To begin, let us look at the point on , as in Figure 4. When is displaced an amount δ toward , the horizontally displaced will encounter too large an inclination from the element of passing through that point, as that element will have the same inclination on as does the indicated solution , and at every position x between the plates, its height u is smaller than that of that solution and thus by (1.1) its curvature is smaller. Thus the inclination of the field element at exceeds that at q, which in turn exceeds that at . Therefore in order to attain the initial slope again, one will have to move upward on the displaced plate . Since , and high enough points on yield , a (unique) such point can always be found.
We observe now that on the original vertical segment of joining I and II, holds (see Figures 1, 2, 4 for notation). Thus the angle is attained at some intermediate point on I between and . We choose δ so that passes through that point. The solution curve then has the identical data on the two vertical plates as does I, and by the uniqueness theorem (see ) must coincide with I.
Looking more closely, we see that by movingcontinuously toward, we obtain a continuous family of solutions joining the plates, with left-hand end points rising in the motion, and such that at some intermediate position strictly between the plates, the solution will coincide with the portion ofIto the right of that point. Once that happens, all further motion ofto the right leads to attracting solutions to which the material of preceding sections applies.
Note that for the given on , this ‘crossover’ behavior occurs for all in the range between and . The ‘crossover position’ between the plates, where the solution curve joins with I and yields zero net force, is determined explicitly from the relation
For any chosen within its range, we obtain a picture of a succession of solution curves, defined over successively shorter subintervals attached to between the plates, each successively closer to I and yielding successively smaller repelling forces, until the solution curve coincides with a non-null portion of I, providing zero force. See Figure 14.
Figure 14. Behavior of solution curves with changing plate separation; contact angles prescribed.
These are still repelling solutions as they continue to cross the x-axis. Nevertheless, there are significant changes from the case just considered, as the initial heights on are negative. The sense of curvature of the solution curves reverses in the negative region, and account must be taken of that change in two senses.
We note first that the range of angles that arises is . To determine the force arising from a given solution curve, we need only determine according to (1.9) the angle of intercept with the x-axis. Using (1.5) separately in the positive and in the negative regions and adding, we obtain
A further change in behavior occurs in that if one moves a small distance δ toward , one finds that to maintain the same initial angle , one must look downward instead of upward as before. As a consequence, the new solution curve lies below the previous one, it will be further from I than was the previous one, and will increase the repelling force rather than decreasing it as above.
Nevertheless, it turns out that the sequence of solutions thus constructed converges to a segment ofIjust as did the previous one. This assertion may at first seem in conflict with the behavior just described; however, one can show that although the solution curves at first diverge from I, their starting points on the intersections with the successive planes actually rise and become positive prior to reaching . Once that happens, the discussion of 5.4 applies again without change, and the corresponding behavior is observed. A complete proof of this behavior appears in .
We see that if the initially chosenlies in the range, then, as the plates are brought together, the repelling force will initially strengthen to a maximum, and then will weaken to a critical separation at which the solution coincides withIand yields zero force, and will finally become attracting with force increasing as the inverse square of the separation, according to (5.2a), (5.2b).
We can characterize these critical configurations explicitly. The maximum repelling force is achieved corresponding to a starting point lying on the x-axis with solution inclined at angle . We thus set , in (1.4), and in (1.5) we let denote the coordinate of the crossing point. We find
5.6 The symmetric curve III
As we move downward through the range , the angle increases from the angle with which II cuts the x-axis to the angle , achieved by III. The curve III itself fails, however, to become attracting with decreasing separation; we see that immediately since, due to its symmetry, it cuts the x-axis for every separation. In Section 4S we have already established upper and lower bounds for the repelling force in this case, notably the non-zero limit as the plates approach each other. The material above together with what is to follow shows that IIIis isolated in this respect; every other configuration with fixed angles on the plates becomes attracting as the plates come together, of magnitude rising to infinity as the inverse square of the separation distance. Thus there is a very striking singular limiting behavior in configurations adjacent to the symmetric one. Physically, this corresponds to liquid going to positive infinity when it is initially above III, and to negative infinity when it is initially below III. It should be of interest to observe this transition experimentally.
We continue to discuss the remaining cases that occur; to this purpose we return to non-dimensional notation.
This is illustrated in Figure 1. A new region appears with in the range . Since all these angles exceed , they cannot reappear on the curve I as happens for initial datum above that of III, and thus the convergence to a segment of I does not recur here. We observe that the range of that appears is identical to the range of ψ on the portion of V between the plates. Since ψ is monotonic on V, there is a unique point on this arc at which the initially chosen value of appears, see Figure 4. Denote by the x-coordinate of that point.
When is moved toward , one finds one must move downward from the initial height in order to achieve again the same initial inclination . Thus the succession of solution curves moves toward V, with the repelling force decreasing. The exotic behavior noted in 5.5 above, with repelling force initially increasing as the plates come together, does not reappear for the region belowIII.
When is situated at , the data of the relevant solution curve at its two endpoints coincide with those of V at those points, and thus the two curves coincide on the interval, with vanishing force. Further approach of toward yields attracting solutions, with forces controlled by (5.2a), (5.2b).
The relevant picture for the initial configuration is now Figure 2. We obtain two new regions for repelling solutions, viz. and .
All solutions are repelling and cross the axis between the plates. The configuration is fully analogous to that of , and analogous considerations apply. See Section 5.5. The repelling forces successively increase to a maximum of , then solutions move to V and proceed to cross over and become attracting.
All solutions are repelling and cross the axis outside the plates. The situation is essentially that of . As the plates move together, the repelling force decreases monotonically to zero and then attracting forces prevail. See Section 5.4.
The situation is now essentially analogous to the initial discussion for curves lying above III. We remark the technical distinction that the minimizing point on the upper barrier arc T lies always between the plates; those for the corresponding lower barriers lie to the right of , although they approach that plate with decreasing separation 2a.
6 Notes added in proof
1. After completing this work, we were informed by John McCuan of an earlier paper  in which some of the material relates closely to the topic of the present study. Our contribution can be regarded as an improvement on Section 4 of , in the sense that we study the question in the context of the fully nonlinear equations, in preference to the linearization adopted in that reference. The particular geometry of the configuration permits us to integrate the equations explicitly in original form, yielding expressions that describe general physical laws. Beyond the evident improvements in precision and detail, we were led to the discovery that the net attracting (repelling) force on the plates is independent of the contact angles that occur on their outer sides; thus the restriction made in  to plates with identical angles on the two sides is superfluous. We find also the general theorem that the net force is repelling or attracting, according as the (extended) solution curve joining the plates in a vertical section does or does not contain a zero for the height on its traverse, the net force being then provided respectively by the elementary formulas (1.9) or (1.8). We obtain additionally a more complete description of the limiting behavior as given plates approach each other (this behavior becomes dramatically singular for solutions close to the symmetric one; see Section 5 of the present work).
2. The exact formal theory was additionally a help for us toward avoiding misleading inferences suggested by the linearization, among them the erroneous statement in  opening the final paragraph on p.819: ‘This result shows that vertical plates…will attract if they have like menisci and otherwise repel…’. In fact (as shown in Section 4AP) for any plate separation and acute angle , the solutions in the non-null subset of for which have unlike menisci at the plates and for these solutions the plates nevertheless attract each other.
The authors declare that they have no competing interests.
Both authors contributed equally in this work, in all respects.
The latter author is indebted to the Mathematische Abteilung der Universität and to the Max-Planck-Institut für Mathematik in den Naturwissenschaften, in Leipzig, for invaluable support during preparation of this work.
Finn, R, Lu, D: Mutual attractions of partially immersed parallel plates. J. Math. Fluid Mech.. 15(2), 273–301 (2013). Publisher Full Text
Siegel, D: Height estimates for capillary surfaces. Pac. J. Math.. 88(2), 471–515 (1980). Publisher Full Text
Laplace, PS: Traité de mécanique céleste, oeuvres complète, vol. 4, Gauthier-Villars, Paris, 1805, Supplément 1, livre X, pp. 771-777. Supplément 2, livre X, pp. 909-945. See also the annotated English translation by N. Bowditch 1839. Chelsea, New York (1966)
Vella, D, Mahadevan, L: The cheerios effect. Am. J. Phys.. 73, 817–825 (2005). Publisher Full Text