University of Twente

Soft elastic interfaces can strongly deform under the influence of external forces, and can even exhibit elastic singularities. Here we discuss two cases where such singularities occur. First, we describe surface creases that form under compression (or swelling) of an elastic medium. Second, we consider the elastocapillary ridges that form when a soft substrate is wetted by a liquid drop. Analytical descriptions are presented and compared to experiments. We reveal that, like for liquid interfaces, the surface tension of the solid is a key factor in shaping the surface, and determines the nature of the singularity.

When two drops come into contact they will rapidly merge and form a single drop. Here we address the coalescence of drops on a substrate, focussing on the initial dynamics just after contact. For very viscous drops we present similarity solutions for the bridge that connects the two drops, the size of which grows linearly with time. Both the dynamics and the self-similar bridge profiles are verified quantitatively by experiments. We then consider the coalescence of water drops, for which viscosity can be neglected and liquid inertia takes over. Once again, we find that experiments display a self-similar dynamics, but now the bridge size grows with a power-law $t^{2/3}$. We provide a scaling theory for this behavior, based on geometric arguments. The main result for both viscous and inertial drops is that the contact angle is important as it determines the geometry of coalescence -- yet, the contact line dynamics appears irrelevant for the early stages of coalescence.

Exponential time integrators are a powerful tool for numerical solution

of time dependent problems. The actions of the matrix functions on vectors,

necessary for exponential integrators, can be efficiently computed by

different elegant numerical techniques, such as Krylov subspaces.

Unfortunately, in some situations the additional work required by

exponential integrators per time step is not paid off because the time step

can not be increased too much due to the accuracy restrictions.

To get around this problem, we propose the so-called time-stepping-free

approach. This approach works for linear ordinary differential equation (ODE)

systems where the time dependent part forms a small-dimensional subspace.

In this case the time dependence can be projected out by block Krylov

methods onto the small, projected ODE system. Thus, there is just one

block Krylov subspace involved and there are no time steps. We refer to

this method as EBK, exponential block Krylov method. The accuracy of EBK

is determined by the Krylov subspace error and the solution accuracy in the

projected ODE system. EBK works for well for linear systems, its extension

to nonlinear problems is an open problem and we discuss possible ways for

such an extension.