CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)
-   OpenFOAM Programming & Development (https://www.cfd-online.com/Forums/openfoam-programming-development/)
-   -   Eliptic Relaxation in interFoam (https://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 09:55.