Abstract:We propose an alternative method to implement the contact angle boundary condition on a solid wall and apply it in hybrid lattice-Boltzmann finite-difference simulations of two-phase flows with immersed surfaces in which the flow equations are solved by the lattice-Boltzmann method and the interface equations are solved by the finite-difference method. Using the hyperbolic tangent profile of the order parameter across an interface in phase-field theory, we were able to obtain its unknown value at a ghost point from the information at only one point in the fluid domain. This is in contrast with existing approaches relying on interpolations involving several points. The special feature allows it to be more easily implemented on immersed surfaces cutting through the grid lines. It was properly incorporated into the framework of the hybrid lattice-Boltzmann finite-difference simulation, and applied successfully for several problems with different levels of complexity. First, the equilibrium shapes of a droplet on a sphere with different contact angles and radii were studied under cylindrical geometry and a good agreement with theoretical predictions was found. Preliminary studies on a three-dimensional droplet spreading on a sphere were also performed and the results agreed well with the corresponding axisymmetric results. Second, the spreading of a two-dimensional drop on an embedded inclined wall with a given contact angle was investigated and the results matched those on a flat wall at the domain boundary under the same condition. Third, capillary filling in a cylindrical tube was studied and the speed of the interface in the tube was found to follow Washburn's law. Fourth, a droplet impacting on a sphere was investigated and several different outcomes were captured depending on the Reynolds number, the viscosity ratio, and the wettability and radius of the sphere. Finally, the proposed method was shown to be capable of studying even more complicated problems involving the interaction between a droplet and multiple objects of different sizes and contact angles.

An Alternative Method to Implement Contact Angle Boundary Condition on Immersed Surfaces for Phase-Field Simulations

An Alternative Method to Implement Contact Angle Boundary Condition and Its Application in Hybrid Lattice-Boltzmann Finite-Difference Simulations of Two-Phase Flows with Immersed Surfaces

Simplified Method for Wetting on Curved Boundaries in Conservative Phase-Field Lattice-Boltzmann Simulation of Two-Phase Flows with Large Density Ratios

Contact-angle Implementation in Multiphase Smoothed Particle Hydrodynamics Simulations

An Energy-Based Equilibrium Contact Angle Boundary Condition on Jagged Surfaces for Phase-Field Methods.

A new approach for the implementation of contact line motion based on the phase-filed lattice Boltzmann method

An Inconsistency in the Formulation of the Phase-Field Method Between Flat and Curved Interfaces and Its Rational Remedy

Comparison of Contact Angle Models in Two‐Phase Flow Simulations Using a Conservative Phase Field Equation

Simplified wetting boundary scheme in phase-field lattice Boltzmann model for wetting phenomena on curved boundaries

A Diffuse-Interface Immersed-Boundary Method for Two-Dimensional Simulation of Flows with Moving Contact Lines on Curved Substrates

Implementation of Contact Angles in Pseudopotential Lattice Boltzmann Simulations with Curved Boundaries

Implementation of contact angles in the pseudopotential lattice Boltzmann simulations with curved boundaries

An Efficient Scheme for A Phase Field Model for the Moving Contact Line Problem with Variable Density and Viscosity

Implementation of contact line motion based on the phase-field lattice Boltzmann method

Applying Contact Angle to a 2D Multiphase Smoothed Particle Hydrodynamics Model

Wetting Boundary Schemes in Modified Phase-Field Lattice Boltzmann Method for Binary Fluids with Large Density Ratios

Numerical Simulation of Water Entry of Wedges Based on the Cip Method

Wetting Boundary Conditions in Numerical Simulation of Binary Fluids by Using Phase-Field Method: Some Comparative Studies and New Development

Implementing Contact Angle Hysteresis in Moving Mesh-Based Two-Phase Flow Numerical Simulations

Enhanced conservative phase field method for moving contact line problems