Finite element simulation of three-dimensional free-surface flow problems

An adaptive finite element algorithm is described for the stable solution of three-dimensional free-surface-flow problems based primarily on the use of node movement. The algorithm also includes a discrete remeshing procedure which enhances its accuracy and robustness. The spatial discretisation allows an isoparametric piecewise-quadratic approximation of the domain geometry for accurate resolution of the curved free surface. The technique is illustrated through an implementation for surface-tension-dominated viscous flows modelled in terms of the Stokes equations with suitable boundary conditions on the deforming free surface. Two three-dimensional test problems are used to demonstrate the performance of the method: a liquid bridge problem and the formation of a fluid droplet

Viscous sintering of unimodal and bimodal cylindrical packings with shrinking pores

Published versio

Conformal mapping methods for interfacial dynamics

The article provides a pedagogical review aimed at graduate students in materials science, physics, and applied mathematics, focusing on recent developments in the subject. Following a brief summary of concepts from complex analysis, the article begins with an overview of continuous conformal-map dynamics. This includes problems of interfacial motion driven by harmonic fields (such as viscous fingering and void electromigration), bi-harmonic fields (such as viscous sintering and elastic pore evolution), and non-harmonic, conformally invariant fields (such as growth by advection-diffusion and electro-deposition). The second part of the article is devoted to iterated conformal maps for analogous problems in stochastic interfacial dynamics (such as diffusion-limited aggregation, dielectric breakdown, brittle fracture, and advection-diffusion-limited aggregation). The third part notes that all of these models can be extended to curved surfaces by an auxilliary conformal mapping from the complex plane, such as stereographic projection to a sphere. The article concludes with an outlook for further research.Comment: 37 pages, 12 (mostly color) figure

On the calculation of normals in free-surface flow problems

The use of boundary-conforming finite element methods is considered for the solution of surface-tension-dominated free-surface flow problems in three dimensions. This class of method is based upon the use of a moving mesh whose velocity is driven by the motion of the free surface, which is in turn determined via a kinematic boundary condition for the normal velocity. The significance of the method used to compute the normal direction at the finite element node points for a C0 piecewise-polynomial free surface is investigated. In particular, it is demonstrated that the concept of mass-consistent normals on an isoparametric quadratic tetrahedral mesh is flawed. In this case an alternative, purely geometric, normal is shown to lead to a far more robust numerical algorithm

Exact solution to a class of functional difference equations with application to a moving contact line flow

A new integral representation for the Barnes double gamma function is derived. This is canonical in the sense that solutions to a class of functional difference equations of first order with trigonometrical coefficients can be expressed in terms of the Barnes function. The integral representation given here makes these solutions very simple to compute. Several well-known difference equations are solved by this method and a treatment of the linear theory for moving contact line flow in an inviscid fluid wedge is given.https://journals.cambridge.org/action/displayAbstract?aid=231870

Adaptive finite element simulation of three-dimensional surface tension dominated free-surface flow problems

An arbitrary Lagrangian--Eulerian finite element method is described for the solution of time-dependent, three-dimensional, free-surface flow problems. Many flows of practical significance involve contact lines, where the free surface meets a solid boundary. This contact line may be pinned to a particular part of the solid but is more typically free to slide in a manner that is characterised by the dynamic contact angle formed by the fluid. We focus on the latter case and use a model that admits spatial variation of the contact angle: thus permitting variable wetting properties to be simulated. The problems are driven by the motion of the fluid free surface (under the action of surface tension and external forces such as gravity) hence the geometry evolves as part of the solution, and mesh adaptivity is required to maintain the quality of the computational mesh for the physical domain. Continuous mesh adaptivity, in the form of a pseudo-elastic mesh movement scheme, is used to move the interior mesh nodes in response to the motion of the fluid's free surface. Periodic, discrete remeshing stages are also used for cases in which the fluid volume has grown, or is sufficiently distorted, by the free-surface motion. Examples are given of a droplet sliding on an inclined uniform plane and of a droplet spreading on a surface with variable wetting properties