Source author record

Stéphane Zaleski

Stéphane Zaleski appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

physics.flu-dyn math.AP physics.comp-ph

Catalog footprint

What is connected

9works

3topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2026arXiv

Applying contact angles based on a continuously moving contact line in 2D VOF simulations

A height-function-based numerical approach is developed for enforcing contact angles on flat and curved solid surfaces within two-dimensional volume-of-fluid simulations. This method incorporates the contact line position into the curvature estimation in the contact line cell, where the interface normal is constrained to the prescribed angle to ensure smooth contact line motion. On flat solid surfaces, the proposed model achieves higher accuracy than the conventional vertical height-function method for enforcing very small and very large contact angles. Moreover, it extends naturally to curved solid surfaces represented by the embedded boundary method, enabling the imposition of arbitrary contact angles without inducing significant spurious currents near the contact line. Numerical validations confirm the accuracy and robustness of the proposed method, demonstrating its applicability to multiphase flows involving complex wall geometries and contact-line dynamics.

preprint2021arXiv

A novel momentum-conserving, mass-momentum consistent method for interfacial flows involving large density contrasts

We propose a novel method for the direct numerical simulation of interfacial flows involving large density contrasts, using a Volume-of-Fluid method. We employ the conservative formulation of the incompressible Navier-Stokes equations for immiscible fluids in order to ensure consistency between the discrete transport of mass and momentum in both fluids. This strategy is implemented on a uniform 3D Cartesian grid with a staggered configuration of primitive variables, wherein a geometrical reconstruction based mass advection is carried out on a grid twice as fine as that for the momentum. The implementation is in the spirit of Rudman (1998) [41], coupled with the extension of the direction-split time integration scheme of Weymouth & Yue (2010) [46] to that of conservative momentum transport. The resulting numerical method ensures discrete consistency between the mass and momentum propagation, while simultaneously enforcing conservative numerical transport to arbitrary levels of precision in 3D. We present several quantitative comparisons with benchmarks from existing literature in order to establish the accuracy of the method, and henceforth demonstrate its stability and robustness in the context of a complex turbulent interfacial flow configuration involving a falling raindrop in air.

preprint2021arXiv

A numerical cough machine

We introduce a simplified model of physiological coughing or sneezing, in the form of a thin liquid layer subject to a rapid (30 m/s) air stream. The setup is simulated using the Volume-Of-Fluid method with octree mesh adaptation, the latter allowing grid sizes small enough to capture the Kolmogorov length scale. The results confirm the trend to an intermediate distribution between a Log-Normal and a Pareto distribution $P(d) \propto d^{-3.3}$ for the distribution of droplet sizes in agreement with a previous re-analysis of experimental results by one of the authors. The mechanism of atomisation does not differ qualitatively from the multiphase mixing layer experiments and simulations. No mechanism for a bimodal distribution, also sometimes observed, is evidenced in these simulations.

preprint2020arXiv

Dynamic wetting failure in curtain coating by the Volume-of-Fluid method

In this paper we investigate dynamic wetting in the curtain coating configuration. The two-phase Navier-Stokes equations are solved by a Volume-of-Fluid method on an adaptive Cartesian mesh. We introduce the Navier boundary condition to regularize the solution at the triple point and remove the implicit numerical slip induced by the cell-centered interface advection. We use a constant contact angle to describe the dynamic contact line. The resolution of the governing equations allows us to predict the substrate velocity at which wetting failure occurs. The model predictions are compared with prior computations and experimental observations.

preprint2020arXiv

Short-term oscillation and falling dynamics for a water drop dripping in quiescent air

