Flow of evaporating, gravity-driven thin liquid films over topography

The effect of topography on the free surface and solvent concentration profiles of an evaporating thin film of liquid flowing down an inclined plane is considered. The liquid is assumed to be composed of a resin dissolved in a volatile solvent with the associated solvent concentration equation derived on the basis of the well-mixed approximation. The dynamics of the film is formulated as a lubrication approximation and the effect of a composition-dependent viscosity is included in the model. The resulting time-dependent, nonlinear, coupled set of governing equations is solved using a full approximation storage multigrid method. The approach is first validated against a closed-form analytical solution for the case of a gravity-driven, evaporating thin film flowing down a flat substrate. Analysis of the results for a range of topography shapes reveal that although a full-width, spanwise topography such as a step-up or a step-down does not affect the composition of the film, the same is no longer true for the case of localized topography, such as a peak or a trough, for which clear nonuniformities of the solvent concentration profile can be observed in the wake of the topography

A numerical investigation of the solution of a class of fourth-order eigenvalue problems

This paper is concerned with the accurate numerical approximation of the spectral properties of the biharmonic operator on various domains in two dimensions. A number of analytic results concerning the eigenfunctions of this operator are summarized and their implications for numerical approximation are discussed. In particular, the asymptotic behaviour of the first eigenfunction is studied since it is known that this has an unbounded number of oscillations when approaching certain types of corners on domain boundaries. Recent computational results of Bjorstad & Tjostheim, using a highly accurate spectral Legendre-Galerkin method, have demonstrated that a number of these sign changes may be accurately computed on a square domain provided sufficient care is taken with the numerical method. We demonstrate that similar accuracy is also achieved using an unstructured finite-element solver which may be applied to problems on domains with arbitrary geometries. A number of results obtained from this mixed finite-element approach are then presented for a variety of domains. These include a family of circular sector regions, for which the oscillatory behaviour is studied as a function of the internal angle, and another family of (symmetric and non-convex) domains, for which the parity of the least eigenfunction is investigated. The paper not only verifies existing asymptotic theory, but also allows us to make a new conjecture concerning the eigenfunctions of the biharmonic operator

A combined experimental and computational fluid dynamics analysis of the dynamics of drop formation

This article presents a complementary experimental and computational investigation of the effect of viscosity and flowrate on the dynamics of drop formation in the dripping mode. In contrast to previous studies, numerical simulations are performed with two popular commercial computational fluid dynamics (CFD) packages, CFX and FLOW-3D, both of which employ the volume of fluid (VOF) method. Comparison with previously published experimental and computational data and new experimental results reported here highlight the capabilities and limitations of the aforementioned packages

A Moving-Mesh Finite Element Method and its Application to the Numerical Solution of Phase-Change Problems

A distributed Lagrangian moving-mesh finite element method is applied to problems involving changes of phase. The algorithm uses a distributed conservation principle to determine nodal mesh velocities, which are then used to move the nodes. The nodal values are obtained from an ALE (Arbitrary Lagrangian-Eulerian) equation, which represents a generalisation of the original algorithm presented in Applied Numerical Mathematics, 54:450–469 (2005). Having described the details of the generalised algorithm it is validated on two test cases from the original paper and is then applied to one-phase and, for the first time, twophase Stefan problems in one and two space dimensions, paying particular attention to the implementation of the interface boundary conditions. Results are presented to demonstrate the accuracy and the effectiveness of the method, including comparisons against analytical solutions where available.

Quantitative phase-field modeling of solidification at high Lewis number

A phase-field model of nonisothermal solidification in dilute binary alloys is used to study the variation of growth velocity, dendrite tip radius, and radius selection parameter as a function of Lewis number at fixed undercooling. By the application of advanced numerical techniques, we have been able to extend the analysis to Lewis numbers of order 10 000, which are realistic for metals. A large variation in the radius selection parameter is found as the Lewis number is increased from 1 to 10 000

Simulations of three-dimensional dendritic growth using a coupled thermo-solutal phase-field model

Using a phase field model, which fully couples the thermal and solute concentration field, we present simulation results in three dimensions of the rapid dendritic solidification of a class of dilute alloys at the meso scale. The key results are the prediction of steady state tip velocity and radius at varying undercooling and thermal diffusivities. Less computationally demanding 2-dimensional results are directly compared with the corresponding 3-dimensional results, where significant quantitative differences emerge. The simulations provide quantitative predictions for the range of thermal and solutal diffusivities considered and show the effectiveness and potential of the computational techniques employed. These results thus provide benchmark 3-dimensional computations, allow direct comparison with underlying analytical theory, and pave the way for further quantitative results

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