Quantifying interface compression (interFoam)

mvdl1996 · February 10, 2020, 09:53

Hello all,

I am using interFoam for several multiphase simulations. In the alpha (volume fraction) equation, there is a constant cAlpha that has to do with interface compression. I noticed that for some cases, cAlpha has a large influence whereas for other cases it does not. I understand why, from the equations, but I would like to quanitfy / visualize this.

In order to do this, I would like to know the value of phic and phir which appear in the alphaEqn.H file (see equations from that file below).
phic = mag(phi/mesh.magSf());
phic = min(interface.cAlpha()*phic, max(phic));
surfaceScalarField phir = phic*interface.nHatf();

I tried using output stream and a swakExpression to accomplish this, but due to its surfaceScalar nature it seems quite hard as it is only defined for each cell surface and not just for each cell like volumeScalars.
Can anyone help me with suggestions how I can quantify or visualize phir and/or phic?

Thanks!

einstein_zee · April 27, 2023, 06:51

Hi there,

I just came across this post and it certainly is too late but just wanted to mention that the same issue was addressed in another thread "Visualisation of a surfaceScalarField with paraview". As Mattijs said you need to use Glyph...

February 10, 2020, 09:53	Quantifying interface compression (interFoam)	#1
mvdl1996 New Member Mike Join Date: Sep 2019 Posts: 2 Rep Power: 0	Hello all, I am using interFoam for several multiphase simulations. In the alpha (volume fraction) equation, there is a constant cAlpha that has to do with interface compression. I noticed that for some cases, cAlpha has a large influence whereas for other cases it does not. I understand why, from the equations, but I would like to quanitfy / visualize this. In order to do this, I would like to know the value of phic and phir which appear in the alphaEqn.H file (see equations from that file below). phic = mag(phi/mesh.magSf()); phic = min(interface.cAlpha()phic, max(phic)); surfaceScalarField phir = phicinterface.nHatf(); I tried using output stream and a swakExpression to accomplish this, but due to its surfaceScalar nature it seems quite hard as it is only defined for each cell surface and not just for each cell like volumeScalars. Can anyone help me with suggestions how I can quantify or visualize phir and/or phic? Thanks!

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
How to have sharp interface in nano scale two phase flow problem using interFoam	hosseinfathi	OpenFOAM Programming & Development	10	April 9, 2017 10:07
Error - Solar absorber - Solar Thermal Radiation	MichaelK	CFX	12	September 1, 2016 06:15
Waterwheel shaped turbine inside a pipe simulation problem	mshahed91	CFX	3	January 10, 2015 12:19
abnormal temperature near interface when adding energy equation to interFoam	houkensjtu	OpenFOAM	4	June 26, 2013 14:40
VOF interface compression and mesh size	idrama	OpenFOAM	1	September 26, 2010 07:36

April 27, 2023, 06:51		#2
einstein_zee Member Hosein Join Date: Nov 2011 Location: Germany Posts: 94 Rep Power: 15	Hi there, I just came across this post and it certainly is too late but just wanted to mention that the same issue was addressed in another thread "Visualisation of a surfaceScalarField with paraview". As Mattijs said you need to use Glyph...