Large amplitude behavior of the Grinfeld instability, Part I: High-order   weakly nonlinear analysis
  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. 
Peter Kohlert, Klaus Kassner and Chaouqi Misbah