CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Programming & Development

Eliptic Relaxation in interFoam

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

Reply
 
LinkBack Thread Tools Display Modes
Old   June 1, 2010, 17:05
Default Eliptic Relaxation in interFoam
  #1
Senior Member
 
santiagomarquezd's Avatar
 
Santiago Marquez Damian
Join Date: Aug 2009
Location: Santa Fe, Santa Fe, Argentina
Posts: 418
Rep Power: 15
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:

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.
__________________
Santiago MÁRQUEZ DAMIÁN, Ph.D.
Post-doctoral Fellow
Research Center for Computational Mechanics (CIMEC) - CONICET/FICH-UNL
T.E.: 54-342-4511594 Ext. 1005
Güemes 3450 - (3000) Santa Fe
Santa Fe - Argentina
http://www.cimec.org.ar
santiagomarquezd is offline   Reply With Quote

Reply

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 08:38.