CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Community Contributions

[isoAdvector] IsoAdvector: A new interface advection scheme for interFoam type calculations

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

Like Tree62Likes

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   July 2, 2019, 10:35
Default setAlphaField trouble with waveAlpha boundary condition
  #61
New Member
 
Marc Batlle
Join Date: Mar 2017
Posts: 4
Rep Power: 4
Marc Batlle is on a distinguished road
Hi Foamers,

Congratulations for your isoAdvector addition Mr. Roenby, which looks it is working really good in water waves environment, and furthermore, congrats on the adaptation for deformable meshes of this interface method.


Problem:
I am facing a problem these days when running the setAlphaField tool, which allows higher precision when defining the water depths, on OpenFOAM versions 1812 and 1906 but which were not appearing on version 1806. The issue appears when adding the "waveAlpha" boundary condition for alpha.water.

- I added the case from the tutorial (of-v1812) streamFunction with the Allrun switched to run setAlphaField instead of setFields

Question:
I am also really interested in the addition of isoAdvector in compressible two phase solvers as compressibleInterFoam. I will be testing the library from HenningScheufler these days, but I expected the realise on the official 1906 version. Did you found any issue in it you could prevent me from?

Thanks again for your dedication and contributions,

Marc
Attached Files
File Type: gz waveExampleStreamFunction.tar.gz (3.3 KB, 3 views)
Marc Batlle is offline   Reply With Quote

Old   July 2, 2019, 11:16
Default
  #62
Member
 
Johan Roenby
Join Date: May 2011
Location: Denmark
Posts: 77
Rep Power: 16
roenby will become famous soon enough
Hi Marc

setAlphaField issue:

setAlphaField complains about the boundary waveAlpha because it is not aware of it because the libwaveModels.so was not included in the Make/options file of the setAlphaField application (should be in the EXE_LIBS list). Good thing is you can simply add the following line to your case controlDict to make the setAlpahField aware of the waveAlpha BC:

libs ( "libwaveModels.so" );

VoFLibrary of Henning Scheufler:
This has now been added as a community contribution:
https://develop.openfoam.com/Communi...on-vof-library
(you might need to log in and/or request access to the project)
It will now be further tested and will hopefully make it into v1912.

Best regards,
Johan
roenby is offline   Reply With Quote

Old   July 16, 2019, 11:46
Default resolving spurrious currents using isoAdvector scheme
  #63
New Member
 
Navid
Join Date: Nov 2015
Posts: 15
Rep Power: 6
navidamin is on a distinguished road
Hi Professor Roenby,
Recently I have been struggling with spurrious currents problem for fluid-gas interfaces, which is produced by the surface tension force term in the momentum equation. I was advised that isoAdvector scheme would resolve the random fluctations that occure at the interface.


So I took the sloshingTank2D case from ESIOpenFOAM tutorials, put gravity and movement of the tank to zero and set a value for surface tension. Of course no spurrious currents are produced for the original case, but when I scaled the geometry to the factor of 0.001(2cmX2cm), the spurrious current on the interface appeared as you can see in the attached image.

So, I was wondering if you think there is anything wrong with my setup? And if not, is there any solutions to this problem?
Attached Images
File Type: jpeg sloshingTank.jpeg (17.1 KB, 32 views)
navidamin is offline   Reply With Quote

Old   November 15, 2019, 04:57
Default IsoAdvector
  #64
New Member
 
Deutschland
Join Date: Jun 2019
Posts: 13
Rep Power: 2
Arghavani is on a distinguished road
Hi everyone,
Does anyone work on the rotation of the slotted disk with IsoAdvector? if so, in which case of interisofoam can I find?
Arghavani is offline   Reply With Quote

Old   November 15, 2019, 05:04
Default
  #65
Member
 
Johan Roenby
Join Date: May 2011
Location: Denmark
Posts: 77
Rep Power: 16
roenby will become famous soon enough
Hi @Arghavani


Test case can be found here:

https://develop.openfoam.com/Develop...idBodyRotation


Youtube video:
https://youtu.be/p03yDzcdv6c


As always: Help the community by filing a bug report if you find something not working properly.


Cheers,
Johan
roenby is offline   Reply With Quote

Old   November 15, 2019, 05:52
Default IsoAdvector
  #66
New Member
 
Deutschland
Join Date: Jun 2019
Posts: 13
Rep Power: 2
Arghavani is on a distinguished road
Thanks a lot, Johan for your reply,

I will check as soon as possible.
I did Young's method for the slotted disk via Matlab, and the volume fraction breaks down so quickly. Do you know that this case in Isoadvector can support the interface capturing and reconstructing the interface?

