CFD Online Logo CFD Online URL
Home > Forums > Software User Forums > OpenFOAM

Cutting cells for constructing alpha=0.5

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

LinkBack Thread Tools Search this Thread Display Modes
Old   May 20, 2020, 13:46
Default Cutting cells for constructing alpha=0.5
New Member
Yağmur GÜLEÇ
Join Date: May 2014
Posts: 5
Rep Power: 8
yagmur_89 is on a distinguished road
Dear Foamers,

I tried to cut the cells to find isoSurface at 0.5 volume fraction for a VOF application using the isoAdvector concept devised by

Roenby, J., Bredmose, H., & Jasak, H. (2016). A computational method for sharp interface advection. Royal Society open science, 3(11), 160405.

I compiled codes from for Foam extend 4.1 version. First, I constructed an icc object isoCutCell icc(mesh,alpha1) using the volume fraction field and mesh itself. Then, to compute which cells are cut by0.5 volume fraction.

label cellStatus = icc.calcSubCell(ci,0.5);

To calculate surface areas and to visualize the concerning variable I created volVectorField called temp to save the surface areas for cells that were cut by 0.5 volume fraction.
label cellStatus = icc.calcSubCell(ci,1);
//label cellStatus=icc.vofCutCell(ci,0.5,1e-3,100);

if (cellStatus == 0)

temp.internalField()[ci]= icc.isoFaceArea();


However, the problem is that the areas saved in temp field is in very irrelevant positions, as seen in Fig in the attachment. What could be the reason? Should I

use the code below

label cellStatus=icc.vofCutCell(ci,0.5,1e-3,100) instead of label cellStatus = icc.calcSubCell(ci,1); ? What are the differences between them?

Thanks in advance,
Best Regards,
Attached Images
File Type: jpg isoCutCells.jpg (30.8 KB, 6 views)
yagmur_89 is offline   Reply With Quote


interface contruction, isoadvector, isocutcell

Thread Tools Search this Thread
Search this Thread:

Advanced Search
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 Off
Pingbacks are On
Refbacks are On

Similar Threads
Thread Thread Starter Forum Replies Last Post
cellZone not taking all the cells inside rahulksoni OpenFOAM Running, Solving & CFD 6 January 25, 2019 00:11
[snappyHexMesh] Problem: after snappyHexMesh, the cells size are not the same kanes OpenFOAM Meshing & Mesh Conversion 0 January 25, 2016 08:06
[snappyHexMesh] No layers in a small gap bobburnquist OpenFOAM Meshing & Mesh Conversion 6 August 26, 2015 09:38
[snappyHexMesh] snappyHexMesh matches wrong cells to CellZone Siegunn OpenFOAM Meshing & Mesh Conversion 4 July 31, 2015 05:10
physical boundary error!! kris Siemens 2 August 3, 2005 00:32

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