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

Cutting cells for constructing alpha=0.5

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   May 20, 2020, 14:05
Post Cutting cells for constructing alpha=0.5
  #1
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 https://github.com/isoadvector 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.
forAll(mesh.cells(),ci)
{
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,0.5); ? What are the differences between them?


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

Old   May 23, 2020, 19:57
Default
  #2
New Member
 
Yağmur GÜLEÇ
Join Date: May 2014
Posts: 5
Rep Power: 8
yagmur_89 is on a distinguished road
I could manage to extract isocells at volume fraction 0.5. I was making mistake. Instead, I should have interpolated as follows:

volPointInterpolation vpi_(alpha_mesh);

scalarField ap_;
ap_ = vpi_.interpolate(alpha1);
isoCutCell icc(alpha_mesh,ap_);
volPointInterpolation vpi_(alpha_mesh);
yagmur_89 is offline   Reply With Quote

Reply

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
Cutting cells for constructing alpha=0.5 yagmur_89 OpenFOAM 0 May 20, 2020 13:46
cellZone not taking all the cells inside rahulksoni OpenFOAM Running, Solving & CFD 6 January 25, 2019 00:11
[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 06:45.