[MCrossover97]

Hi everyone,

I want to implement CSF formulation for modeling surface tension in my momentum equation and for some reason, I don't want to use built-in interFoam modules. One formulation of CSF can be found in this paper.

The general form of the momentum equation as well as the CSF formulation is attached. At first, I came up with something like this as Fc:

sigma * fvc::laplacian(fvc::grad(Sl)/mag(fvc::grad(Sl)),Sl)

But this doesn't work. The problem is that the laplacian doesn't accept fvc::grad(Sl) and I should modify it. Can you please help me implement it in the momentum equation? Can you also tell me if this form is even a correct form of implementing such a formula with respect to OpenFOAM programming standards?