[rcastilla]

Hi,

I have working with high weber number flows (capillary flows) and I am finding som important parasitic currents in the interface.

I have read that Brackbill (1992) J. Comput. Phys. 100 335-354, suggested to weight the surface tension force with the density in the cell in order to reduce this non-real currents.

I have modified this term in interFoam,

fvc::interpolate(interface.sigmaK())*fvc::snGrad(a lpha1)

multiplying by

fvc::interpolate(rho/(twoPhaseProperties.rho1()-twoPhaseProperties.r
ho2())

The parasitic currents have been appreciably reduced. I am surprised than I so simple thing has not been before implemented in interFoam. Maybe there is some drawback that I don't know?

Best regards