[Arif14]

Dear CFD Community,


I am attempting to reproduce the results presented in an article (please find attached a PDF file containing some information and screenshots). In brief, the authors used a lattice Boltzmann method to model a liquid plug rupture in a channel. Periodic boundary conditions were applied, and a pressure drop allowed the liquid plug to move until rupture.

I am trying to replicate these results in OpenFOAM using the VOF method instead (interFoam). I created the same geometry, applied periodic boundary conditions using cyclic patches, and implemented a fixed jump to simulate the pressure drop.
However, I am struggling to achieve the same results, specifically the rupture time of the plug. Despite double-checking the dimensions, material parameters (density, viscosity, surface tension), and refining the mesh, rupture has not occurred. Instead, the liquid plug remains wider compared to the results presented in the article.
As a CFD beginner, I am seeking guidance and advice from this community:

1. Could the discrepancy be due to differences in the methods (VOF vs. lattice Boltzmann)? From my understanding, VOF is a macroscopic method and may face challenges with interface tracking, especially during topology changes like rupture.
2. Alternatively, is this problem "simple" enough that the issue could lie in my specific model implementation in OpenFOAM?

Thank you very much for your time and support! Any advice or direction would be greatly appreciated.

Sincerely,
Arif