CFD Online Logo CFD Online URL
Home > Forums > OpenFOAM Programming & Development

Eliptic Relaxation in interFoam

Register Blogs Members List Search Today's Posts Mark Forums Read

LinkBack Thread Tools Display Modes
Old   June 1, 2010, 17:05
Default Eliptic Relaxation in interFoam
Senior Member
santiagomarquezd's Avatar
Santiago Marquez Damian
Join Date: Aug 2009
Location: Santa Fe, Santa Fe, Argentina
Posts: 438
Rep Power: 17
santiagomarquezd will become famous soon enough
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();
00104     // Cell gradient of alpha
00105     volVectorField gradAlpha = fvc::grad(alpha1_);
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());
00113     // Face unit interface normal
00114     surfaceVectorField nHatfv = gradAlphaf/(mag(gradAlphaf) + deltaN_);
00115     correctContactAngle(nHatfv.boundaryField());
00117     // Face unit interface normal flux
00118     nHatf_ = nHatfv & Sf;
here alpha1_ field is that was read from /0 directory without any smoothing proccess.
Research Scientist
Research Center for Computational Methods (CIMEC) - CONICET/UNL
Tel: 54-342-4511594 Int. 7032
Colectora Ruta Nac. 168 / Paraje El Pozo
(3000) Santa Fe - Argentina.
santiagomarquezd is offline   Reply With Quote


Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On

Similar Threads
Thread Thread Starter Forum Replies Last Post
InterFoam stops after deltaT goes to 1e14 francesco_b OpenFOAM Running, Solving & CFD 8 July 31, 2013 02:29
Purpose of relaxation factors Mohsin FLUENT 5 April 30, 2010 11:57
Open Channel Flow using InterFoam type solver sxhdhi OpenFOAM Running, Solving & CFD 3 May 5, 2009 21:58
Relaxation and convergence sammi Phoenics 0 March 20, 2008 04:32
Question on adjusting relaxation factor CFD Rookie Main CFD Forum 3 January 26, 2004 15:37

All times are GMT -4. The time now is 15:16.