This study investigates the formation and evolution of fishbone patterns in oblique impinging liquid microjets through high-speed imaging experiments and numerical simulations. The results identify periodic oscillations in the upper region of the liquid sheet as the primary mechanism driving fishbone instabilities, which induce rim disturbances and lead to bifurcations into diverse fishbone morphologies. Transitions between stable and unstable flow patterns are systematically mapped across varying Weber numbers and impingement angles, providing a comprehensive framework for understanding this interfacial dynamics. Two critical transitions – marking the onset and disappearance of fishbone patterns – are characterised, offering insights into the underlying physics governing the stability and instability of these flow structures.

A new arbitrary Lagrangian–Eulerian (ALE) formulation for Navier–Stokes flow on self-evolving surfaces is presented. It is based on a general curvilinear surface parameterisation that describes the motion of the ALE frame. Its in-plane part becomes fully arbitrary, while its out-of-plane part follows the material motion of the surface. This allows for the description of flows on deforming surfaces using only surface meshes. The unknown fields are the fluid density or pressure, the fluid velocity and the surface motion, where the latter two share the same normal velocity. The corresponding field equations are the continuity equation or area-incompressibility constraint, the surface Navier–Stokes equations and suitable surface mesh equations. Particularly advantageous are mesh equations based on membrane elasticity. The presentation focuses on the coupled set of strong and weak form equations, and presents several manufactured steady and transient solutions. These solutions are used together with numerical simulations to illustrate and discuss the properties of the proposed new ALE formulation. They also serve as basis for the development and verification of corresponding computational methods. The new formulation allows for a detailed study of fluidic membranes such as soap films, capillary menisci and lipid bilayers.

A liquid film flowing down a fibre becomes unstable, leading to the formation of droplets that travel downstream. The droplet spacing and speed depend on the flow rate for a given nozzle and fibre radii. We show that fibre morphology further modifies the droplet spacing. In particular, we study the effect of the size of the beads in a granular chain on the evolution of the film thickness. We show that, when the size of the bead exceeds a critical value, the selection mechanism for instability modes is modified from regularly spaced droplets to coarsening by droplet merging. Droplet formation for flow over a single bead on the fibre is modified successively over subsequent beads in the downstream. Further, we show that if the perturbation in the flow produced by the bead is introduced as a velocity perturbation at the nozzle inlet, the formation of droplets on the fibre is qualitatively similar to that for the bead.

The rupture of a liquid film, where a thin liquid layer between two other fluids breaks and forms holes, commonly occurs in both natural phenomena and industrial applications. The post-rupture dynamics, from initial hole formation to the complete collapse of the film, are crucial because they govern droplet formation, which plays a significant role in many applications such as disease transmission, aerosol formation, spray drying nanodrugs, oil spill remediation, inkjet printing and spray coating. While single-hole rupture has been extensively studied, the dynamics of multiple-hole ruptures, especially the interactions between neighbouring holes, are less well understood. Here, this study reveals that when two holes ‘meet’ on a curved film, the film evolves into a spinning twisted ribbon before breaking into droplets, distinctly different from what occurs on flat films. We explain the formation and evolution of the spinning twisted ribbon, including its geometry, orbits, corrugations and ligaments, and compare the experimental observations with models. We compare and contrast this phenomena with its counterpart on planar films. While our experiments are based on the multiple-hole ruptures in corona splash, the underlying principles are likely applicable to other systems. This study sheds light on understanding and controlling droplet formation in multiple-hole rupture, improving public health, climate science and various industrial applications.

Ultra-thin liquid sheets generated by impinging two liquid jets are crucial high-repetition-rate targets for laser ion acceleration and ultra-fast physics, and serve widely as barrier-free samples for structural biochemistry. The impact of liquid viscosity on sheet thickness should be comprehended fully to exploit its potential. Here, we demonstrate experimentally that viscosity significantly influences thickness distribution, while surface tension primarily governs shape. We propose a thickness model based on momentum exchange and mass transport within the radial flow, which agrees well with the experiments. These results provide deeper insights into the behaviour of liquid sheets and enable accurate thickness control for various applications, including atomization nozzles and laser-driven particle sources.

