[kalaikdr]

Hi,

I am trying to implement the CLSVOF scheme in interFoam. as given to Albadawi et al, International Journal of Multiphase Flow, 2013.

In his paper, it is given that the redistancing of the level set function is solved by using the equation

d(phi)/dt' = sign(phi)*(1-|grad(phi)|)

where t' is an artificial time step chosen as delta t' = 0.1*delta x

when I tried it on interFoam as

solve
(
fvm::ddt(lsPhi)
+(sign(lsPhi) * mag(fvc::grad(lsPhi)))
-sign(lsPhi)
);

Foam compiles, but crashes on running due to dimension mismatch.

Any idea how to set it right? Does t' have the units of distance, as it is based on the mesh size?

Also, how do I solve this in a loop of artificial time until it converges to the steady state solution?