Search

In previous work, we have performed amplitude expansions of the continuum equations for the Grinfeld instability and carried them to high orders. Nevertheless, the approach turned out to be restricted to relatively small amplitudes. In this article, we use a variational approach in terms of multi-cycloid curves instead. Besides its higher precision at given order, the method has the advantages of giving a transparent physical meaning to the appearance of cusp singularities and of not being restricted to interfaces representable as single-valued functions. Using a single cycloid as ansatz function, the entire calculation can be performed analytically, which gives a good qualitative overview of the system. Taking into account several but few cycloid modes, we obtain remarkably good quantitative agreement with previous numerical calculations. With a few more modes taken into consideration, we improve on the accuracy of those calculations. Our approach extends them to situations involving gravity effects. Results on the shape of steady-state solutions are presented at both large stresses and amplitudes. In addition, their stability is investigated., Comment: subm. to EPJ B

Amplitude expansions are used to determine steady states of a semi-infinite solid subject to the Grinfeld instability in systems with a fixed (wave)length. We present two methods to obtain high-order weakly nonlinear results. Using the system size as a control parameter, we circumvent the problem that there is no instability threshold for an extended system in the absence of gravity. This way, the case without gravity becomes accessible to a weakly nonlinear treatment. The dependence of the branch structure of solution space on the level of gravity (or density difference) is exhibited. In the zero-gravity limit, we recover the solution branch obtained by Spencer and Meiron. A transition from a supercritical to a subcritical bifurcation is observed as gravity is increased or the nonhydrostatic stress is decreased at fixed gravity. At given values of the system parameters, we find a discrete, possibly infinite, set of solution branches. This is reminiscent of dendritic or eutectic growth, where similar solution sets exist, of which only a particular one is linearly stable. Despite the high order of our expansions, the approach is restricted to relatively small nondimensional amplitudes ($\lesssim 0.2$), a disadvantage we can overcome by a variational approach that will be discussed in a companion paper. At the critical point, we find that not only the first Landau coefficient is negative but all of them up to the highest amplitude order (15) we could compute so far., Comment: This article has been withdrawn by the authors. The authors determined that the results claimed were not as complete as originally thought, and the project has ended so no future work on this is planned

A phase-field approach describing the dynamics of a strained solid in contact with its melt is developed. By rigorous asymptotic analysis we show that the sharp-interface limit of this model recovers the continuum model equations for the Grinfeld instability. Moreover, we use our approach to derive hitherto unknown sharp-interface equations for a situation including a field of body forces. The numerical utility of the phase-field approach is demonstrated by comparison with a sharp-interface simulation. We then investigate the dynamics of extended systems within the phase-field model which contains an inherent lower length cutoff, thus avoiding cusp singularities. It is found that a periodic array of grooves generically evolves into a superstructure which arises from a series of imperfect period doublings. For wavenumbers close to the fastest-growing mode of the linear instability, the first period doubling can be obtained analytically. Both the dynamics of an initially periodic array and a random initial structure can be described as a coarsening process with winning grooves temporarily accelerating whereas losing ones decelerate and even reverse their direction of motion., Comment: 40 pages, 12 figures, submitted to Phys. Rev. E