In typical atomic force microscopy (AFM) measurements, the AFM probe, mounted on a compliant cantilever, is brought into close proximity to the test substrate. At this range, interfacial attractive van der Waals (vdW) forces can deflect the cantilever by pulling the probe, often causing the probe to suddenly jump into contact with the substrate. For deformable substrates such as gels or bio-tissues, the attraction-induced substrate deformation can further reduce the gap beneath the probe, which can increase the vdW force and hence trigger jump-to-contact prematurely. Since soft gels and tissues are frequently tested in liquid environments, where surface tension and the approaching dynamics of the probe can significantly influence deformation behaviour, this study examines the statics and dynamics of jump-to-contact on elastic substrates incorporating the effect of solid surface tension. We first discuss the theoretical setting for the static problem, deriving perturbation solutions for limiting cases of small and large solid surface tension. Notably, even under conditions of large solid surface tension, elasticity remains critical, as far-field elastic forces are required to smooth surface deformations in a convergent manner. Recognising that practical experiments are inherently dynamic, we also analyse the role of hydrodynamic pressure, which can delay the premature jump-to-contact. Our analysis focuses on identifying the conditions under which dynamic effects are negligible, enabling the simple analytical solutions in the static problem to reliably interpret AFM experimental results.

Viscous fingering instabilities, common in confined environments such as porous media or Hele-Shaw cells, surprisingly also occur in unconfined, non-porous settings as revealed by recent experiments. These novel instabilities involve free-surface flows of dissimilar viscosity. We demonstrate that such a free-surface flow, involving a thin film of viscous fluid spreading over a substrate that is prewetted with a fluid of higher viscosity, is susceptible to a similar type of novel viscous fingering instability. Such flows are relevant to a range of geophysical, industrial and physiological applications from the small scales of thin-film coating applications and nasal drug delivery to the large scales of lava flows. In developing a theoretical framework, we assume that the intruding layer and the liquid film over which it flows are both long and thin, the effects of inertia and surface tension are negligible, and both layers are driven by gravity and resisted by viscous shear stress so that the principles of lubrication theory hold. We investigate the stability of axisymmetric similarity solutions, describing the base flow, by examining the growth of small-amplitude non-axisymmetric perturbations. We characterise regions of instability across parameter space and find that these instabilities emerge above a critical viscosity ratio. That is, a fluid of low viscosity intruding into another fluid of sufficiently high viscosity is susceptible to instability, akin to traditional viscous fingering in a porous medium. We identify the mechanism of instability, compare with other frontal instabilities and demonstrate that high enough density differences suppress the instability completely.

The phenomenon of bulge evolution under the action of gravity on shallow water is prevalent both in natural occurrences and engineering industries. However, despite its ubiquity, its physical process remains largely unexplored. The evolution of bulge contains two fundamental physical processes: collapse and propagation. The collapse process can be further divided into two sub-processes: squeezing process and diffusion process. Based on the weakly nonlinear shallow water assumption with the classical perturbation method, the governing equations controlling the surface elevations in the diffusion process and the propagation process have been theoretically derived, where a bulge-induced surface pressure is modeled for the propagation process. Moreover, their scaling laws for the decay of wave height are also established, which have been validated by direct numerical simulation results. The derived scaling laws for wave height attenuation of bulge evolution provide profound insights, which hold the potential to applications in the engineering industry.

We investigate experimentally the effect of salinity and atmospheric humidity on the drainage and lifetime of thin liquid films motivated by conditions relevant to air–sea exchanges. We show that the drainage is independent of humidity and that the effect of a change in salinity is reflected only through the associated change in viscosity. On the other hand, film lifetime displays a strong dependence on humidity, with more than a tenfold increase between low and high humidities: from a few seconds to tens of minutes. Mixing the air surrounding the film also has a very important effect on lifetime, modifying its distribution and reducing the mean lifetime of the film. From estimations of the evaporation rate, we are able to derive scaling laws that describe well the evolution of lifetime with a change of humidity. Observations of the black film, close to the top where the film ruptures, reveal that this region is very sensitive to local humidity conditions.

