Abstract:The Immersed Boundary Method (IBM) is one of the popular one-fluid mixed Eulerian-Lagrangian methods to simulate motion of droplets. While the treatment of a moving complex boundary is an extremely time consuming and formidable task in a traditional boundary-fitted fluid solver, the one-fluid methods provide a relatively easier way to track moving interfaces on a fixed Cartesian grid since the regeneration of a mesh system that conforms to the interface at every time step can be avoided. In the IBM, a series of connected Lagrangian markers are used to represent a fluid-fluid interface and the boundary condition is enforced by adding a forcing term to the Navier-Stokes equations. It is known that the IBM suffers two problems. One is spontaneous generation of unphysical kinetic energy, which is known as parasitic currents, and the other is spurious reconstruction of interface. These two problems need to be solved for useful long-time-scale simulations of droplets with high density ratio and large surface tension. This work detects that the discrete delta function is the cause of unphysical parasitic currents. Specifically, the irrotational condition is not preserved when the common discrete delta function is used to spread the surface tension from Lagrangian markers to Cartesian grid cells. To solve this problem, a new scheme that preserves the irrotational condition is proposed to remove the spurious currents. Furthermore, for a smooth reconstruction of an interface, a B-spline fitting by least squares is adopted to relocate the Lagrangian markers. The conventional and new interpolation schemes are implemented in a multigrid finite volume Direct Numerical Simulation (DNS) solver and are subjected to standard test cases. It is confirmed that the unphysical parasitic currents are substantially reduced in the new scheme and droplet's surface fluctuation is eliminated.

A Discrete Droplet Method for Modelling Thin Film Flows

A Lagrangian meshfree model for solidification of liquid thin-films

Numerical Analysis of the Impact of Two Droplets with a Liquid Film Using an Incompressible SPH Method

A Front-Tracking Method for Simulating Interfacial Flows with Particles and Soluble Surfactants

Modelling two-dimensional droplet rebound off deep fluid baths

An implicit surface tension model for the analysis of droplet dynamics

Computing dynamics of thin films via large scale GPU-based simulations

A three-dimensional numerical model for the motion of liquid drops by the particle finite element method

The flow and solidification of a thin fluid film on an arbitrary three-dimensional surface

A smoothed dissipative particle dynamics methodology for wall-bounded domains

Droplet formation simulation using mixed finite elements

Adaptive mesh refinement for simulation of thin film flows

Mesoscopic Simulation of a Thinning Liquid Bridge Using the Dissipative Particle Dynamics Method.

Application of Dissipative Particle Dynamics to Fluid Flow Simulation

An Enhanced Semi-Implicit Particle Method for Simulating the Flow of Droplets with Free Surfaces

Modified Finite Difference Methods for Reynold Equation with Film Thickness Discontinuity

A one-dimensional mathematical model for shear-induced droplet formation in co-flowing fluids

New Immersed Boundary Method with Irrotational Discrete Delta Vector for Droplet Simulations with Large Density ratio

Application of the thin-film equations in modelling of Marangoni flow patterns amongst surfactant source and drain locations

Swimming droplet in 1D geometries, an active Bretherton problem

A simulation method for the wetting dynamics of liquid droplets on deformable membranes