[santiagomarquezd]

Hi FOAMers!

I'm working in interFoam again, this time with some more experience... I was reading Rusche's Thesis and if I'm not wrong, Smoothing of alpha field before nhatf calculation, is no longer done? Or never was used really in the code?

From interfaceProperties.C (1.6) we have:

00102     const surfaceVectorField& Sf = mesh.Sf();
00103 
00104     // Cell gradient of alpha
00105     volVectorField gradAlpha = fvc::grad(alpha1_);
00106 
00107     // Interpolated face-gradient of alpha
00108     surfaceVectorField gradAlphaf = fvc::interpolate(gradAlpha);
00109     //gradAlphaf -=
00110     //    (mesh.Sf()/mesh.magSf())
00111     //   *(fvc::snGrad(alpha1_) - (mesh.Sf() & gradAlphaf)/mesh.magSf());
00112 
00113     // Face unit interface normal
00114     surfaceVectorField nHatfv = gradAlphaf/(mag(gradAlphaf) + deltaN_);
00115     correctContactAngle(nHatfv.boundaryField());
00116 
00117     // Face unit interface normal flux
00118     nHatf_ = nHatfv & Sf;

here alpha1_ field is that was read from /0 directory without any smoothing proccess.