kind regards,
Arghavan
Arghavani is offline   Reply With Quote

Old   November 15, 2019, 06:38
Default
  #67
Member
 
Johan Roenby
Join Date: May 2011
Location: Denmark
Posts: 77
Rep Power: 16
roenby will become famous soon enough
Yes it can. The linked video gives visual evidence. For alpha bounding and volume conservation details, run the case and look in the log file where this info is displayed
roenby is offline   Reply With Quote

Old   November 15, 2019, 06:45
Default
  #68
New Member
 
Deutschland
Join Date: Jun 2019
Posts: 13
Rep Power: 2
Arghavani is on a distinguished road
Thanks, Johan,

yes, you are right. I am going to do that.

Best,
Arghavan
Arghavani is offline   Reply With Quote

Old   November 15, 2019, 09:02
Default
  #69
New Member
 
Deutschland
Join Date: Jun 2019
Posts: 13
Rep Power: 2
Arghavani is on a distinguished road
Hey Johan,
I ran already the simulation as you told. But I got the error : ./Allrun: 7: ./Allrun: restore0Dir: not found
I download Isoadvector last weak in my tutorial folder but still is not working and in each log files I have these errors : ./Allrun: 93: ./Allrun: interIsoFoam: not found
FOAM FATAL IO ERROR:
keyword centre is undefined in dictionary "/home/izadshenas/notchedDiscInSolidBodyRotation/system/setAlphaFieldDict"

file: /home/izadshenas/notchedDiscInSolidBodyRotation/system/setAlphaFieldDict from line 18 to line 22.

From function const Foam::entry& Foam::dictionary::lookupEntry(const Foam::word&, bool, bool) const
in file db/dictionary/dictionary.C at line 566.

FOAM exiting


Do you think that these error s are related to the installation?



kind regards,
Arghavan
Arghavani is offline   Reply With Quote

Old   November 18, 2019, 08:30
Default IsoAdvector
  #70
New Member
 
Deutschland
Join Date: Jun 2019
Posts: 13
Rep Power: 2
Arghavani is on a distinguished road
Hi Johan,

for some other projects, I have OpenFoam of5x on my computer at work and until now I didn't need the Isoadvector and InterIsoFoam solvers but now I want to do another project which needs theses solvers but I am not sure that the of5x version of OpenFoam has InterisoFoam solver or not?

kind regards,
Arghavan
Arghavani is offline   Reply With Quote

Old   December 10, 2019, 07:22
Default
  #71
Member
 
Johan Roenby
Join Date: May 2011
Location: Denmark
Posts: 77
Rep Power: 16
roenby will become famous soon enough
IsoAdvector from github.com/isoAdvector works with OpenFOAM-5.x.

The isoAdvector solver in that repo is called interFlow - not interIsoFoam. If you have the Foundation version of OpenFOAM loaded (the one from openfoam.org) that would explain the error "./Allrun: 93: ./Allrun: interIsoFoam: not found" in your previous post.

interIsoFoam is the name of the corresponding solver using isoAdvector in the ESI-OpenCFD version of OpenFOAM (openfoam.com). That OpenFOAM version has a restore0Dir function to restore the 0 directory from the 0.orig directory. The Foundation version of OpenFOAM (openfoam.org) does not have that. Hence the first error in your previous post: "./Allrun: 7: ./Allrun: restore0Dir: not found".

Quote:
Originally Posted by Arghavani View Post
Hi Johan,

for some other projects, I have OpenFoam of5x on my computer at work and until now I didn't need the Isoadvector and InterIsoFoam solvers but now I want to do another project which needs theses solvers but I am not sure that the of5x version of OpenFoam has InterisoFoam solver or not?

kind regards,
Arghavan
roenby is offline   Reply With Quote

Reply

Tags
interface, interfoam, isoadvector, multiphase, unstructured mesh, vof

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
[Other] simulation of closing the gate using moving mesh simin_ds OpenFOAM Meshing & Mesh Conversion 8 April 12, 2019 06:49
rhoPimpleFoam hardship petrus OpenFOAM Running, Solving & CFD 0 October 7, 2016 03:41
Wrong flow in ratating domain problem Sanyo CFX 17 August 15, 2015 07:20
interFoam/kOmegaSST tank filling with printStackError/Mules simpomann OpenFOAM Running, Solving & CFD 3 February 17, 2014 18:06
T Junction Stability ignacio OpenFOAM Running, Solving & CFD 5 May 2, 2013 11:44


All times are GMT -4. The time now is 12:20.