cellular structures@xcite appear in a wide range of natural phenomena , and have puzzled and fascinated scientists for decades@xcite .
they can be generally described as packings of space - filling cells of roughly polygonal shape , separated by thin interfaces to which a surface energy is associated .
they arise as a consequence of competition between domains under the constraint of space - filling , and there is often some mechanism , such as the migration of a conserved quantity across the interface , which makes them evolve in time , i.e. coarsen . +
foams such as those obtained by shaking soapy water are the examples of cellular structures closest to our everyday experience .
two - dimensional foams may be obtained by confining a soap froth between two closely spaced glass plates . in spite of their apparent simplicity , foams display much of the phenomenology appearing in coarsening cellular structures .
foams have been the subject of interest since their relevance in the problem of grain growth was pointed out by smith@xcite . despite the attention they have received , the understanding of their dynamical properties proved to be a tricky problem .
even some of the most basic issues , such as their asymptotic scaling properties , has been a matter of debate until recently ( see references above ) .
numerically exact simulation procedures@xcite , as well as analytical@xcite and numerical@xcite approximations , have been useful in understanding the dynamics of ideal foams , but the system sizes accessible to available simulation procedures are strongly limited ( for a complete account see refs.@xcite ) . + a satisfactory theoretical framework for the description of foams has not yet been achieved .
the main result of this work consists in establishing a _ rigorous connection _ between foams and partitions(vp ) .
this connection provides a set of fundamental degrees of freedom for the foam ( the source s locations in space ) and therefore constitutes a step towards the above mentioned theoretical understanding . on the other hand , an efficient method for the numerical simulation of ideal foams
would also be highly desirable , and we propose that such a method could be obtained by exploiting the equivalence between foams and vp reported in this work .
more precisely , we establish the following _ correspondence _ between equilibrated two - dimensional foams ( ef ) and a generalization of vps , the sectional multiplicative partition ( smvp ) : + _ given an arbitrary ef , it is always possible to find sources in space such that a smvp with respect to them exactly gives this foam . _ + in section  [ sec : tessellations ] the vp and its generalizations are reviewed .
the simplest generalizations of the vp concept are the sectional partition and the multiplicative partition .
these correspond to adding a constant to the ( square ) distance and to multiplying it by a constant respectively . at the end of this section a combination of
both is introduced , the sectional multiplicative partition .
we will use this partition in order to describe foams .
+ the demonstration of the above mentioned correspondence between foams and smvps is divided in two parts ( sections [ sec : smvprecognition ] and [ sec : equilibrium ] ) for clarity . in section  [ sec : smvprecognition ] , the _ recognition problem _ for smvps is solved .
the recognition problem consists in giving the sufficient geometrical conditions that a circular partition has to satisfy in order to be a smvp .
we will see that if a circular partition has all of its vertices _ aligned _ and admits an _ oriented reciprocal figure _
, then it is a smvp . in other words
, it is always possible to find sources in three - dimensional space that give this circular partition as a smvp . in this section
we also describe the procedure to _ find _ the sources when these conditions are met .
some material that is needed in this section is described in the appendix .
+ section  [ sec : equilibrium ] deals with the equilibrium conditions for foams and its geometrical consequences .
we start by writing the force and pressure equilibrium conditions in compact form in section  [ sec : eq_equations ] . in sections
[ sec : eq_align ] and [ sec : eq_reciprocal ] it is shown that an equilibrated foam has all of its vertices _ aligned _ and admits an _ oriented reciprocal figure _ respectively .
therefore all equilibrated foams satisfy the conditions required in section  [ sec : smvprecognition ] , and can be described as smvps . + finally in section  [ sec : conclusions ] the implications of this result are discussed , and some perspectives for future work are advanced .
given a set of @xmath2 point - like _ sources _ @xmath3 in @xmath4-space , the _ partition _ ( vp ) or _ tessellation _ of space with respect to @xmath3 is a classification of space into _ cells _
@xmath5 defined such that @xmath6 if @xmath7 is closer to @xmath8 than to any other source  @xcite .
@xmath9 where @xmath10 denotes the distance between @xmath7 and source @xmath8 .
this construction is also known under the names : wigner - seitz cells , dirichlet tessellations , thiessen polygons . in two dimensions , the dual lattice of a partition
is called delaunay triangulation .
cell @xmath5 can be seen as the `` region of influence '' of @xmath8 , were the sources competing for some spatially distributed resource .
it is common to use these constructions , and their generalizations , as _ approximate _ models for cellular structures occurring in nature@xcite . +
two neighboring cells @xmath5 and @xmath11 are delimited by an _ interface _ @xmath12 of points @xmath7 equidistant from @xmath8 and @xmath13 . @xmath14
these interfaces are @xmath15-dimensional hyperplanes normal to @xmath16 . in two dimensions @xmath5 are convex polygons and the interfaces @xmath12 are straight lines ( fig .
[ fig : voronoi_partition ] ) .
three interfaces @xmath12 , @xmath17 and @xmath18 meet at a _
@xmath19 , which is equidistant from @xmath8 , @xmath13 and @xmath20 .
this means that @xmath19 is the center of a circle through @xmath8 , @xmath13 and @xmath20 .
vertices of higher multiplicity are possible for particular locations of the sources .
for example a fourfold vertex would exist if four sources lay on the same circle .
we will ignore the existence of these particular configurations , that is we will assume _ generic _ source locations . under this assumption ,
all vertices are of multiplicity three .
+ the concept defining partitions is _
equidistance_. a simple way to generalize them is changing the way in which distances to the sources are measured .
one lets each source measure the distance to points @xmath7 according to its own `` rule '' .
the interface @xmath12 is then the set of points @xmath7 for which two neighbors `` claim '' the same distance .
two simple ways to do this are : adding an arbitrary constant to the square of the distance ( sectional partition  @xcite ) , and , multiplying the distance by a constant ( multiplicative partition  @xcite ) .      a _ sectional partition _ ( svp )
is defined  @xcite as a @xmath21-dimensional cut of a higher dimensional partition .
the sources @xmath22 of the sectional partition are defined as the projections of the original sources @xmath3 onto this lower dimensional @xmath21-space .
cell @xmath5 associated to source @xmath0 is the intersection of cell @xmath23 ( associated to @xmath8 ) with this @xmath21-dimensional hyperplane . for example , take sources @xmath3 in three dimensions ( fig .
[ fig : svpview ] ) and construct a vp with three - dimensional cells @xmath24 and plane interfaces @xmath25 .
now project the sources @xmath26 onto the @xmath27 plane ( which we call the @xmath28 plane ) to obtain the projected sources @xmath29 , which have associated `` heights '' @xmath30 .
assign to each source @xmath0 the intersection @xmath5 of @xmath23 with @xmath28 .
this defines the a sectional partition of @xmath28 with respect to sources @xmath22 with heights @xmath30 .
cells @xmath5 are thus defined as @xmath31 where @xmath10 is the distance , on the @xmath28 plane , between @xmath0 and point @xmath7 .
+ in the same way , interfaces @xmath12 are defined as @xmath32 properties of these partitions in two dimensions are :    * interfaces are straight lines and ( generically ) meet at triple vertices .
( in dimension @xmath33 , vertices have multiplicity @xmath34 ) . * interface @xmath12 is normal to @xmath1 , but not in general equidistant from @xmath0 and @xmath35 . *
the partition is unchanged if all @xmath30 are changed according to : @xmath36 for @xmath37 an arbitrary real .
an important difference with the partition , is that in this case there may be sources to which no cell is associated .
this happens if @xmath23 is not cut by @xmath28 .
it is easy to see that svps are _ equivalent _ to laguerre partitions  @xcite .
+ appendix  [ sec : rprecognition ] discusses the _ recognition problem _ for this kind of partitions , that is giving the sufficient conditions for an arbitrary rectilinear partition to be a svp .
we will use similar concepts in section  [ sec : smvprecognition ] in order to solve the recognition problem for sectional multiplicative partitions . + as an example of the application of sectional partitions , consider the case of sources @xmath8 whose location in space changes in time , giving rise to a time - dependent partition of @xmath28 .
if for example source @xmath38 moves away from @xmath28 , the two - dimensional cell @xmath39 will shrink and finally disappear when @xmath40 no longer cuts @xmath28 . therefore the number of cells can vary without changing the number of sources .
thus sectional partitions can be used as dynamical models for crystal growth  @xcite , and other processes in which some cells disappear or are created .      the multiplicative partition ( mvp )  @xcite is defined by assigning to each source @xmath8 a positive _ intensity _ @xmath41 , and defining the _ multiplicative distance _ @xmath42 .
the cell @xmath5 associated to source @xmath8 is then the set of points @xmath7 closer to @xmath8 ( in terms of this multiplicative distance ) than to any other source @xmath43 interfaces @xmath12 are hyperspherical surfaces ( circle arcs in two dimensions ) satisfying @xmath44 in two dimensions , the circular interface of a mvp with two sources is one of the apollonius circles ( see for example  @xcite ) of those two points . in arbitrary dimensions
let @xmath45 be the distance between two sources @xmath8 and @xmath13 , and @xmath46 . without loss of generality
we take @xmath47 .
this means that @xmath13 has the smallest intensity and therefore @xmath11 will be the interior of a hypersphere , while @xmath5 will be its exterior .
the following properties are satisfied in any dimension .
* @xmath8 is contained in @xmath5 .
* @xmath12 has radius , and its center @xmath48 is on the straight line going through @xmath8 and @xmath13 . *
@xmath48 is located at a distance @xmath49 from @xmath8 , that is , it never lays between @xmath8 and @xmath13 .
three sources give rise to a tessellation like the one shown in fig .
[ fig : mvp_vertex_2d ] .
the exterior of the two `` bubbles '' is the infinite cell associated to the source with the largest intensity , @xmath38 in this example .
for some choices of @xmath50 the interfaces will not intersect . in this case
the tessellation is simply a pair of circles , each containing one of the sources with smaller intensities , while the exterior of these two circles is the cell of the source with larger intensity .
mvps can be interpreted again as `` areas of influence '' of sources @xmath8 , but now each source has a different strength . + consider a mvp of two dimensional space with respect to a set of sources @xmath3 . if three interfaces @xmath12 , @xmath17 and @xmath18 meet at a vertex ( again , vertices of higher multiplicity are only possible for particular configurations , which we ignore here ) , then centers @xmath48 , @xmath51 and @xmath52 lay on a line .
the reason why this is so is simple . if two interfaces @xmath12 and @xmath17 have a common point @xmath19 , then their _ continuations _ must meet again at another point @xmath53 , the _ conjugate vertex_. but then the third interface @xmath54 must also go through this point since @xmath55 and @xmath56 implies @xmath57 . therefore
two points exist ( @xmath58 and @xmath59 ) at which all three circles intersect
. then the centers of these circles must be on a line ( fig .
[ fig : mvp_vertex_2d ] ) .
+ let us define a _
circular partition_(cp ) of two - dimensional space to be a classification of space in cells delimited by arbitrary circle arcs that meet at triple points called vertices .
no three centers of these arbitrarily defined interfaces will in general be on a line .
we will say that a vertex of a cp is _ aligned _ if the centers of the three interfaces defining it are on a line .
as we saw , all vertices of a mvp are aligned in 2d .
therefore for each triple vertex of a mvp in 2d , the three sources and the three centers form the configuration  @xcite @xmath60 of projective geometry . since @xmath48 can not be between @xmath0 and @xmath35 , the three centers can only be on one of the two external segments of the configuration .
+ in a three - dimensional mvp , four cells meet at each vertex @xmath61 , giving rise to six spherical interfaces @xmath62 .
the centers of these interfaces are aligned in triplets so that the six centers also form the configuration @xmath60 , this time in three - dimensional space .
this has the implication that the six centers are necessarily on the same plane .
+ the mvp was first introduced by boots  @xcite to describe areas of influence in geography .
ash and bolker  @xcite also studied the recognition problem , that is , under which conditions a cp is a mvp .
the visual similarity between this kind of partitions and two - dimensional foams is striking , and suggests the idea to find a connection between them .
obviously the centers of the films of a foam must be aligned for each triple vertex , if the foam is to be described as a mvp , since mvps are aligned .
one finds  @xcite that this alignment condition is indeed satisfied by all vertices of arbitrary equilibrated foams . despite this ( which is a necessary but not sufficient condition for a foam to be a mvp )
, two - dimensional mvps can not describe all possible equilibrated foams in two dimensions .
this work shows the reason of this limitation : in order to describe arbitrary equilibrated two - dimensional foams we must introduce the sectional variant of a multiplicative partition . in other words , instead of confining the sources to the plane , we let them exist in a three - dimensional space and obtain the foam as _ a two - dimensional cut of a three - dimensional mvp_.      we will restrict the description to the case of a plane cut of a three dimensional mvp .
the generalization to a @xmath21-dimensional cut of an @xmath4-dimensional mvp is straightforward .
+ a plane cut , with a plane @xmath63 , of a three - dimensional mvp defines a sectional multiplicative partition ( smvp ) of @xmath28 ( see figs .  [
fig : smvp_vertex ] and [ fig : smvp_upper ] ) .
the sources @xmath0 of this smvp are the normal projections onto @xmath28 of the original sources @xmath8 , and may be seen to have as attributes both an intensity @xmath41 and a height @xmath30 .
cells @xmath64 associated to sources @xmath0 are defined as @xmath65 in the same way , interfaces @xmath12 are circle arcs satisfying @xmath66 an example of a smvp with three sources is shown in figs .
[ fig : smvp_vertex ] and [ fig : smvp_upper ] .
we see that there are two vertices on @xmath28 at which the three interfaces meet . in a general case ( for example in a partition with respect to many sources like in fig .
[ fig : mvp2d - many ] ) , if vertex @xmath19 exists , then the continuations of interfaces @xmath12 , @xmath17 and @xmath18 meet at a _ conjugate _ vertex @xmath53 also .
this implies the alignment of centers @xmath48 , @xmath51 and @xmath52 , which could also be deduced from the alignment of @xmath67 , @xmath68 and @xmath69 in three - dimensions ( fig .
[ fig : smvp_vertex ] ) .
+ we notice that the smvp is equivalent to a multiplicative laguerre partition , since a svp is equivalent to a laguerre partition  @xcite .
as is usual in sectional partitions , in the smvp there may be sources @xmath0 with no associated cell , those for which the corresponding three - dimensional `` bubble '' @xmath23 does not cut @xmath28 .
+ it was seen in section  [ sec : mvp ] that the spherical interface @xmath25 in three dimensions is cut normally by the segment @xmath16 . as a projective consequence of this ,
the straight line containing @xmath0 and @xmath35 on @xmath28 is normal to the circular interface @xmath12 . in other words ,
@xmath0 , @xmath35 and @xmath48 are on the same line . furthermore and since @xmath67 never lays between @xmath8 and @xmath13 in 3d , we notice that @xmath48 is always outside the segment @xmath1 .
in this section the sufficient conditions for a cp to be a smvp are given . we first introduce the notion of oriented reciprocal figure for cps , and then proceed to demonstrate that an aligned cp admitting such a reciprocal figure is a smvp .
we will now generalize the concept of reciprocal figure as appropriate for circular partitions .
we say that a graph @xmath70 made of points @xmath0 connected by edges @xmath71 forms a reciprocal figure for a cp @xmath72 if ,    * points @xmath0 are in correspondence with the cells @xmath5 of @xmath73 .
* edges @xmath71 are in correspondence with the interfaces @xmath12 of @xmath73 .
* for each @xmath12 in @xmath73 , points @xmath0 , @xmath35 and @xmath48 lay on the same line .
consider two regions @xmath5 and @xmath11 separated by a circular interface @xmath12 .
a reciprocal figure @xmath70 will be said to respect _ orientation _
if for each @xmath12 in @xmath73 :    * @xmath48 is not between @xmath0 and @xmath35 . * starting from @xmath48 and traveling along @xmath74 to infinity , points @xmath0 and @xmath35 are found in the same order as regions @xmath5 and @xmath11 are .
we saw already that all smvps are aligned ( all vertices satisfy the alignment condition introduced in sec .  [
sec : mvp ] ) . on the other hand , it is clear that the sources @xmath0 of a smvp form a reciprocal figure that satisfies orientation .
therefore @xmath75 figure  [ fig : convexity ] shows four possible partitions that share the same centers @xmath76 , @xmath77 and @xmath78 . all four admit reciprocal figures .
but only in cases a ) and c ) a reciprocal figure that satisfies orientation can be constructed .
therefore neither b ) nor c ) can be a smvp .
the reason is that the orientation condition can not be satisfied if any of the vertices is not _ convex_. a vertex is _ convex _ if all the internal angles formed by the interfaces are less than @xmath79 .
convexity of all vertices is clearly a necessary condition for a circular partition to be a smvp .
+ we are now ready to give the _ sufficient _ conditions for a circular partition to be a smvp .
we will show in section  [ sec : find_sources ] that if a circular partition @xmath73 is aligned and admits a reciprocal figure that satisfies orientation , then it is a smvp .
@xmath80      given a circular partition @xmath73 of a region of two dimensional space , all whose ( triple ) vertices are aligned , we show here that , if @xmath73 admits a reciprocal figure that satisfies orientation , then @xmath73 is a smvp with respect to sources @xmath8 located somewhere above points @xmath0 .
our demonstration is constructive , that is we explicitly show how the sources and intensities are determined .
we will use for this purpose the inversion transformation ( appendix  [ sec : inversion ] ) to _ straighten _ each vertex in turn , that is , transform a circular vertex into a rectilinear one .
this rectilinear vertex is one of a of a svp .
sources are located in this _ straight _ representation ( appendix  [ sec : rprecognition ] ) , and then transformed back to the original system .
the intensities @xmath41 of the corresponding smvp are obtained in this back - transformation , since the smvp is invariant under inversion .
+ let @xmath58 be an aligned vertex on which three interfaces @xmath81 meet , and @xmath82 the three points of the reciprocal figure associated to the three bubbles sharing the vertex , as in fig .
[ fig : reciproca1 ] . the conjugate vertex @xmath59 is obtained as the intersection point of the _ continuations _ of the interfaces , which happens at a unique point because of alignment .
let furthermore @xmath83 be the straight line through @xmath59 and @xmath0 .
we start by showing that    [ th : recircles ] a monoparametric family of triplets of circles @xmath84 exists , which has the following properties :    * @xmath85 , @xmath86 and @xmath87 intercept each other on @xmath59 . *
@xmath88 . *
the intersection point @xmath89 between @xmath90 and @xmath91 lays on @xmath83 .
the non - trivial content of the theorem is the fact that the three intersection points @xmath89 of these normal circles @xmath90 lay on the lines @xmath83 .
for example take an arbitrary point @xmath92 on @xmath93 and draw a circle @xmath85 through @xmath92 , @xmath59 and normal to @xmath94 .
let @xmath95 be the intersection of @xmath85 with @xmath96 .
draw now a circle @xmath86 trough @xmath95 , @xmath59 and normal to @xmath97 .
let @xmath98 be its intersection with @xmath99 .
then the circle @xmath87 through @xmath59 , @xmath92 and @xmath98 is normal to @xmath100 .
+ * proof . *
the demonstration is done by first performing an inversion with center @xmath59 and arbitrary radius , whereupon the interfaces @xmath12 are transformed into straight lines @xmath101 meeting at @xmath102 ( fig .
[ fig : reciproca2 ] ) .
this transformed system will be referred to as the `` straight '' representation of that vertex .
primed names are used in the straight representation , with the exception of @xmath59 , whose original location is kept in fig .
[ fig : reciproca2 ] ( @xmath59 itself is mapped to @xmath103 by the inversion ) .
lines originally not going through @xmath59 are transformed into circles through @xmath59 , as is the case of @xmath104 , which form the reciprocal figure in the original system .
circles @xmath90 are transformed into straight lines @xmath105 going through @xmath106 and @xmath107 , respectively on @xmath83 and @xmath108 , which are invariant .
+ the inversion transformation preserves angles , therefore circle @xmath109 is normal to the now straight interface @xmath101 .
this means that its center @xmath110 must lay on @xmath101 .
consider the figure formed by @xmath102 , @xmath111 , @xmath112 and @xmath113 , and the six lines joining them .
maxwell  @xcite showed that a figure made of four points joined by six lines _
always _ has a reciprocal figure ( to see this just consider the centers of the four circles going through three of these points ) .
we will now identify it on fig .
[ fig : reciproca2 ] as the figure formed by @xmath59 , @xmath114 , @xmath115 and @xmath116 .
first we notice that circles @xmath109 and @xmath117 intercept each other at two points @xmath59 and @xmath118 , both on @xmath108 .
therefore lines @xmath119 are normal to @xmath108 , and we have identified the first three lines of the reciprocal figure .
the other three lines @xmath120 of the reciprocal figure have to be normal to the segments @xmath121 , so they are normal to the interfaces @xmath101 .
since the inversion preserves angles , this means that @xmath88 in the original system ( fig .
[ fig : reciproca1 ] ) .
@xmath122 + 0.2 cm we see that the existence of triplets of circles @xmath123 with the above mentioned properties is a consequence of the existence , in the straight representation , of a figure @xmath124 , which is the reciprocal of @xmath125 . in other words ,
theorem  [ th : recircles ] means that in the straight representation , triangle @xmath126 forms a reciprocal figure for the rectilinear vertex .
clearly this reciprocal figure satisfies orientation if the original figure satisfies orientation in the circular system .
therefore as discussed in appendix  [ sec : rprecognition ] , it is always possible to find three sources @xmath127 , located at heights @xmath128 above points @xmath106 , such that this rectilinear vertex is a svp with respect to them .
once these sources @xmath129 are known in the straight representation , a second inversion around the same point @xmath59 takes us back to the original system and provides the original locations @xmath8 of the sources .
+ if @xmath130 is the first vertex to be considered in the system , we have one degree of freedom in the determination of the circles @xmath90 , or equivalently of points @xmath89 .
once in the straight representation , there is one more degree of freedom : the determination of _ one _ of the heights @xmath128 . we will now show that , if this @xmath128 is chosen appropriately , the three back - transformed sources @xmath8 are located above their respective points @xmath0 . + consider a plane @xmath131 that is normal to @xmath28 and contains @xmath83 , as shown in fig .
[ fig : above ] .
draw a circle @xmath132 through @xmath89 , @xmath59 and normal to @xmath28 .
let @xmath133 be the intersection of this circle with the normal @xmath134 to @xmath28 through @xmath0 .
now under an inversion with center @xmath59 , point @xmath135 is above @xmath106 , since circle @xmath132 is transformed in a straight line @xmath136 , normal to @xmath28 .
straight line @xmath134 normal to @xmath28 through @xmath0 is now a circle @xmath137 through @xmath59 , @xmath135 and normal to @xmath28 .
its intersection with @xmath28 determines the image @xmath138 .
now we now that source @xmath129 is located somewhere on @xmath136 .
the sources can be displaced vertically ( simultaneously ) according according to what we saw in appendix  [ sec : rprecognition ] , but not independently . fixing the height @xmath128 of one of them determines the other two uniquely .
what we want to demonstrate is that if one of the sources @xmath129 coincides with its point @xmath135 , then the other two also do . in order to do this ,
take for example the sphere @xmath139 through @xmath59 , @xmath140 , @xmath141 and normal to @xmath28 . in the straight representation , interface @xmath142 of the partition with respect to sources @xmath140 and
@xmath141 is a plane , equidistant from @xmath140 and @xmath141 , and normal to @xmath143 .
this means that @xmath142 is normal to @xmath139 . as a consequence ,
the intersection @xmath144(a circle ) of @xmath139 with @xmath28 is normal to the intersection @xmath145(a straight line ) of @xmath146 with @xmath28 .
now assume that @xmath147 .
this means that @xmath148 since @xmath149 is in this case contained in @xmath139 .
but then the circle @xmath144 is the circle @xmath150 ( fig .
[ fig : reciproca2 ] ) that goes through @xmath151 , @xmath152 and @xmath59 and is normal to @xmath145 .
this implies that @xmath152 is also in @xmath139 , which in turn implies @xmath153 .
the same reasoning can be repeated for the pair of sources @xmath154 and @xmath155 .
+ we have thus shown that all sources @xmath8 are , in the back - transformed , original system , located above the points @xmath0 of the reciprocal figure .
in other words , that the given reciprocal figure is the projection of the sources @xmath8 on the @xmath28 plane .
+ now we have to determine the heights @xmath30 and intensities @xmath41 .
the heights @xmath30 of the sources can be found found by back - transforming the sources @xmath129 from the straight system .
but in practice this is not necessary . by looking at fig .
[ fig : above ] we notice that triangle is rectangle .
therefore @xmath156 which suffices to locate the sources @xmath8 with the sole knowledge of points @xmath0 and @xmath89 in the original system . knowing now the spatial locations of the sources @xmath8 in three dimensions , the intensities @xmath41 follow immediately since the vertex @xmath59 is equidistant  @xcite from the three sources .
this means @xmath157 here @xmath158 is an arbitrary positive constant .
alternatively we can get the same result by using the transformation properties ( eq .  [ eq : dist_transf ] ) of the @xmath41 s and the fact that the intensities are all equal in the straight system . +
we have used the inversion transformation to demonstrate that the given vertex is a smvp with respect to sources @xmath8 located above @xmath0 , and showed how the sources @xmath8 can be located without the need to explicitly perform an inversion for each vertex .
all we need is to know the location of points @xmath89 , which we are free to choose for the first vertex under consideration .
this first choice determines all subsequent @xmath89 points , and therefore all sources .
it is clear from fig .
[ fig : above ] that @xmath159 has to be satisfied .
therefore one has to choose the triplet of circles @xmath160 such that this condition is verified for the three points @xmath89 .
if one of the @xmath89 coincides with @xmath0 , this means that the corresponding height @xmath30 is zero , i.e. the source @xmath8 lays on the @xmath28 plane .
a point @xmath89 closer to @xmath59 than the corresponding @xmath0 is not acceptable .
+ the general procedure to determine the sources can now be described .
assume we are given a circular tessellation which satisfies the required conditions of alignment and convexity , and that we are able to find , or are given , a reciprocal figure @xmath0 for it .
we start from one of the cells of the partition , as in fig .
[ fig : bubblesources ] .
take an arbitrary vertex to start with , for example vertex @xmath161 , and determine tentative locations for point @xmath162 on @xmath163 .
this fixes @xmath92 and @xmath95 as discussed in section  [ sec : find_sources ] .
after fixing @xmath162 the rest of the construction is uniquely determined .
the locations @xmath164 determine @xmath165 and @xmath166 , and the corresponding intensities @xmath167 and @xmath168 , through equations ( [ eq : zi ] ) and ( [ eq : ai ] ) . now proceed to the neighboring vertex @xmath169 .
two of the sources , @xmath170 and @xmath171 are already known , only the height @xmath172 of @xmath173 has to be found .
first one has to find the _
locations of @xmath162 and @xmath95 when defined from vertex @xmath174 . in order to do this ,
lines @xmath175 from vertex @xmath174 are drawn , and @xmath162 and @xmath95 are located using the fact that the heights @xmath176 and @xmath166 are already known .
equation ( [ eq : zi ] ) implies    @xmath177    once the new @xmath162 and @xmath95 are known , the next step is to draw circles @xmath178 and @xmath86 through them and normal to the respective interfaces .
their intersection gives the location of @xmath98 ( this intersection always occurs on @xmath179 , as theorem [ th : recircles ] shows ) .
this determines @xmath172 and @xmath180 using equations ( [ eq : zi ] ) and ( [ eq : ai ] ) .
this procedure of triangulation is repeated until all sources are located .
eventually it may happen that a point @xmath89 is found to be closer to the conjugate vertex @xmath59 than the corresponding @xmath0 .
this is not acceptable and means that the _ tentative _ starting position of @xmath162 has to be changed .
it must be shifted away from @xmath181 and the whole procedure repeated .
it is easy to see that taking a large enough value of @xmath162 always solves this problem .
+ one could ask whether this construction can be closed self - consistently .
for example after determining @xmath182 from @xmath170 and @xmath183 in our example of fig .
[ fig : bubblesources ] , one can go on with the procedure as if @xmath38 were unknown .
would the position of @xmath38 determined by @xmath170 and @xmath182 be the same one as found initially ?
alternatively : if we used @xmath170 and @xmath38 to determine @xmath182 instead of going around the bubble in the opposite sense , would its position be the same as found after going around the bubble ?
the answer is yes because of unicity .
as discussed in appendix  [ sec : rprecognition ] , the height of one of the sources of a vertex determines the other two uniquely .
this means that @xmath184 are all uniquely determined by @xmath170 .
in this section we show that a two dimensional foam in equilibrium satisfies all conditions required for it to be a smvp . in order to do this
let us first write down the equilibrium conditions for an arbitrary vertex of the foam , in compact form .
we will consider the case of foams with arbitrary surface tensions , and also allow forces to act on the foam s vertices .
let @xmath185 be the location of a vertex @xmath58 at which three interfaces @xmath186 , @xmath187 and @xmath188 meet , as shown in fig .
[ fig : eq_sign_convention ] .
each interface @xmath189 is a circle arc with radius @xmath190 and center @xmath191 .
it produces on @xmath58 , due to its surface tension @xmath192 , a force @xmath193 of modulus @xmath194 in the direction of the tangent to the film at @xmath185 .
these forces @xmath195 can therefore be written as @xmath196 where @xmath197 and @xmath198 .
the asterisk stands here for a counterclockwise rotation in @xmath199 , so that @xmath200 .
+ let us more generally assume that an external force @xmath201 acts on @xmath58 .
equilibrium of all forces acting on the vertex implies @xmath202 there is still one more condition , which is related to pressure equilibrium around the vertex .
the pressure drop across an interface can be written as @xmath203 we will adopt the convention for @xmath204 to be positive if the pressure decreases when crossing the interface in the counterclockwise sense of rotation around @xmath58 .
this is of course related to a convention for the signs of the @xmath190 . in fig .
[ fig : eq_sign_convention ] @xmath205 and @xmath206 are negative according to this convention .
notice that pressure jumps and radii have different signs when considered from the two opposite vertices joined by a film .
+ the fact that the total accumulated pressure drop around a vertex has to be zero implies then @xmath207 this is the pressure equilibrium condition for the vertex .
equations ( [ forceeq ] ) and ( [ pressureeq ] ) are satisfied if the vertex is equilibrated , and are together equivalent to @xmath208 where @xmath209 is an _
arbitrary _ point .
this is what we will call the equilibrium condition for the vertex , and encloses both force and pressure equilibrium .
the alignment of the centers of a two - dimensional foam was already known long time ago by plateau  @xcite for equal surface tensions , and in the case of small self - standing clusters of 2 and 3 bubbles . for the case of a cluster of two bubbles and zero load
it is a trivial consequence of symmetry  @xcite .
the demonstration for the case of clusters of three bubbles is referred to by boys as being``rather long and difficult ''  @xcite . + it is not difficult to see that the alignment of the centers is in no way a property of clusters , and also not restricted to vertices with equal surface tensions and zero loads , but a general consequence of equilibrium .
we will find that under very general conditions such as arbitrary surface tensions and external loads applied on the vertices , if a vertex is equilibrated then the centers of the three arcs converging to it lay on a line .
+ consider equation  ( [ eq : aligned ] ) , and assume for a moment that the external load is zero .
then @xmath210 and @xmath211 lay on the same line , as can be seen by taking @xmath212 for example .
thus all vertices in equilibrium are aligned if no external force is applied on it .
this alignment property is even true under non - zero load conditions , if the force is perpendicular to the line of centers . a load satisfying such condition
will be called a _
normal load_. the alignment condition has the geometrical consequence that the interfaces @xmath189 , when continued , meet each other again at a unique point @xmath59 which we called the _ conjugate vertex _  @xcite of @xmath58 .
this also means that the vertex @xmath58 could be physically realized as a self - standing cluster of two bubbles ( by continuing its interfaces ) , and is a kind of a `` separability '' condition for the static equilibrium conditions , in the sense that each vertex of a foam might as well be that of an isolated cluster of two bubbles .
thus    _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ a vertex in equilibrium under a normal load is aligned .
_ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _    or , equivalently    _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ a vertex in equilibrium under a normal load has a conjugate vertex . _
_ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _    next we would like to consider the existence of a reciprocal figure , since this is also a condition that has to be satisfied in order for a foam to be a smvp . this condition must be separately considered .
the reader may easily build examples of bubbles with @xmath4 neighbors , all of whose @xmath4 vertices are aligned , but yet do not admit a reciprocal figure .
the reason is that the attempted reciprocal figure will not in general `` close '' around that bubble , the same case as described in appendix  [ sec : rprecognition ] .
we will now show that if all vertices of a foam are equilibrated , then an oriented reciprocal figure exists for it .
we start by considering a bubble and show that the reciprocal figure can be found for it . the figure for the whole foam
can then be formed by patching together those of neighboring bubbles .
consider a bubble with @xmath4 neighbors @xmath213 , as shown in figure [ fig : bubblea ] . at each vertex @xmath214 of this bubble ,
three films @xmath215 , @xmath216 and @xmath217 meet .
interface @xmath215 separates the central bubble from its neighbor @xmath218 , while interface @xmath217 is the limit between neighbors @xmath218 and @xmath219 . we will assume vertices @xmath214 to be in equilibrium under arbitrary _
loads @xmath220 .
these loads we can generally write as @xmath221 that is , we have decomposed the external load of vertex @xmath214 in the _ dependent _ basis formed by the three centers of the films meeting at the vertex . because these centers are aligned , and the load is normal to the line of centers , the following condition is always satisfied by the coefficients .
@xmath222 the equilibrium condition ( [ eq : aligned ] ) now reads @xmath223 and holds for _ arbitrary _
@xmath224 s .
the sign of @xmath225 is determined by the sign convention at vertex @xmath214 , therefore @xmath226 must appear with a minus sign in the equilibrium equation for vertex @xmath218 .
+ since we use a dependent basis ( equation .
( [ eq : loadgeneral ] ) ) the coefficients @xmath227 ( the `` representation '' of the load ) are not uniquely determined for a given load .
there is instead a monoparametric family of coefficients @xmath227 , all satisfying both ( [ eq : loadgeneral ] ) and ( [ eq : nullsum ] ) , for each load @xmath228 .
we will use these degrees of freedom to choose a representation @xmath229 that satisfies @xmath230 notice that this condition relates the coefficients of two consecutive loads .
we now show that such a representation always exists .
we start by making the degree of freedom in the representation of the loads explicit . for arbitrary @xmath231
, we add the null vector ( [ eq : geneq ] ) multiplied by @xmath231 to the load ( [ eq : loadgeneral ] ) and get @xmath232 where    @xmath233    condition ( [ eq : ring ] ) then implies @xmath234 which has always a solution in the generic case .
+ using this representation of the loads we can rewrite the equilibrium condition ( [ eq : geneq ] ) as @xmath235 where    @xmath236    this amounts to hiding the loads in a redefinition of the surface tensions and pressures : @xmath237 .
as we have shown , this can always be done for aligned loads .
we notice that the coefficients @xmath238 satisfy @xmath239 as can be verified by subtracting eq .
( [ eq : geneq1 ] ) with two different values of @xmath240 , and adding up the result for @xmath241 .
+ more generally the fact that all vertices of the bubble are in equilibrium has the consequence that @xmath242 for @xmath243 arbitrary .
condition ( [ eq : vectornullaround ] ) can be generalized to any closed path in the foam .
the sum is in that case over all films cut by the closed path .
the smallest such path , enclosing a vertex , gives the vertex equilibrium condition ( [ eq : aligned ] ) .
now that we have written the bubble equilibrium condition in the convenient form ( [ eq : geneq1 ] ) , we will write an algebraic condition for the existence of a reciprocal figure . + the reciprocal figure was defined to be a set of points @xmath244 , each associated to a bubble ( but not necessarily contained in it ) , such that the straight line passing through @xmath245 and @xmath246 is normal to the interface @xmath12 between bubbles @xmath247 and @xmath248 .
this means that @xmath249 , @xmath245 and @xmath246 are on the same line , as in fig .
[ fig : bubblec ] .
we will require that @xmath250 and @xmath251 be arbitrary ( with the only condition that @xmath250 , @xmath251 are @xmath252 are on a line ) .
this will allow us to patch together the reciprocal figures of neighboring bubbles to form that of the whole foam .
this is equivalent to the translation and dilatation degrees of freedom existent in the definition of a reciprocal figure for a svp .
for arbitrary @xmath250 , take @xmath251 anywhere on the line @xmath253 .
all other points are now uniquely determined .
@xmath254 is located in the intersection of @xmath255 and @xmath256 , next @xmath257 is found as the intersection of @xmath258 and @xmath259 , etc .
when the figure is closed with the last point @xmath260 , there is an extra condition , since it has to be the intersection of _ three _ lines : @xmath261 , @xmath262 , @xmath263 .
the construction of the reciprocal figure is thus an _ overdetermined _ problem , and would not in general be possible if the centers @xmath48 were arbitrarily located .
let us now write the conditions for this reciprocal figure to close , in an algebraic form .
the points @xmath0 of the reciprocal figure are determined by a set of coefficients @xmath264 satisfying the following conditions .
@xmath265    substituting [ rfcond1 ] into [ rfcond2 ] we get @xmath266 we will now show that the equilibrium conditions ( [ eq : geneq1 ] ) ensure that ( [ eq : rfcond3 ] ) always has a solution , and therefore that a reciprocal figure exists .
comparison with equation ( [ eq : geneq1 ] ) lets us conclude that a solution will exist for @xmath250 arbitrary , if we are able to find a set of coefficients @xmath264 that satisfy @xmath267 this is equivalent to    @xmath268    it is not difficult to see that these equations are satisfied if the @xmath269 are related by @xmath270 starting from an _ arbitrary _ @xmath271 , this recursion relation gives us the following values of @xmath272 .
once all are known , equation  ( [ eq : rel1 ] ) provides the values of the @xmath273 , which in turn determine @xmath274 using ( [ rfcond1 ] ) .
the condition that the figure can be closed is @xmath275 , and is equivalent to @xmath276 .
this condition is satisfied because ( [ eq : nullaround ] ) holds , and therefore is a consequence of equilibrium .
+ if we were given an arbitrary circular partition , it would not in general be possible to find a reciprocal figure for it .
the fact that this cp is an equilibrated foam imposes geometrical constraints on it ensuring , for example , that it admits a reciprocal figure .
+ we have thus shown that , for an arbitrary equilibrated bubble , it is always possible to find a reciprocal figure .
we can arbitrarily fix @xmath250 since @xmath224 in ( [ eq : geneq1 ] ) can be arbitrarily chosen , and we can also choose the `` scale '' @xmath277 of the drawing at will , since the starting value @xmath271 , that fixes this scale , is arbitrary .
therefore the reciprocal figures of neighboring bubbles can be patched together to form a reciprocal figure for the whole foam .
each selection of a starting point @xmath250 and a `` scale '' @xmath277 determines the other points uniquely , therefore there are three degrees of freedom in the determination of the reciprocal figure .
once @xmath250 and @xmath251 are chosen , all other points are found as intersections of two lines passing through the centers @xmath48 and one already existing point @xmath0 .
now we have to show that it is always possible to find , among all possible reciprocal figures , at least one that satisfies orientation as defined in section  [ sec : smvp ] .
in the first place , if all the surface tensions are positive then all films will be under traction and therefore the vertices will be convex .
if there are non - convex vertices in the foam ( which would happen if some of the films are compressed instead of stretched ) then we know that it is not possible to satisfy orientation ( fig .
[ fig : convexity ] ) .
positive surface tensions is thus a necessary condition  @xcite for the foam to be a smvp , although their modulus can be arbitrary for each film . + the orientation condition could fail in the first place because a center @xmath48 lays in between two points @xmath0 and @xmath35 .
it is always possible to avoid this by choosing @xmath251 close enough to @xmath250 . in this way
all following points @xmath0 are confined within a ( arbitrarily chosen ) small region of space where no center is located .
this ensures that no center @xmath48 lays between @xmath0 and @xmath35 . + now regarding the second part of the orientation condition ( sec .
[ sec : smvp ] ) , consider a vertex @xmath130 , which is convex and aligned .
two orientations of the triangle @xmath278 are possible , as shown in figures [ fig : convexity]a and [ fig : convexity]c .
notice that both constitute reciprocal figures for both vertices , but only one of them satisfies orientation in each case . in our construction of the reciprocal figure for the whole foam
, we can decide the orientation of the initial triangle , choosing the one that respects orientation .
the question is now if the correct orientation of this starting triangle ensures that of all subsequent ones , whose locations are determined by @xmath250 and @xmath251 . to demonstrate that this is indeed the case
, we notice that the triangle of figure [ fig : convexity]a , if considered as a reciprocal figure for the vertex of [ fig : convexity]c , has all three segments @xmath1 wrongly oriented .
the point we want to make is that , if the vertex is convex , there are only two possibilities : either all pairs @xmath279 satisfy orientation , or all are wrongly oriented .
then if one of the pairs of a triangle forming part of a reciprocal figure is correctly oriented , the other two must necessarily also be .
this demonstrates that if the starting pair @xmath280 is chosen with the correct orientation , then all subsequent triangles must be correctly oriented since they share at least a pair of sources with one preexisting triangle . therefore in order to ensure correct orientation of the whole figure
it is enough to correctly choose the orientation of the first pair .
a dissection of space into cells separated by circular interfaces that meet at triple vertices is called a circular partition ( cp ) .
a two - dimensional foam therefore defines a circular partition of two - dimensional space .
the equilibrium conditions for the foam impose geometrical constraints on this cp .
we have here shown that the cp defined by an equilibrated foam is _ aligned _ and admits an _ oriented reciprocal figure_. this result is valid in general for heterogeneous foams , each of whose films may have an arbitrary ( positive ) surface tension , and even if loads are applied on the vertices , with the sole requirement of equilibrium .
we have seen in sec .
[ sec : find_sources ] , that any cp satisfying the conditions of alignement and existence of oriented reciprocal figure is a sectional multiplicative partition ( smvp ) .
a smvp is a plane cut of a multiplicative partition , thus two - dimensional foams are plane cuts of three - dimensional foams , these last being a multiplicative partition with respect to sources in three dimensions . therefore given an arbitrary equilibrated two - dimensional foam , it is always possible to find sources @xmath3 in three - dimensional space , and amplitudes @xmath41 such that the given foam is a smvp with respect to those sources .
+ a first implication of this correspondence is the identification of a new set of degrees of freedom for the foam : the intensities @xmath41 and locations of the sources @xmath8 in three dimensions .
this allows a more natural description of a foam , than the one that is done in terms of films and vertices . when a foam is interpreted as a tessellation of space with respect to some sources , we see the foam s films and vertices are secondary constructions , and their evolution is a consequence of that of the sources .
the dynamical description is conceptually simpler using the smvp interpretation .
for example the processes of neighbor switching ( @xmath281 ) and cell disappearance ( @xmath282 ) are described in a unified manner ( fig .
[ fig : process ] ) .
both are due to the fact that a fourfold vertex in 3d crosses the projection plane @xmath28 .
depending on the spatial orientation of the vertex with respect to @xmath28 , this is seen as a @xmath281 or @xmath282 process . + an evolving foam can be seen as a particular instance of a dynamical random lattice@xcite , in which the evolution of a cellular structure is fixed by assigning a given dynamics to the sources of a mathematically defined tessellation . in the case of foams
the dynamics is usually fixed by gas diffusion across the membranes .
alternatively other dynamical evolution rules may also be interesting , but in any case one has to make the translation to obtain the source s dynamics .
the next step is then to find , for a given proposed dynamic evolution for the foam , the corresponding dynamics@xcite for the sources and intensities .
+ foams are usually studied inside a bounded region or cage , which imposes the only constraint that the films be normal to its boundary .
boundaries of this kind do not affect the fact that the foam is a smvp .
we have shown that the result holds with the sole requirement of vertex equilibrium .
but it is obvious that smvps are always closed on themselves forming a self standing cluster , that is , there are no open film ends .
this implies that even bounded foams _ must _ be a region of a larger self - standing cluster of bubbles that is closed on itself , and everywhere equilibrated .
the point is not trivial in that it ensures that all films ending at the boundary can be continued , eventually forming new ( phantom ) vertices , and that the resulting foam will have all its vertices in equilibrium .
we see then that there is no fundamental difference between bounded and self - standing foams , since all foams are regions of a self - standing cluster .
this does not mean that the boundaries have no effect , which would of course be false .
if the foam s dynamic inside the cage is for example produced by gas diffusion across the films , the evolution of the `` phantom '' bubbles existing as continuations of the physical foam outside the cage , will _ not _ follow this dynamic , but a different one , which is determined by the constraint that the films be normal to the cage s boundaries .
+ in the field of joint - bar structures , an old result due to maxwell@xcite states that if a lattice accepts a reciprocal figure then it can support a self - stress , and conversely .
more recently ash and bolker@xcite have shown that the existence of a reciprocal figure is sufficient condition for the lattice to be a sectional partition . in this case
there is the additional requirement that all vertices are convex , therefore all stresses in the lattice must be of the same sign , and the lattice can be an equilibrated spider web .
a chain of results that span a century allow us to see equilibrated spider webs as svps , and conversely .
the alert reader may have noticed that equilibrated foams can be seen as a kind of `` generalized '' spider webs , in which the pressure difference between cells is the new ingredient , and is equilibrated by the curvature of the interfaces .
it therefore turns out as no surprise that these generalized spider webs ( foams ) are equivalent to an appropriate generalization of svps , namely smvps , which include a multiplicative constant that gives rise to curved interfaces .
+ i have benefited from discussions with k.  lauritsen , h.  herrmann , h.  flyvbjerg , d.  le caer and d.  weaire .
i thank d.  weaire for reference @xcite , and d.  le caer for making me aware of plateau s results and sending me references @xcite and @xcite .
i am also grateful to h.  flyvbjerg and d.  stauffer for their efforts to help me make this work more readable . +
a classification of space into cells @xmath5 separated by rectilinear interfaces @xmath12 is called a _ rectilinear partition _ of space . given a rectilinear partition @xmath73 of the plane , a _ reciprocal figure _ @xmath283is a planar graph composed of sites @xmath247 joined by edges @xmath71 satisfying  @xcite :    * sites @xmath247 of @xmath70 are in one - to - one correspondence with cells @xmath5 of @xmath73 .
* edges @xmath71 of @xmath70 are in one - to - one correspondence with interfaces @xmath12 of @xmath73 .
* edges @xmath71 of @xmath70 are * normal * to interfaces @xmath12 of @xmath73 .    clearly a reciprocal figure is defined up to arbitrary global dilatations and translations , since angles are not changed by them .
therefore , if @xmath73 admits a reciprocal figure , there are three degrees of freedom in its determination  @xcite .
an arbitrary partition @xmath73 will _ not _ in general admit a reciprocal figure .
we can see this with a simple example .
draw an arbitrary polygonal cell @xmath284 with @xmath4 faces , and take arbitrary rectilinear interfaces between its @xmath4 neighbors(fig .
[ fig : closing ] ) .
now take an arbitrary point @xmath285 on the plane to start with , an assign it to the central cell .
this starting point is arbitrary since a reciprocal figure is defined up to arbitrary translations .
the other @xmath4 points @xmath286 associated to the external cells must be somewhere on the @xmath4 rays @xmath190 stemming from @xmath285 and normal to the faces @xmath287 of @xmath284 .
the global length scale is also arbitrary so that say @xmath288 can be freely chosen .
then we choose @xmath289 somewhere on @xmath205 .
now point @xmath290 is determined as the intersection of ray @xmath291 with the normal to face @xmath94 going through @xmath289 .
this can be repeated to obtain all @xmath4 external points , but in general the figure will not _ close _ , that is , segment @xmath292 will not be normal to interface @xmath293 , whose orientation is arbitrary .
+ the planar graph formed by joining the sources @xmath0 of a svp with edges @xmath71 , one for each nonempty interface @xmath12 , constitutes a reciprocal figure for the svp .
therefore every svp has a reciprocal figure .
@xmath294 recently ash and bolker have shown that the existence of a reciprocal figure satisfying _ orientation _
@xcite is also a _ sufficient _
condition for @xmath73 to be a svp .
a reciprocal figure @xmath70 satisfies orientation if for each bond @xmath295 , the sites @xmath247 and @xmath248 are oriented in the same way as cells @xmath5 and @xmath11 are with respect to @xmath12 .
@xmath296 this result solves the recognition problem for svps .
the orientation condition can only be satisfied if all vertices of @xmath73 are _
convex_. we will say that a vertex is convex if the internal angles formed by the interfaces are all smaller than @xmath79 .
figure  [ fig : orientation ] shows two partitions with three cells .
one of them is convex , the other is not . in the second case
it is not possible to find a reciprocal figure that satisfies orientation , and therefore it can not be a svp . + given a rectilinear partition @xmath297 , and an oriented reciprocal figure for it ( as in fig .
[ fig : recfig ] ) , it is always possible@xcite to find sources @xmath3 in three dimensions , located at heights @xmath30 above the points @xmath0 , such that @xmath297 is the section with @xmath28 of a three - dimensional vp with respect to @xmath3 . the procedure to determine the heights @xmath30 can be easily described .
vertices @xmath19 of @xmath297 will be equidistant from sources @xmath8 , @xmath13 , @xmath20 ( see fig .  [
fig : svplocate_sources ] ) . start from an arbitrary vertex , say @xmath19 , and draw a spherical surface of _ arbitrary _ radius @xmath298 with center at that vertex .
now define sources @xmath8 , @xmath13 and @xmath20 as the intersections of this surface with the verticals ( normals to @xmath28 ) through @xmath0 , @xmath35 and @xmath299 respectively .
next go to a neighboring vertex , which shares two sources with this one .
let us call it @xmath300 . for this vertex
, only source @xmath301 has to be located since @xmath8 and @xmath13 are known . draw a spherical surface with center @xmath300 and containing @xmath8 and @xmath13 .
both will be simultaneously contained , since interface @xmath12 , on which @xmath300 is located , is equidistant from @xmath8 and @xmath13 .
the radius @xmath302 of this surface is determined by the locations of @xmath8 and @xmath13 , which in turn is fixed by @xmath303 .
its intersection with the vertical @xmath304 through @xmath305 determines @xmath301 .
if this spherical surface does not intersect @xmath304 , just choose a larger value of @xmath303 and start all over again ( from the initial vertex ) .
the construction proceeds in this manner until all sources have been determined .
as mentioned , the initial value of @xmath303 is tentative , in the sense that it may have to be modified ( increased ) if at some point during the procedure , a normal line is not cut by the corresponding spherical surface from the vertex .
it is easy to see that increasing the value of the starting radius @xmath303 is always enough to solve this problem .
+ there is thus one degree of freedom in this construction ( @xmath303 ) .
we conclude that , given a two - dimensional partition @xmath73 that admits a reciprocal figure , there is a four - parametric family of source locations such that @xmath73 is a svp with respect to them .
three degrees of freedom come from the determination of the reciprocal figure itself ( since a dilatation and/or translation of a reciprocal figure is again a reciprocal figure ) and the last one from @xmath303 .
this last degree of freedom results from the fact that a svp is invariant if all heights are changed according to @xmath306 with @xmath218 arbitrary ( see equations ( [ eq : svpcell ] ) and ( [ eq : svpinterface ] ) ) .
+ reciprocal figures were first studied by maxwell  @xcite in relation with the rigidity of bar - joint frameworks in the plane .
maxwell pointed out that frameworks that have a reciprocal figure are able to support a self - stress , and conversely .
the reason is that the edges of the reciprocal figure can be taken to represent forces transmitted by the edges of the original framework ( rotated by @xmath307 ) .
since these edges form closed polygons , the existence of a reciprocal figure implies the existence of an equilibrated set of internal stresses in the absence of external load .
the addition of the orientation condition ( a condition not required by maxwell s definition of a reciprocal figure ) has the statical consequence that all signs of the stresses are equal , for example all traction or all compression .
it is clear that no equilibrium is possible in the case of fig .
[ fig : orientation]a if all three stresses are to have the same sign . fig .
[ fig : orientation]b on the other hand , can be in equilibrium under compression or traction on all three interfaces .
+ the conclusion is that every svp is an equilibrium configuration for a spider web  @xcite , and conversely , each such equilibrium configuration is a svp . +
the existence of a reciprocal figure has also projective consequences , which have been studied by crapo  @xcite and whiteley  @xcite .
we briefly describe here a geometric transformation called _ inversion _  @xcite
. we will find it extremely useful for our purpose of discussing circular partitions .
an inversion with radius @xmath21 around a point @xmath308 located at @xmath309 transforms a point @xmath310 located at @xmath311 into a point @xmath312 at @xmath313 satisfying @xmath314 where @xmath315 .
the sphere of radius @xmath21 and centered at @xmath308 is invariant under this transformation , while the inside and outside of this sphere are interchanged .
obviously this transformation is self - inverse .
let us now describe some of the properties of this transformation in two dimensions  @xcite .
most of them apply trivially in higher dimensions .
* circles not through @xmath308 are transformed in circles not trough @xmath308 . *
circles through @xmath308 are transformed into straight lines not through @xmath308 . *
straight lines not through @xmath308 are transformed in circles through @xmath308 .
* straight lines through @xmath308 are invariant .
* angles are preserved ( in modulus ) by the inversion .    given two points @xmath251 and @xmath254 at distances @xmath205 and @xmath291 from the inversion center , the distance @xmath316 between them transforms as @xmath317 using this result we can easily see that the inversion is a `` symmetry '' of a mvp in any dimension if the intensities are also appropriately transformed  @xcite .
more precisely , a mvp of @xmath318 with respect to sources @xmath22 with intensities @xmath50 is transformed by an inversion into a mvp with respect to @xmath319 with intensities @xmath320 , where the new intensities satisfy @xmath321 here @xmath322 is an arbitrary prefactor , the same for all @xmath41 s , and @xmath190 is the distance between source @xmath247 and the inversion center @xmath308 . to see this
it suffices to demonstrate that if @xmath323 then after an inversion , @xmath324 , which is easily done using ( [ eq : mvpinterface ] ) , ( [ eq : dist_transf ] ) and ( [ eq : transforma ] ) .
the inversion transformation is of course also a symmetry of the smvp ( section  [ sec : smvp ] ) , if the inversion center @xmath308 is on the cutting plane @xmath28 , since in this case the inversion leaves this plane invariant .
+ if the inversion center @xmath308 happens to be located on an interface @xmath12 of a mvp ( initially a spherical surface ) , the transformed interface @xmath101 is a plane not through @xmath308 .
the resulting interface thus corresponds to a partition with respect to sources @xmath247 and @xmath248 in their new locations .
therefore the transformed intensities @xmath325 and @xmath326 have to be equal after the inversion , which is verified using ( [ eq : transforma ] ) @xmath327 in the same way a smvp with respect to two sources @xmath247 and @xmath248 is transformed into a svp if @xmath328 .
the intensities @xmath41 are transformed according to ( [ eq : transforma ] ) , where @xmath190 is the distance between @xmath308 and source @xmath8 in three - dimensional space .
the way in which heights @xmath30 transform is easily found using ( [ eq : inversion ] ) . notice that if the inversion center coincides with the location of a conjugate vertex @xmath59 , then the transformed partition has a rectilinear vertex since three interfaces are simultaneously transformed into straight lines .
we will use this property of the inversion transformation when solving the recognition problem for smvps in section  [ sec : find_sources ] .
n.  rivier , _ on the structure of random tissues or froths , and their evolution .
_ , philosophical magazine b * 47 * ( 1983 ) , l45-l49 .
d.  weaire and n.  rivier , _ soap , cells and statistics
random patterns in two dimensions _ , contemp .
* 25 * ( 1984 ) , 59 - 99 .
j.  glazier and d.  weaire , _ the kinetics of cellular patterns _ , j.phys.:condens .
matter * 4 * , ( 1992 ) , 1867 - 1894 .
j.  stavans , _ the evolution of cellular patterns _ , rep
.  prog .  phys .
* 56 * ( 1993 ) , 733 - 789 . c.  s.  smith ,
_ the shape of things _ , sci .  am .
* 190 * ( 1954 ) , 58 . c.  s.  smith , _ grain shapes and other metallurgical applications of topology _ , in : metal interfaces , am .
metals , cleveland , 1952 .
d.  weaire and j.  kermode , _ computer simulation of a two - dimensional soap froth i. method and motivation _ , philosophical magazine b
* 48 * ( 1983 ) , 245 - 259 .
+ j.  kermode and d.  weaire , _
2d - froth : a program for the investigation of 2-dimensional froths _ , comput
.  phys .
* 60 * ( 1990 ) , 75 - 109 .
+ t.  herdtle and h.  aref , _ numerical experiments on a two - dimensional foam _ , j.  fluid mech .
* 241 * ( 1992 ) , 233 - 260 .
h.  flyvbjerg and c.  jeppesen , _ a solvable model for coarsening soap froths and other domain boundary networks in two dimensions _ , physica scripta
* t38 * , ( 1991 ) , 49 - 54 . + h.  flyvbjerg , _ dynamics of soap froth _ , physica a * 194*,(1993 ) , 298 . + h.  flyvbjerg , _ model for coarsening froths and foams _ , phys .  rev .
e * 47 * , ( 1993 ) , 4037 .
k.  kawazaki , t.  nagai and k.  nakashima , phil .  mag .
b * 60 * , ( 1989 ) , 399 .
+ k.  kawasaki , physica a * 163 * , ( 1990 ) , 59 . , a.  okabe , b.  boots and k.  sugihara , wiley 1992
. h.  honda , _ description of cellular patterns by dirichlet domains : the two - dimensional case _
, j.  theor .
* 72 * , ( 1978 ) , 523 . v.
icke and r. van de weygaert , fragmenting the universe .
i. statistics of two - dimensional voronoi foams , astron
. astrophys . *
184 * , ( 1987 ) , 16 - 32 .
hiroshi imai , masao iri and kazuo murota , _ diagram in the laguerre geometry and its applications .
_ , siam j.  comput . * 14 * , ( 1985 ) , 93 . b.  n.  boots , _ weighting thiessen polygons _
, _ economic geography _ ( 1979 ) , 248 - 259 . , c.  vernon boys , ny : dover 1959 , pp . 120 - 127 .
+ reprinted in : _ the world of mathematics _ , tempus books of microsoft press 1988 vol .
ii , pp . 883 - 886 . , darcy wentworth thompson , cambridge university press 1961 , page 96 .
maxwell , j.  c. _ on reciprocal figures and diagrams of forces _ , _ phil .  mag .
series _ * 4 * ( 1864),250 - 261 .
maxwell , j.  c. _ on reciprocal figures , frames , and diagrams of forces _ , _ trans .
royal soc .
edinburgh _ * 26 * ( 1869 - 72),1 - 40 . peter f.  ash and ethan d.  bolker , _ generalized dirichlet tessellations _
, _ geometriae dedicata _ * 20 * ( 1986),209 - 243
. henry crapo , _ structural rigidity _
, _ structural topology _ * 1 * ( 1979 ) , 26 - 45 . walter whiteley , _ realizability of polyhedra _ , structural topology , * 1 * ( 1979 ) , 46 - 58 . h.  telley , _ modelisation et simulation bidimensionelle de la croissance des polycrystaux _ , unpublished , thesis report n. 780 .
ecole polytechnique federale de lausanne , ( 1989 )
. this is called `` phantom vertex '' in @xcite .
, h.  s.  m.  coxeter , second edition , wiley and sons 1989 .
a configuration @xmath329 is formed by @xmath330 points and @xmath331 lines , such that each line is adjacent to @xmath79 points , and each point is adjacent to @xmath332 lines .
see for example , _ geometry and the imagination _ , d.  hilbert and s.  cohn - vossen , chelsea publishing co. , new york 1952 , chapter iii . c.  moukarzel,_voronoi foams _
, physica * a 199 * ( 1993),19 - 30 . k.  lauritsen , c.  moukarzel and h.  j.  herrmann , _ statistical physics and mechanics of random lattices _ , _ j.  physique  i _ ( france ) * 3 * ( 1993),1941 . c.  moukarzel , to be published
. within our generalized notion of distance , that is , including the multiplicative constants .
although it is statically unstable , a foam all of whose films are under compression would also be geometrically possible .