Deformation occurs in a thin liquid film when it is subjected to a non-uniform electric field, which is referred to as the electrohydrodynamic patterning. Due to the development of a non-uniform electrical force along the surface, the film would evolve into microstructures/nanostructures. In this work, a linear and a nonlinear model are proposed to thoroughly investigate the steady state (i.e. equilibrium state) of the electrohydrodynamic deformation of thin liquid film. It is found that the deformation is closely dependent on the electric Bond number BoE. Interestingly, when BoE is larger than a critical value, the film would be deformed remarkably and get in contact with the top template. To model the ‘contact’ between the liquid film and the solid template, the disjoining pressure is incorporated into the numerical model. From the nonlinear numerical model, a hysteresis deformation is revealed, i.e. the film may have different equilibrium states depending on whether the voltage is increased or decreased. To analyse the stability of these multiple equilibrium states, the Lyapunov functional is employed to characterise the system’s free energy. According to the Lyapunov functional analysis, at most three equilibrium states can be formed. Among them, one is stable, another is metastable and the third one is unstable. Finally, the model is extended to study the three-dimensional deformation of the electrohydrodynamic patterning.

The influence of parametric forcing on a viscoelastic fluid layer, in both gravitationally stable and unstable configurations, is investigated via linear stability analysis. When such a layer is vertically oscillated beyond a threshold amplitude, large interface deflections are caused by Faraday instability. Viscosity and elasticity affect the damping rate of momentary disturbances with arbitrary wavelength, thereby altering the threshold and temporal response of this instability. In gravitationally stable configurations, calculations show that increased elasticity can either stabilize or destabilize the viscoelastic system. In weakly elastic liquids, higher elasticity increases damping, raising the threshold for Faraday instability, whereas the opposite is observed in strongly elastic liquids. While oscillatory instability occurs in Newtonian fluids for all gravity levels, we find that parametric forcing below a critical frequency will cause a monotonic instability for viscoelastic systems at microgravity. Importantly, in gravitationally unstable configurations, parametric forcing above this frequency stabilizes viscoelastic fluids, until the occurrence of a second critical frequency. This result contrasts with the case of Newtonian liquids, where under the same conditions, forcing stabilizes a system for all frequencies below a single critical frequency. Analytical expressions are obtained under the assumption of long wavelength disturbances predicting the damping rate of momentary disturbances as well as the range of parameters that lead to a monotonic response under parametric forcing.

A linear stability analysis of a soluble surfactant-laden liquid film flowing down a compliant substrate is performed. Our purpose is to expand the prior studies (Carpenter and Garrad 1985 J. Fluid Mech. 155, 465–510; Alexander et al., 2020 J. Fluid Mech. 900, A40) by incorporating a soluble surfactant into the flow configuration. As a result, we formulate the Orr–Sommerfeld-type boundary value problem and solve it analytically by using the long-wave series expansion as well as numerically by using the Chebyshev spectral collocation method in an arbitrary wavenumber regime for infinitesimal disturbances. The long-wave result reveals that surface instability is stabilized in the presence of a surfactant, whereas it is destabilized in the presence of a compliant substrate. These opposing impacts suggest an analytical relationship between parameters associated with the soluble surfactant and compliant wall, ensuring the same critical Reynolds number for the emergence of surface instability corresponding to both surfactant-laden film flow over a compliant wall and surfactant-free film flow over a non-compliant wall. In the arbitrary wavenumber regime, along with the surface mode, we identify two additional modes based on their distinct phase speeds. Specifically, the wall mode emerges in the finite wavenumber regime, while the shear mode emerges only when the Reynolds number is large. As the surfactant Marangoni number increases, the wall mode destabilizes, resulting in a different outcome from the surface mode. Moreover, increasing the value of the ratio of adsorption and desorption rate constants stabilizes surface instability but destabilizes wall mode instability. As a result, we perceive that the soluble surfactant-laden film flow is linearly more unstable than the insoluble one due to surface instability but linearly more stable than the insoluble one due to wall mode instability. Additionally, we see a peculiar behaviour of base surface surfactant concentration on the primary instability. In fact, it has a specific value depending on adsorption and desorption rate constants below which surface instability stabilizes but wall mode instability destabilizes, whereas above which an opposite phenomenon occurs. Finally, in the high-Reynolds-number regime, we can suppress shear mode instability by raising the surfactant Marangoni number and the ratio of adsorption and desorption rate constants when the angle of inclination is sufficiently small. Unlike surface instability, the base surface surfactant concentration exhibits both stabilizing and destabilizing influences on shear mode instability.

