CFD Online Logo CFD Online URL
Home > Forums > OpenFOAM

Signs used in advecting volume fractions in interfoam

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

LinkBack Thread Tools Display Modes
Old   June 29, 2017, 16:42
Default Signs used in advecting volume fractions in interfoam
New Member
Join Date: May 2017
Posts: 10
Rep Power: 2
Katt is on a distinguished road

I am aware of the standard equation used to advect volume fractions and about the derivation regarding "interface compression" that interfoam solver uses.

Coming to the code in "alphaEqn.H", I have the following two questions:
surfaceScalarField phic = mag(phi/mesh.magSf());                      
phic = min(interface.cAlpha()*phic, max(phic));                      
 surfaceScalarField phir = phic*interface.nHatf();                        
   for (int aCorr=0; aCorr<nAlphaCorr; aCorr++)                      
      surfaceScalarField phiAlpha =   
                                  fvc::flux ( phi,  
                                  alphaScheme  )
                               + fvc::flux ( -fvc::flux(-phir, scalar(1) - alpha1, alpharScheme),
                                              alpharScheme );
Q1) I don't really understand the sign convention used for the compression term.
In the relative flux term: -fvc::flux(-phir, ---). So, arent you finally adding an additional convective term to the first term of "phi*alpha1"?

Q2) If C-alpha =1 and considering the relative velocity acts as a shock front against the original convective term, wouldnt both the flux term and the relative flux term balance each other and cancel out leading to no advection?

I do realise my imagination is wrong as I do see correct flow behaviour. Can anyone enlighten me.

Katt is offline   Reply With Quote


advection, interface compression, interfoam

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
alphaEqn.H in twoPhaseEulerFoam cheng1988sjtu OpenFOAM Bugs 15 May 1, 2016 16:12
interFoam: No finite volume options present mehdinekooee OpenFOAM Running, Solving & CFD 4 July 15, 2014 11:05
InterFoam channel Flow runnig very slowly Wokl OpenFOAM Running, Solving & CFD 4 January 9, 2014 03:52
volume fractions in homogeneous flow model? Matthias Rohner CFX 1 August 26, 2007 18:42
Higher order for phase volume fractions? George Bergantz CFX 2 December 11, 2000 12:55

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