Michael C. Dallaston

Accurate series solutions for gravity-driven Stokes waves

In the past, high order series expansion techniques have been used to study the nonlinear equations that govern the form of periodic Stokes waves moving steadily on the surface of an inviscid fluid. In the present study, two such series solutions are recomputed using exact arithmetic, eliminating any loss of accuracy due to accumulation of round-off error, allowing a much greater number of terms to be found with confidence. It is shown that a higher order behavior of the series generated by the solution casts doubt over arguments that rely on estimating the series’ radius of convergence. Further, the exact nature of the series is used to shed light on the unusual nature of convergence of higher order Pade approximants near the highest wave. Finally, it is concluded that, provided exact values are used in the series, these Pade approximants prove very effective in successfully predicting three turning points in both the dispersion relation and the total energy.

New exact solutions for Hele-Shaw flow in doubly connected regions

Radial Hele–Shaw flows are treated analytically using conformal mapping techniques. The geometry of interest has a doubly connected annular region of viscous fluid surrounding an inviscid bubble that is either expanding or contracting due to a pressure difference caused by injection or suction of the inviscid fluid. The zero-surface-tension problem is ill-posed for both bubble expansion and contraction, as both scenarios involve viscous fluid displacing inviscid fluid. Exact solutions are derived by tracking the location of singularities and critical points in the analytic continuation of the mapping function. We show that by treating the critical points, it is easy to observe finite-time blow-up, and the evolution equations may be written in exact form using complex residues. We present solutions that start with cusps on one interface and end with cusps on the other, as well as solutions that have the bubble contracting to a point. For the latter solutions, the bubble approaches an ellipse in shape at ex...

Self-similar finite-Time singularity formation in degenerate parabolic equations arising in thin-film flows

A thin liquid film coating a planar horizontal substrate may be unstable to perturbations in the film thickness due to unfavourable intermolecular interactions between the liquid and the substrate, which may lead to finite-time rupture. The self-similar nature of the rupture has been studied before by utilising the standard lubrication approximation along with the Derjaguin (or disjoining) pressure formalism used to account for the intermolecular interactions, and a particular form of the disjoining pressure with exponent n  =  3 has been used, namely, Π(h)∝−1/h3where h is the film thickness. In the present study, we use a numerical continuation method to compute discrete solutions to self-similar rupture for a general disjoining pressure exponent n (not necessarily equal to 3), which has not been previously performed. We focus on axisymmetric point-rupture solutions and show for the first time that pairs of solution branches merge as n decreases, starting at nc≈1.485. We verify that this observation also holds true for plane-symmetric line-rupture solutions for which the critical value turns out to be slightly larger than for the axisymmetric case, nplanec≈1.499. Computation of the full time-dependent problem also demonstrates the loss of stable similarity solutions and the subsequent onset of cascading, increasingly small structures.

Corner and finger formation in Hele-Shaw flow with kinetic undercooling regularisation

We examine the effect of a kinetic undercooling condition on the evolution of a free boundary in Hele--Shaw flow, in both bubble and channel geometries. We present analytical and numerical evidence that the bubble boundary is unstable and may develop one or more corners in finite time, for both expansion and contraction cases. This loss of regularity is interesting because it occurs regardless of whether the less viscous fluid is displacing the more viscous fluid, or vice versa. We show that small contracting bubbles are described to leading order by a well-studied geometric flow rule. Exact solutions to this asymptotic problem continue past the corner formation until the bubble contracts to a point as a slit in the limit. Lastly, we consider the evolving boundary with kinetic undercooling in a Saffman--Taylor channel geometry. The boundary may either form corners in finite time, or evolve to a single long finger travelling at constant speed, depending on the strength of kinetic undercooling. We demonstrate these two different behaviours numerically. For the travelling finger, we present results of a numerical solution method similar to that used to demonstrate the selection of discrete fingers by surface tension. With kinetic undercooling, a continuum of corner-free travelling fingers exists for any finger width above a critical value, which goes to zero as the kinetic undercooling vanishes. We have not been able to compute the discrete family of analytic solutions, predicted by previous asymptotic analysis, because the numerical scheme cannot distinguish between solutions characterised by analytic fingers and those which are corner-free but non-analytic.

Discrete self-similarity in interfacial hydrodynamics and the formation of iterated structures

The formation of iterated structures, such as satellite and subsatellite drops, filaments, and bubbles, is a common feature in interfacial hydrodynamics. Here we undertake a computational and theoretical study of their origin in the case of thin films of viscous fluids that are destabilized by long-range molecular or other forces. We demonstrate that iterated structures appear as a consequence of discrete self-similarity, where certain patterns repeat themselves, subject to rescaling, periodically in a logarithmic time scale. The result is an infinite sequence of ridges and filaments with similarity properties. The character of these discretely self-similar solutions as the result of a Hopf bifurcation from ordinarily self-similar solutions is also described.

A curve shortening flow rule for closed embedded plane curves with a prescribed rate of change in enclosed area

Motivated by a problem from fluid mechanics, we consider a generalization of the standard curve shortening flow problem for a closed embedded plane curve such that the area enclosed by the curve is forced to decrease at a prescribed rate. Using formal asymptotic and numerical techniques, we derive possible extinction shapes as the curve contracts to a point, dependent on the rate of decreasing area; we find there is a wider class of extinction shapes than for standard curve shortening, for which initially simple closed curves are always asymptotically circular. We also provide numerical evidence that self-intersection is possible for non-convex initial conditions, distinguishing between pinch-off and coalescence of the curve interior.

A tunable high-pass filter for simple and inexpensive size-segregation of sub-10-nm nanoparticles

Recent advanced in the fields of nanotechnology and atmospheric sciences underline the increasing need for sizing sub-10-nm aerosol particles in a simple yet efficient way. In this article, we develop, experimentally test and model the performance of a High-Pass Electrical Mobility Filter (HP-EMF) that can be used for sizing nanoparticles suspended in gaseous media. Experimental measurements of the penetration of nanoparticles having diameters down to ca 1nm through the HP-EMF are compared with predictions by an analytic, a semi-empirical and a numerical model. The results show that the HP-EMF effectively filters nanoparticles below a threshold diameter with an extremely high level of sizing performance, while it is easier to use compared to existing nanoparticle sizing techniques through design simplifications. What is more, the HP-EMF is an inexpensive and compact tool, making it an enabling technology for a variety of applications ranging from nanomaterial synthesis to distributed monitoring of atmospheric nanoparticles.

Free-boundary models of a meltwater conduit

We analyse the cross-sectional evolution of an englacial meltwater conduit that contracts due to inward creep of the surrounding ice and expands due to melting. Making use of theoretical methods from free-boundary problems in Stokes flow and Hele–Shaw squeeze flow we construct an exact solution to the coupled problem of external viscous creep and internal heating, in which we adopt a Newtonian approximation for ice flow and an idealized uniform heat source in the conduit. This problem provides an interesting variant on standard free-boundary problems, coupling different internal and external problems through the kinematic condition at the interface. The boundary in the exact solution takes the form of an ellipse that may contract or expand (depending on the magnitudes of effective pressure and heating rate) around fixed focal points. Linear stability analysis reveals that without the melting this solution is unstable to perturbations in the shape. Melting can stabilize the interface unless the aspect ratio is too small; in that case, instabilities grow largest at the thin ends of the ellipse. The predictions are corroborated with numerical solutions using boundary integral techniques. Finally, a number of extensions to the idealized model are considered, showing that a contracting circular conduit is unstable to all modes of perturbation if melting occurs at a uniform rate around the boundary, or if the ice is modelled as a shear-thinning fluid.