santiagomarquezd |
June 1, 2010 17:05 |
Eliptic Relaxation in interFoam
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:
Code:
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.
|