Interface growth in the channel geometry and tripolar Loewner evolutions
A class of Laplacian growth models in the channel geometry is studied using the formalism of tripolar Loewner evolutions, in which three points, namely, the channel corners and infinity, are kept fixed. Initially, the problem of fingered growth, where growth takes place only at the tips of slit-like fingers, is revisited and a class of exact exact solutions of the corresponding Loewner equation is presented for the case of stationary driving functions. A model for interface growth is then formulated in terms of a generalized tripolar Loewner equation and several examples are presented, including interfaces with multiple tips as well as multiple growing interfaces. The model exhibits interesting dynamical features, such as tip and finger competition.
pacs:68.70.+w, 05.65.+b, 61.43.Hv, 47.54.–r
The Loewner equation describes a rather general class of growth processes in two dimensions where a curve starts from a given point on the boundary of a domain in the complex -plane and grows into the interior of . More specifically, the Loewner equation loewner () is a first-order differential equation for the conformal mapping from the ‘physical domain,’ consisting of the region minus the curve, onto a ‘mathematical domain’ represented by itself. The specific form of the Loewner equation depends both on the domain and on the point where the trace ends. The most studied cases are the chordal and radial Loewner equations review3 (), where in the former case the trace ends on a given point on the boundary, while in the latter it ends on a point in the interior. Other geometries such as the dipolar case BB2005 () and the channel geometry poloneses () have also been considered in the literature. The Loewner equation depends explicitly on a driving function, here denoted by , which is the image of the growing tip under the mapping . An important development in the field was the discovery by Schramm schramm () that when is a Brownian motion the resulting Loewner evolution describes the scaling limit of certain statistical mechanics models. This result spurred great interest in the so-called stochastic Loewner equation (SLE) and the subject has now been widely reviewed, see, e.g. review1 (); review2 (); cardy (); BB ().
Although the SLE has attracted most of the attention lately, the deterministic version of the Loewner equation remains of considerable interest both in a purely mathematical context kadanoff_jsp () and in connection with applications to growth processes selander (); makarov (); poloneses () and integrable systems theory zabrodin2 (). Indeed, the deterministic Loewner equation has been used to study Laplacian fingered-growth in both the chordal and radial cases selander (); makarov () as well as in the channel geometry poloneses (). In this class of models, growth takes place only at the tips of slit-like fingers and the driving function has to follow a specific time evolution in order to ensure that the tip grows along gradient lines of the corresponding Laplacian field. Although this thin-finger model poloneses () was able to reproduce some of the qualitative behavior seen in experiments combustion (), treating the fingers as infinitesimally thin is a rather severe approximation.
More recently, the growth of “extended fingers” us_pre () or “fat slits” zabrodin1 (), meaning a domain encircled by an interface with endpoints on the real axis, was considered within the formalism of deterministic Loewner evolutions in the upper half-plane. In particular, the interface growth model considered in Ref. us_pre () was shown to exhibit interesting dynamical features akin to the finger competition mechanism observed in actual Laplacian growth pelce (). Many growth processes, however, take place in a bounded domain, e.g., within a channel pelce (), for which a Loewner-equation approach was still lacking. To develop a formalism based on Loewner evolutions to describe the growth of an interface in the channel geometry is thus the main goal of the present study. Because the Loewner function fixes the points (the channel corners) and (the channel ‘open’ end), we shall refer to growth processes in the channel geometry as tripolar Loewner evolutions, to distinguish it from the dipolar case where two points are kept fixed BB2005 ().
We begin our analysis by revisiting in Sec. II the problem of fingered growth in the channel geometry, where a curve (finger) grows from a point on the real axis into a semi-infinite channel. We present a brief derivation of the corresponding tripolar Loewner and report a new class of exact solutions for the case of multiple fingers with stationary driving functions. In Sec. III we then discuss a general class of models where a domain, bounded by an interface, grows into a channel starting from a segment on the real axis. In our model, the growth rate is specified at a certain number of special points on the interface, referred to as tips and troughs, which then determine the growth rate at the other points of the interface according to a specific growth rule (formulated in terms of a polygonal curve in the mathematical plane). A generalized tripolar Loewner equation for this problem is derived and several examples are given in which the interface evolution is obtained by a direct numerical integration of the Loewner equation. Our main results and conclusions are summarized in Sec. IV.
Ii Loewner Evolution for Fingers in the Channel
ii.1 Tripolar Loewner Equation
In order to set the stage for the remainder of the paper and to establish the relevant notation, we begin our discussion by considering the simplest Loewner evolution in the channel geometry, namely, that in which a single curve grows from a point on the real axis into a semi-infinite channel whose side walls are placed at and the bottom wall is at ; see Fig. 1. Here the relevant domain, , for the Loewner evolution is
The curve at time is denoted by and its growing tip is labeled by . Now let be the conformal mapping that maps the ‘physical domain,’ corresponding to the upper half-channel in the -plane minus the curve , onto the upper half-channel in an auxiliary complex -plane, referred to as the ‘mathematical plane,’ i.e., we have
with the curve tip being mapped to a point on the real axis in the -plane; see Fig. 1.
The mapping function is required to satisfy the so-called ‘hydrodynamic condition’ at infinity, namely,
where is a real-valued, monotonously increasing function of time satisfying the condition . Condition (1), which implies , together with the requirement , uniquely determines (up to a reparametrization of the time ) the mapping . The time parametrization is specified by fixing the function , which in turn is related to the growing rate of the curve; see below. The Loewner function also satisfies the initial condition
since we start with an empty channel.
We consider the growth process to be such that growth occurs only the tip of the curve and in such a way that the accrued portion from time to time , where is an infinitesimal time interval, is mapped under to a vertical slit of (infinitesimal) height in the mathematical -plane; see Fig. 1. For convenience of notation, we shall represent the mathematical plane at time as the complex -plane, so that the mapping maps the physical domain at time onto the upper half-channel in the -plane, with the new tip being mapped to the point .
In order to derive the Loewner equation for the problem formulated above it is necessary to write in terms of , i.e., , where is the conformal mapping between the mathematical - and - planes, and then consider the limit . Alternatively, one can transform the growth problem defined in the channel geometry into a corresponding Loewner evolution in the upper half-plane. To do that, we apply the following transformations:
so that the channels in the -, -, and - planes are mapped to the upper half-plane in the auxiliary complex -, -, and - planes, respectively. Let us also introduce introduce the following notation:
Next, consider the mapping from the upper half--plane onto the upper half--plane with a vertical slit, so that one can write
The slit mapping can be easily computed (see Ref. poloneses () for details), and after taking the limit , one obtains the following Loewner equation
and the growth factor is defined by
From the boundary condition (1) it follows that . Without loss of generality we shall take in this section, implying that . Note that the Loewner equation (11) indeed fixes the points as well as the point at infinity, in the sense of Eq. (1). Hence we shall refer to Eq. (11) as the tripolar Loewner equation, in analogy with the dipolar case BB2005 () where only the points are kept fixed. The tripolar Loewner equation in the upper half-plane given in Eq. (7) also fixes the points and the point at infinity. More specifically, the boundary condition at infinity reads: , for .
The Loewner equation given in Eq. (7) can be readily extended to the case of multiple curves growing simultaneously in the channel, yielding
We shall assume for simplicity that the growth factors are all constant, i.e., , subjected to the condition , as follows from our choice above for . Other growth models have been considered poloneses () where the normal velocity is of the form , with the growth factors being given by , where is the inverse of . In this case, the analysis of the Loewner equation is much more difficult because of the implicit dependence of the growth factors on the Loewner function . By considering the simpler (but still very interesting) case where the growth factors are constant in time, we are able both to find exact solutions for multifingers and to easily integrate the Loewner equation (13) on the computer.
ii.2 Exact solutions
An exact solution for the Loewner function for a single curve in the channel geometry can be found explicitly for a stationary driving function, i.e., , corresponding to a finger growing vertically along the channel centerline poloneses (). Here we extend such stationary solutions to the case of multifingers. For simplicity we consider only the case of two fingers, but the procedure outlined below applies, in principle, to any number of fingers (although the calculations become increasingly more difficult).
For two fingers, Eq. (14) becomes
whose stationary solution (fixed point), and , are given by
Note that if , one has and , and the solution corresponds to two identical fingers growing vertically. This symmetrical two-finger solution is, of course, related to the single-finger solution described in Ref. poloneses () by a simple transformation, since in this case we can view each finger as growing along the centerline of a channel with half the width of the original channel. Thus, we are mainly interested here in asymmetrical solutions where . In the remainder of this section we shall omit, for convenience of notation, the upperscripts from and .
The Loewner equation (13) for two fingers reads
where we have used the fact that . This equation can be easily integrated, yielding an implicit solution for the mapping in the form
The tip trajectories , for , can now be obtained by setting and in Eq. (19), where . One then obtains
where the roots for and must be chosen in the upper half-plane. The asymptotic trajectories in the limit can be found by writing and considering the limit . After a straightforward calculation one obtains
We thus see that for large times the two fingers approach the vertical axes
respectively, and both of them move asymptotically with the constant speed .
An example of such a solution is shown in Fig. 2(a), where we used and , which from (17) implies and . In this case, the fingers should approach the vertical axes and , as predicted by (23), and this is indeed observed in the figure. Note that this asymptotic behavior holds for any initial conditions and , on account of the fact that the fixed point given in (17) is stable. This is illustrated in Fig. 2(b) where we show a solution with initial conditions and . To generate the curves shown in Fig. 2 we integrated the Loewner equation given in Eq. (13) using the numerical scheme described in us_pre ().
ii.3 Tripolar Loewner Chains
In the spirit of chordal Loewner chains BB (), it is clear that the tripolar Loewner equation (13) can be extended to describe the growth of more general domains in the channel geometry, not restricted to a curve or a set of curves. The corresponding tripolar Loewner chain thus reads
where the density of singularities can be viewed as a measure of the growth rate at a point at the boundary of the growing set that is the preimage of under . Since the shape of the growing domain is fully encoded in the map , equation (25) specifies the growth model once the density is known. As discussed elsewhere us_pre (), the formalism of Loewner chains is of limited practical use, except when the density is a sum of Dirac -functions, in which case we are back to multiple growing curves. In this context, it is important to consider more explicit models to describe the growth of a domain bounded by an interface. One such model was first introduced in Ref. us_pre () for the upper half-plane, and in the next section it is extended to the channel geometry.
Iii Loewner Evolutions for Interfaces in the Channel
iii.1 Tripolar Loewner Equation for Interfaces
Here we consider the problem of an interface starting from a segment on the real axis in the -plane and growing into the upper half-channel, as indicated in Fig. 3. In our growth model we assume, for simplicity, that the growing interface has certain special points, referred to as tips and troughs, where the growth rate is a local maximum and a local minimum, respectively, while the interface endpoints and remain fixed. Let us denote by the interface at time and by the growing region delimited by and the real axis. It is assumed that the curve is simple so that the physical domain is simply connected, hence is a hull. As before, the Loewner function maps the physical domain in the -plane onto the upper half-channel in the mathematical plane -plane, that is,
Under the action of , the interface is mapped to an interval on the real axis in the -plane, with the images of the tips, troughs and end points being denoted by , , where is the total number of such special points; see Fig. 3. The mapping function is again required to satisfy satisfy the hydrodynamic condition (1) and the initial condition (2). Here, however, we shall not fix the function a priori.
The growth dynamics in our model is specified by requiring that tips and troughs grow along gradient lines in such a way that the interface at time , for infinitesimal , is mapped under to a polygonal curve, as shown in Fig. 3. The exterior angle at the -th vertex of our polygonal curve is denoted by , with the convention that if the angle is greater than the corresponding parameter is negative. In other words, the parameter ’s are negative for tips and positive for troughs and end points. From a trivial trigonometric relation it follows that
If we denote by the heigth of the th vertex of the polygonal curve in the -plane, then as the parameters and all go to zero, and in this limit one finds the following relation:
As before, the mathematical plane at time is represented by the complex -plane, so that under the mapping , the interface tips, troughs, and end points at time are mapped to the points , as indicated in Fig. 3. To derive the tripolar Loewner equation for the growing interface we first need to obtain the mapping , so that we can write , and then take the limit . As discussed in Sec. II, it is more convenient however to work with the tripolar Loewner equation in the upper half-plane by applying the transformations given in Eq. (3). The corresponding transformed domains are shown in Fig. 4. In the -plane, we obtain (in the limit of ) a polygonal curve whose vertices are located at the points , where and the heigth in the -plane is given by
Next consider the mapping, , from the upper half--plane onto the upper half--plane with a polygonal cutout region; see Fig. 4. Since the domain in the -plane can be viewed as a degenerate polygon, the mapping function can be easily obtained from the Schwarz-Christoffel transformation CKP (), which yields
The constants and are determined so as to guarantee that , and hence we can rewrite Eq. (29) as
The integral in Eq. (30) cannot, in general, be performed explicitly. It is thus convenient to expand the integrand in powers of the infinitesimal quantities and then proceed with the integration term by term. Doing this up to the first order in the ’s, one finds
where the ‘tilded’ growth factors are defined by
The time evolution of the driving functions is determined by requiring that the points be taken under the mapping to the points . In other words, we impose that
Applying the procedure used above to derive the Loewner equation Eq. (33), it immediately follows that the differential equations governing the dynamics of the driving functions can be obtained by simply setting in Eq. (33), which yields
The tilded growth factors defined in Eq. (36) can be related to the growth factors, , defined for the original channel geometry by virtue of trivial trigonometric relations. After some calculation, one finds
In view of the relations above, we can express the parameters and associated with the end points in terms of the other growing factors , so that our growth model is completely specified by prescribing the functions at the tips and troughs of the interface. Thus, once the growth factors are given, one can compute the tilded growth factors from Eq. (39), proceed with the integration of the Loewner equation given in Eq. (33), and then invert the transformation (5) to obtain the Loewner mapping . Next we discuss some examples of interface growth described by the above model.
For an interface with a single tip (i.e., ), the factors and can be written in terms of , so that we can set without loss of generality, since this amount to a mere rescaling of the time coordinate. Let us consider first the case in which the interface is symmetric with respect to the channel centerline: and . This imply, in particular, that . Introducing the parameter , the Loewner equation (33) then becomes
with the dynamics for being given by
In this case, one can show from Eq. (39) that becomes
where we have used the fact that .
An example of a symmetric interface in a channel is shown in Fig. 5. In this figure, the solid curves represent the interface at various times , starting from up to with a time separation of between successive curves. To generate the curves shown in Fig. 5 we integrated the Loewner equation (43) backwards in time from ‘terminal conditions’ , for , to obtain the respective initial values from which we determine the corresponding points on the interface: . One sees from the figure that as time proceeds the interface expands and tends to occupy the entire channel. It is worth mentioning, however, that it is hard to go past the latest time shown in Fig. 5 because the function becomes very close to unity, rendering the numerical integration of the Loewner equation very difficult after this point.
In the case of an asymmetric interface, we need to integrate the full Loewner equation given in Eq. (33) with . In Fig. 6 we show a solution where the solid curves represent the interface at various times (see caption for details), while the dashed curve indicates the path traced by the tip . As seen in Fig. 6, the tip is “repelled” by its image with respect to the closest channel wall, so that an asymptotically symmetrical shape is expected for sufficiently long time. (It is difficult, however, to integrate the Loewner equation for larger times for reasons already mentioned.)
Let us now illustrate some of the patterns that emerge from our growth model for interfaces with two tips and one trough (i.e., ), in which case tip competition can arise. Here we shall consider for simplicity only cases in which We expect nevertheless that the qualitative behavior seen in this simpler situation (see below) should also be valid in more general cases. As a first example, we show in Fig. 7 a symmetrical growing interface, where the dashed lines represent the trajectories of the tips and the trough. In this figure, we started with symmetrical initial conditions, namely, , , and , and chose the growth factors of the two tips to be the same, , so as to preserve the initial symmetry, with the trough growth factor being .
If we break the symmetry, either by choosing an asymmetric initial condition or using different growth factors for the tips, then inevitably of the tips will grow faster and “screen” the other tip. For instance, in Fig. 8 we show a situation where the growth factors are the same as in Fig. 7 but the initial segment from which the interface grows is off-center, namely, , , , , and . In this case, the tip closest to the right channel wall grows slower and falls behind the other tip. Here, however, the screening is partial in the sense that the slower tip continues to grow but with a velocity that is a fraction of that of the faster tip. (Total screening, whereby the slowest tip eventually stops growing altogether, can be observed if the faster tip has a sufficiently large growth factor.)
The generic Loewner equation given in Eq. (33) can also describe the problem of multiple growing interfaces in the channel geometry. In this case each interface , for , where is the number of distinct interfaces, will be mapped under to a corresponding interval on the real axis in the -plane. Similarly, each advanced interface is mapped by to a polygonal curve in the -plane. One can readily convince oneself that the generic Loewner evolution defined by (33) and (38) applies to this case as well, where is now the total number of vertices corresponding to the sum of the number of vertices for each interface. Note that for each interface we have conditions analogous to (41) and (42), but involving only the growth factors of the respective interface.
In Fig. 9 we show numerical solutions for two growing symmetrical interfaces, with one tip each, where the interfaces are the images of one another with respect to the axis. Here the interfaces start to grow from the intervals and , respectively, with the corresponding tips starting at symmetrical points . (The interfaces have the same growing factors, meaning that .) Notice that, as time goes by, the inner sides of the two interfaces move towards one another leaving a narrow channel between them. For sufficiently large time, the width of such a channel becomes infinitesimally small so that for all practical purposes the resulting evolution will look like a single, symmetrical interface.
If the symmetry between the two interfaces is broken, then one of them will eventually move ahead of the other one. For instance, in Fig. 10 we show the case in which the growth factors are the same as in Fig. 9 but the initial positions of the interfaces are no longer symmetric with respect to the channel centerline. In this case, the interface closer to the centerline is the “winner,” in the sense that its tips moves ahead of the tip of the second interface. This competition between the two growing interfaces resembles the so-called “shadowing effect” in fingering phenomena, whereby the longer fingers grow faster and hinder the growth of the shorter ones in their vicinities.
We have discussed a class of growth models in the channel geometry where the growth dynamics is described in terms of tripolar Loewner evolutions. In the tripolar geometry, the Loewner mapping fixes the points and the point at infinity. The growing domain starts from a point (in the case of a curve) or from an interval (in the case of a domain bounded by an interface) on the real axis, ending at infinity (for infinite time). In particular, a novel class of exact solutions of the tripolar Loewner equation for multiple curves was obtained for the case of stationary driving functions. As for the more difficult problem of interface growth in the channel geometry, we introduced a generalized tripolar Loewner equation and presented several numerical solutions, including the case of an interface with multiple tips as well as the case of multiple interfaces. It was argued that the behavior seen in the model is reminiscent of the phenomenon of tip/finger competition observed in fingering experiments.
A natural extension of the work presented here would be to consider the problem of tripolar stochastic Loewner equation (SLE) and investigate its connection with statistical mechanics models and conformal field theory. Although stochastic growth processes were not considered in this paper, it is nonetheless hoped that a better understanding of deterministic Loewner evolutions in the channel geometry, which was the main thrust of our study, may be of some help in tackling the more difficult problem of tripolar SLEs. Another possible direction for future work would be to consider more general geometries, such as “quadripolar” Loewner evolutions where four points are kept fixed. This geometry would be of interest in connection, for instance, with growth processes in doubly periodic domains.
Acknowledgements.This work was supported in part by the Brazilian agencies FINEP, CNPq, and FACEPE and by the special programs PRONEX and CT-PETRO.
- Löwner, K.: Untersuchungen über schlichte konforme Abbildungen des Einheitskreises. I. Math. Ann. 89, 103–121 (1923).
- Lawler, G.: Introduction to the stochastic Loewner evolution. In: Random Walks and Geometry, pp. 261–294. Walter de Gruyter, Berlin (2004).
- Bauer, M., Bernard, D.: Dipolar stochastic Loewner evolutions, J. Stat. Mech. P03001 (2005). ArXiv:math-ph/0411038.
- Gubiec, T., Szymczak, P.: Fingered growth in the channel geometry: A Loewner-equation approach. Phys. Rev. E 77, 041602-1–041602-12 (2008).
- Schramm, O.: Scaling limits of loop-erased random walks and uniform spanning trees. Isr. J. Math. 118, 221–288 (2000).
- Gruzberg, I. A., Kadanoff, L. P.: The Loewner equation: Maps and shapes. J. Stat. Phys. 114, 1183–1198 (2004).
- Kager, W., Nienhuis, B.: A guide to stochastic Loewner evolution and its applications. J. Stat. Phys. 115, 1149–1229 (2004).
- Cardy, J.: SLE for theoretical physicists. Ann. Phys. 318, 81–118 (2005).
- Bauer, M., Bernard, D.: 2D growth processes: SLE and Loewner chains. Phys. Rep. 432, 115–221 (2006).
- Kager, W., Nienhuis, B., Kadanoff, L. P.: Exact solutions for Loewner evolutions. J. Stat. Phys. 115, 805–822 (2004).
- Selander, G: Two deterministic growth models related to diffusionn limited aggregation. Ph.D. thesis, Royal Institute of Technology, Stockholm (1999).
- Carleson, L., Makarov, N.: Laplacian path models. J. Anal. Math. 87, 103–150 (2002).
- Zabrodin, A.: Growth processes related to the dispersionless Lax equations. Physica D 235, 101–108 (2007).
- Zik, O., Moses, E.: Fingering instability in combustion: An extended view. Phys. Rev. E 60, 518–531 (1999).
- Durán, M. A., Vasconcelos, G. L.: Interface growth in two dimensions: A Loewner-equation approach. Phys. Rev. E 82, 031601-1–03161-8 (2010).
- Zabrodin, A.: Growth of fat slits and dispersionless KP hierarchy. J. Phys. A: Math. Theor. 42, 085206 (2009).
- P. Pelcé, Dynamics of curved fronts (Academic Press, San Diego, 1988).
- G. F. Carrier, M. Krook, and C. E. Pearson, Functions of a complex variable: theory and technique (Hod Books, Ithaca, 1983).