Source author record

Stéphane Popinet

Stéphane Popinet 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 cond-mat.soft

Catalog footprint

What is connected

8works

2topics

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.

preprint2022arXiv

High-resolution direct simulation of deep water breaking waves: transition to turbulence, bubbles and droplet production

[Abridged]We present high-resolution three-dimensional direct numerical simulations of breaking waves solving the two-phase Navier-Stokes equations. We investigate the role of the Reynolds and Bond numbers on the energy, bubble and droplet statistics of strong plunging breakers, and explore the asymptotic regimes at high Reynolds and Bond numbers to be compared with laboratory breaking waves. Energetically, the breaking wave transitions from laminar to three-dimensional turbulent flow on a timescale that depends on the turbulent Reynolds number up to a limiting value of $Re_λ\sim 100$, consistent with the mixing transition observed in other canonical turbulent flows. We characterize the role of capillary effects on the impacting jet and ingested main cavity shape and subsequent fragmentation process. We confirm two regimes in the bubble size distribution, separated by the Hinze scale $r_H$. We extend the buoyant-energetic scaling of Deike et al. (2016) to account for the cavity shape and its scale separation from the Hinze scale. We show resolved bubbles up to one order of magnitude below the Hinze scale and observe a good collapse of the numerical data compared to laboratory breaking waves (Deane and Stokes 2002). We resolve droplet statistics at high Bond number and our data show good agreement with recent experiments (Erinin et al., 2009) in various statistics. We discuss velocity distributions for the droplets, finding ejection velocities up to four times the phase speed of the wave, which are produced during the most intense splashing events of the breaking process.

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

Droplet migration in a Hele-Shaw cell: Effect of the lubrication film on the droplet dynamics

Droplet migration in a Hele--Shaw cell is a fundamental multiphase flow problem which is crucial for many microfluidics applications. We focus on the regime at low capillary number and three-dimensional direct numerical simulations are performed to investigate the problem. In order to reduce the computational cost, an adaptive mesh is employed and high mesh resolution is only used near the interface. Parametric studies are performed on the droplet horizontal radius and the capillary number. For droplets with an horizontal radius larger than half the channel height the droplet overfills the channel and exhibits a pancake shape. A lubrication film is formed between the droplet and the wall and particular attention is paid to the effect of the lubrication film on the droplet velocity. The computed velocity of the pancake droplet is shown to be lower than the average inflow velocity, which is in agreement with experimental measurements. The numerical results show that both the strong shear induced by the lubrication film and the three-dimensional flow structure contribute to the low mobility of the droplet. In this low-migration-velocity scenario the interfacial flow in the droplet reference frame moves toward the rear on the top and reverses direction moving to the front from the two side edges. The velocity of the pancake droplet and the thickness of the lubrication film are observed to decrease with capillary number. The droplet velocity and its dependence on capillary number cannot be captured by the classic Hele--Shaw equations, since the depth-averaged approximation neglects the effect of the lubrication film.

preprint2015arXiv

Droplet in micro-channels: A numerical approach using an adaptive two phase flow solver

We propose a numerical approach to study the mechanics of a flowing bubble in a constraint micro channel. Using an open source two phase flow solver (Gerris, gfs.sourceforge.net) we compute solutions of the bubble dynamics (i.e. shape and terminal velocity) induced by the interaction between the bubble movement, the Laplace pressure variation, and the lubrication film near the channel wall. Quantitative and qualitative results are presented and compared against both theory and experimental data for small Capillary numbers. We discuss the technical issues of explicit integration methods on small Capillary numbers computations, and the possibility of adding Van der Walls forces to give a more precise picture of the Droplet-based microfluidic problem.

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.

preprint2012arXiv

The granular silo as a continuum plastic flow: the hour-glass vs the clepsydra

The granular silo is one of the many interesting illustrations of the thixotropic property of granular matter: a rapid flow develops at the outlet, propagating upwards through a dense shear flow while material at the bottom corners of the container remains static. For large enough outlets, the discharge flow is continuous; however, by contrast with the clepsydra for which the flow velocity depends on the height of fluid left in the container, the discharge rate of granular silos is constant. Implementing a plastic rheology in a 2D Navier-Stokes solver (following the mu(I)-rheology or a constant friction), we simulate the continuum counterpart of the granular silo. Doing so, we obtain a constant flow rate during the discharge and recover the Beverloo scaling independently of the initial filling height of the silo. We show that lowering the value of the coefficient of friction leads to a transition toward a different behavior, similar to that of a viscous fluid, and where the filling height becomes active in the discharge process. The pressure field shows that large enough values of the coefficient of friction ($\simeq$ 0.3) allow for a low-pressure cavity to form above the outlet, and can thus explain the Beverloo scaling. In conclusion, the difference between the discharge of a hourglass and a clepsydra seems to reside in the existence or not of a plastic yield stress.