We investigate the dynamics of close-contact melting (CCM) on ‘gas-trapped’ hydrophobic surfaces, with specific focus on the effects of geometrical confinement and the liquid–air meniscus below the liquid film. By employing dual-series and perturbation methods under the assumption of small meniscus deflections, we obtain numerical solutions for the effective slip lengths associated with velocity $\lambda$ and temperature $\lambda _t$ fields, across various values of aspect ratio $\Lambda$ (defined as the ratio of the film thickness $h$ to the structure’s periodic length $l$) and gas–liquid fraction $\phi$. Asymptotic solutions of $\lambda$ and $\lambda _t$ for $\Lambda \ll 1$ and $\Lambda \gg 1$ are derived and summarised for different surface structures, interface shapes and $\Lambda$, which reveal a different trend of $\lambda$ for $\Lambda \ll 1$ and depending on the presence of a meniscus. In the context of constant-pressure CCM, our results indicate that longitudinal grooves can enhance heat transfer under the effects of confinement and a meniscus when $\Lambda \lesssim 0.1$ and $\phi \lt 1 - 0.5^{2/3} \approx 0.37$. For gravity-driven CCM, the parameters of $l$ and $\phi$ determine whether the melting rate is enhanced, reduced or nearly unaffected. We construct a phase diagram based on the parameter matrix $(\log _{10} l, \phi )$ to delineate these three regimes. Lastly, we derive two asymptotic solutions for predicting the variation in time of the unmelted solid height.

Surfactant transport is central to a diverse range of natural phenomena with numerous practical applications in physics and engineering. Surprisingly, this process remains relatively poorly understood at the molecular scale. Here, we use non-equilibrium molecular dynamics (NEMD) simulations to study the spreading of sodium dodecyl sulphate on a thin film of liquid water. The molecular form of the control volume is extended to a coordinate system moving with the liquid–vapour interface to track surfactant spreading. We use this to compare the NEMD results to the continuum description of surfactant transport on an interface. By including the molecular details in the continuum model, we establish that the transport equation preserves substantial accuracy in capturing the underlying physics. Moreover, the relative importance of the different mechanisms involved in the transport process is identified. Consequently, we derive a novel exact molecular equation for surfactant transport along a deforming surface. Close agreement between the two conceptually different approaches, i.e. NEMD simulations and the numerical solution of the continuum equation, is found as measured by the surfactant concentration profiles, and the time dependence of the so-called spreading length. The current study focuses on a relatively simple specific solvent–surfactant system, and the observed agreement with the continuum model may not arise for more complicated industrially relevant surfactants and anti-foaming agents. In such cases, the continuum approach may fail to predict accompanying phase transitions, which can still be captured through the NEMD framework.

We consider the dynamics of a liquid film with a pinned contact line (for example, a drop), as described by the one-dimensional, surface-tension-driven thin-film equation $h_t + (h^n h_{xxx})_x = 0$, where $h(x,t)$ is the thickness of the film. The case $n=3$ corresponds to a film on a solid substrate. We derive an evolution equation for the contact angle $\theta (t)$, which couples to the shape of the film. Starting from a regular initial condition $h_0(x)$, we investigate the dynamics of the drop both analytically and numerically, focusing on the contact angle. For short times $t\ll 1$, and if $n\ne 3$, the contact angle changes according to a power law $\displaystyle t^{\frac {n-2}{4-n}}$. In the critical case $n=3$, the dynamics become non-local, and $\dot {\theta }$ is now of order $\displaystyle {\rm{e}}^{-3/(2t^{1/3})}$. This implies that, for $n=3$, the standard contact line problem with prescribed contact angle is ill posed. In the long time limit, the solution relaxes exponentially towards equilibrium.

We study the dynamics of a thin liquid sheet that flows upwards along the sides of a vertically aligned, impacting plate. Upon impact of the vertical solid plate onto a liquid pool, the liquid film is ejected and subsequently continues to flow over the solid surface while the plate enters the water. With increasing impact velocity, the liquid film is observed to rise up faster and higher. We focus on the time evolution of the liquid film height and the thickness of its upper rim and discuss their dynamics in detail. Similar to findings in previous literature on sheet fragmentation during drop impact, we find the rim thickness to be governed by the local instantaneous capillary number based on gravity and the deceleration of the liquid sheet, showing that the retraction of the rim is primarily due to capillarity. In contrast, for the liquid film height, we demonstrate that the viscous dissipation in the thin boundary layer is the dominant factor for the vertical deceleration of the liquid sheet, by modelling the time evolution of the film height and showing that the influences of capillarity, gravity and deceleration due to the air phase are all negligible compared with the viscous term. Finally, we introduce characteristic viscous time and length scales based on the initial rim thickness and show that the maximum height of the film and the corresponding time can be determined from these viscous scales.