The short-term transient falling dynamics of a dripping water drop in quiescent air has been investigated through both simulation and experiment. The focus is on the short term behavior and the time range considered covers about eight dominant second-mode oscillations of the drop after it is formed. Due to the small fluid inertia the growth of the drop is quasi-static and is well captured by the static pendant drop theory. Nevertheless, the pinching dynamics and the resulting post-formation state of the drop trigger a nonlinear oscillation when the drop falls. The initial shape of the drop when it is just formed is decomposed into spherical harmonic modes. The pinching dynamics such as interface overturning introduces small-scale variation on the drop contour, which in turn contributes to the finite amplitudes of the higher-order modes. Furthermore, the initial kinetic energy when the droplet is just formed is as important as the initial surface energy contained in the drop shape, and is found to amplify the initial oscillation amplitude and to induce a phase shift in the oscillation of all the modes. By incorporating both the initial surface and kinetic energy, the linear model for a free drop oscillation yields very good predictions for the second and third modes. The mode amplitude spectra show both the primary frequencies that are consistent with the Lamb's theory and the secondary frequencies arising from different modes due to nonlinear inter-mode coupling. The complex transient flow inside and outside the drop is induced by the interaction between the falling motion and the nonlinear oscillation. The streamlines indicate that the internal flow is substantially different from the Hill vortex for a falling drop without oscillation. The temporal evolutions of both the internal flow and the wake morphology follow the dominant second oscillation mode.

preprint2016arXiv

Vortices catapult droplets in atomization

A droplet ejection mechanism in planar two-phase mixing layers is examined. Any disturbance on the gas-liquid interface grows into a Kelvin-Helmholtz wave, and the wave crest forms a thin liquid film that flaps as the wave grows downstream. Increasing the gas speed, it is observed that the film breaks up into droplets which are eventually thrown into the gas stream at large angles. In a flow where most of the momentum is in the horizontal direction, it is surprising to observe these large ejection angles. Our experiments and simulations show that a recirculation region grows downstream of the wave and leads to vortex shedding similar to the wake of a backward-facing step. The ejection mechanism results from the interaction between the liquid film and the vortex shedding sequence: a recirculation zone appears in the wake of the wave and a liquid film emerges from the wave crest; the recirculation region detaches into a vortex and the gas flow over the wave momentarily reattaches due to the departure of the vortex; this reattached flow pushes the liquid film down; by now, a new recirculation vortex is being created in the wake of the wave - just where the liquid film is now located; the liquid film is blown up from below by the newly formed recirculation vortex in a manner similar to a bag-breakup event; the resulting droplets are catapulted by the recirculation vortex.

preprint2015arXiv

Droplet impact on a thin liquid film: anatomy of the splash

We investigate the dynamics of drop impact on a thin liquid film at short times in order to identify the mechanisms of splashing formation. Using numerical simulations and scaling analysis, we show that the splashing formation depends both on the inertial dynamics of the liquid and the cushioning of the gas. Two asymptotic regimes are identified, characterized by a new dimensionless number $J$: when the gas cushioning is weak, the jet is formed after a sequence of bubbles are entrapped and the jet speed is mostly selected by the Reynolds number of the impact. On the other hand, when the air cushioning is important, the lubrication of the gas beneath the drop and the liquid film controls the dynamics, leading to a single bubble entrapment and a weaker jet velocity.

preprint2014arXiv

Scaling properties of viscous fingering

We present a study of viscous fingering using the Volume Of Fluid method and a central injection geometry, assuming a Laplacian field and a simple surface tension law. As in experiments we see branched structures resulting from the Saffman-Taylor instability. We find that the area $A$ of a viscous-fingering cluster varies as a simple power law $A \sim L^α$ of its interface length $L$. Our results are compared to previously published simulations in which the viscosity of the invading fluid is vanishing. We find differences in exponent $α$ and in the appearance of detached droplets and bubbles.

preprint2009arXiv

Impact of a Viscous Liquid Drop

We simulate the impact of a viscous liquid drop onto a smooth dry solid surface. As in experiments, when ambient air effects are negligible, impact flattens the falling drop without producing a splash. The no-slip boundary condition at the wall produces a boundary layer inside the liquid. Later, the flattening surface of the drop traces out the boundary layer. As a result, the eventual shape of the drop is a "pancake" of uniform thickness except at the rim, where surface tension effects are significant. The thickness of the pancake is simply the height where the drop surface first collides with the boundary layer.