CFD Online Discussion Forums

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   OpenFOAM Programming & Development (http://www.cfd-online.com/Forums/openfoam-programming-development/)
-   -   Eliptic Relaxation in interFoam (http://www.cfd-online.com/Forums/openfoam-programming-development/76685-eliptic-relaxation-interfoam.html)

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.


All times are GMT -4. The time now is 00:18.