High-power laser systems require thin films with extremely low absorption. Ultra-low-absorption films are often fabricated via ion beam sputtering, which is costly and slow. This study analyzes the impact of doping titanium and annealing on the absorption characteristics of thin films, focusing on composition and structure. The results indicate that the primary factor influencing absorption is composition. Suppressing the presence of electrons or holes that do not form stable chemical bonds can significantly reduce absorption; for amorphous thin films, the structural influence on absorption is relatively minor. Thus, composition control is crucial for fabricating ultra-low-absorption films, while the deposition method is secondary. Ion beam-assisted electron-beam evaporation, which is relatively seldom used for fabricating low-absorption films, was employed to produce high-reflectivity films. After annealing, the absorption at 1064 nm reached 1.70 parts per million. This method offers a cost-effective and rapid approach for fabricating ultra-low-absorption films.

We investigate the sliding dynamics of a millimetre-sized particle trapped in a horizontal soap film. Once released, the particle moves toward the centre of the film in damped oscillations. We study experimentally and model the forces acting on the particle, and evidence the key role of the mass of the film on the shape of the film and particle dynamics. Not only is the gravitational distortion of the film measurable, it completely determines the force responsible for the motion of the particle – the catenoid-like deformation induced by the particle has negligible effect on the dynamics. Surprisingly, this is expected for all film sizes as long as the particle radius remains much smaller than the film width. We also measure the friction force, and show that ambient air and the film contribute almost equally to the friction. The theoretical model that we propose predicts exactly the friction coefficient as long as inertial effects can be neglected in air (for the smallest and slowest particles). The fit between theory and experiments sets an upper boundary $\eta _s \leqslant 10^{-8}$ Pa s m for the surface viscosity, in excellent agreement with recent interfacial microrheology measurements.

We study the dynamics of fracture deflation following hydraulic fracturing of an infinite elastic solid, with fluid removal from a narrow conduit at the centre. This process involves coupled lubricating flow and elastic deformation, now subject to appropriate descriptions of fluid removal through the conduit towards the ambient, driven by elastic stresses and extraction/suction. When the influence of material toughness is negligible, the dynamics is found to be governed by two dimensionless parameters that describe the relative influence of elasticity-driven backflow ($\Pi _c$) and ambient-pressure-driven backflow ($\Pi _e$), respectively. We also found that the fracture’s thickness eventually approaches zero at the centre, while the fracture evolves into a self-similar shape of the dipole type that conserves the dipole moment $M$. The fracture’s front continues to elongate according to $x_f \propto t^{1/9}$, while the total fluid volume within the fracture decreases according to $V \propto t^{-1/9}$. The model and solutions might find use in practical problems to estimate the rate of backflow and effective permeability of a fractured reservoir once pressure is released.

Spin coating is the process of generating a uniform coating film on a substrate by centrifugal forces during rotation. In the framework of lubrication theory, we investigate the axisymmetric film evolution and contact-line dynamics in spin coating on a partially wetting substrate. The contact-line singularity is regularized by imposing a Navier slip model. The interface morphology and the contact-line movement are obtained by numerical solution and asymptotic analysis of the lubrication equation. The results show that the evolution of the liquid film can be classified into two modes, depending on the rotational speed. At lower speeds, the film eventually reaches an equilibrium state, and we provide a theoretical description of how the equilibrium state can be approached through matched asymptotic expansions. At higher speeds, the film exhibits two or three distinct regions: a uniform thinning film in the central region, an annular ridge near the contact line, and a possible Landau–Levich–Derjaguin-type (LLD-type) film in between that has not been reported previously. In particular, the LLD-type film occurs only at speeds slightly higher than the critical value for the existence of the equilibrium state, and leads to the decoupling of the uniform film and the ridge. It is found that the evolution of the ridge can be well described by a two-dimensional quasi-steady analysis. As a result, the ridge volume approaches a constant and cannot be neglected to predict the variation of the contact-line radius. The long-time behaviours of the film thickness and the contact radius agree with derived asymptotic solutions.