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

PLIC in InterFOAM???

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

Like Tree5Likes
  • 1 Post By tomislav_maric
  • 2 Post By tomislav_maric
  • 1 Post By popinet
  • 1 Post By wavefunction

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   June 20, 2011, 18:01
Default PLIC in InterFOAM???
  #1
New Member
 
lakshay
Join Date: Apr 2010
Posts: 4
Rep Power: 15
lakshay26 is on a distinguished road
Hi,

Does InterFOAM like Fluent also uses PLIC. also which part of the solver code can confirm this. Need urgent help.

Thanks
lakshay26 is offline   Reply With Quote

Old   June 21, 2011, 04:13
Default
  #2
Senior Member
 
Suresh kumar Kannan
Join Date: Mar 2009
Location: Luxembourg, Luxembourg, Luxembourg
Posts: 129
Rep Power: 17
kumar is on a distinguished road
No interFOAM does not use PLIC like Fluent. OpenFOAM uses a scheme called interFAce compression scheme. You can look into

~/OpenFOAM/OpenFOAM-1.6.x/applications/solvers/multiphase/interFoam for the details and then for further details of compression scheme look into MULES . Search for it in the OpenFOAM directories. It will give you the details.

regards
K.Suresh kumar
kumar is offline   Reply With Quote

Old   September 15, 2011, 06:40
Default
  #3
Senior Member
 
Robert Castilla
Join Date: Apr 2009
Location: Spain
Posts: 108
Rep Power: 16
rcastilla is on a distinguished road
Im also interested in the implementation of PLIC in interFoam. Is there any work done in this direction? The interface compression procedure gives bad results form small domains.

Regards

Robert
rcastilla is offline   Reply With Quote

Old   September 15, 2011, 07:22
Default
  #4
Senior Member
 
Andrea Ferrari
Join Date: Dec 2010
Posts: 319
Rep Power: 16
Andrea_85 is on a distinguished road
Quote:
Originally Posted by rcastilla View Post
Im also interested in the implementation of PLIC in interFoam. Is there any work done in this direction? The interface compression procedure gives bad results form small domains.

Regards

Robert
I am also interested in some algorithm that makes the interface sharp. From what i known PLIC algorithm is normally used on very simple mesh (all hexas) and may be difficult to adapt to a completely tet mesh. Would be interisting if some "expert" tell us if it is possible or not the implementation in openFoam.
In particular what i found from my experience using interFoam are very high spurious velocities at the interface when the domain size becomes smaller and this tends to affect the results. For this type of problem, considerable progress has been made in recent times but, in order to implement one of these techniques are necessary geometrical information about the interface, which are not available in interFoam.
So would be really helpful knowing in which direction would be better to move.

Best

andrea
Andrea_85 is offline   Reply With Quote

Old   September 16, 2011, 05:11
Default
  #5
Senior Member
 
Tomislav Maric
Join Date: Mar 2009
Location: Darmstadt, Germany
Posts: 284
Blog Entries: 5
Rep Power: 21
tomislav_maric is on a distinguished road
Hi everyone,

I'm currently working on the "Geometric" Volume of Fluid solver in OpenFOAM (PLIC reconstruction + unsplit geometrical advection) as a part of my PhD topic. It's going quite well so far and the code will be released into the public domain as soon as it is possible (when it runs, validated, in parallel).

T.
Mahmoud_aboukhedr likes this.
tomislav_maric is offline   Reply With Quote

Old   September 16, 2011, 05:32
Default
  #6
Senior Member
 
Tomislav Maric
Join Date: Mar 2009
Location: Darmstadt, Germany
Posts: 284
Blog Entries: 5
Rep Power: 21
tomislav_maric is on a distinguished road
Quote:
Originally Posted by Andrea_85 View Post
I am also interested in some algorithm that makes the interface sharp. From what i known PLIC algorithm is normally used on very simple mesh (all hexas) and may be difficult to adapt to a completely tet mesh. Would be interisting if some "expert" tell us if it is possible or not the implementation in openFoam.
In particular what i found from my experience using interFoam are very high spurious velocities at the interface when the domain size becomes smaller and this tends to affect the results. For this type of problem, considerable progress has been made in recent times but, in order to implement one of these techniques are necessary geometrical information about the interface, which are not available in interFoam.
So would be really helpful knowing in which direction would be better to move.

Best

andrea
The problem with the true interface capturing method, if it is based on purely geometrical solution to the volume fraction equation, is the fact that all of the developments so far have been done on Carthesian equidistant meshes. I'm trying my best to develop the method using OpenFOAM with support for arbitrary unstructured mesh. It is, well, difficulut, for various reasons.

Adaptive Mesh Refinement has been coupled with the geometrical VoF method in an open source flow solver called Gerris which has remarkable results in this area, but cannot deal with complex geometries because of the nature of the mesh it is built upon. If you need an outstanding DNS solutions for jet breakup, wave breaking, various kinds of interface breakup and coalescence dealt with on very low spatial scales, but for simple geometries, this is the way to go right now. If I manage to implement a true geometrical VoF on top of the polyhedral mesh in OF, then, of course, please do use my code in stead.

If you are interested in killing spurious currents, check out the articles on balanced force algorithm from Francois et al, Renardy and Rendardy and, the best one in my view, Popinet et al. Popinet shows how the combination of the height functions, octree based local AMR, and properly numerically balanced contributions of the pressure gradient and the curvature balanced force), reduce spurious currents to machine tolerance. All of the advancements rely, of course, on the true geometrical captured interface (not numerically compressed + iso-surface).

T.
Julien and Pirlu like this.
tomislav_maric is offline   Reply With Quote

Old   September 16, 2011, 05:55
Default
  #7
Senior Member
 
Robert Castilla
Join Date: Apr 2009
Location: Spain
Posts: 108
Rep Power: 16
rcastilla is on a distinguished road
Hi, Tomislav,

do you know this paper?

X. Yang and James, A.J. "Analytic relations for reconstructing piecewise linear interfaces in triangular and tetrahedral grids" Journal of Comp. Physics, 214 (2006) 41-54

If not, I hope you will find it useful.

regards

Robert
rcastilla is offline   Reply With Quote

Old   September 16, 2011, 06:05
Default
  #8
Senior Member
 
Tomislav Maric
Join Date: Mar 2009
Location: Darmstadt, Germany
Posts: 284
Blog Entries: 5
Rep Power: 21
tomislav_maric is on a distinguished road
Quote:
Originally Posted by rcastilla View Post
Hi, Tomislav,

do you know this paper?

X. Yang and James, A.J. "Analytic relations for reconstructing piecewise linear interfaces in triangular and tetrahedral grids" Journal of Comp. Physics, 214 (2006) 41-54

If not, I hope you will find it useful.

regards

Robert
Hi Robert,

yep, I know that, and I did find it useful! Thanks a lot, anyway! Are you working on Geometrical VoF method as well?

T.
tomislav_maric is offline   Reply With Quote

Old   September 16, 2011, 06:45
Default
  #9
Senior Member
 
Andrea Ferrari
Join Date: Dec 2010
Posts: 319
Rep Power: 16
Andrea_85 is on a distinguished road
Quote:
Originally Posted by tomislav_maric View Post
The problem with the true interface capturing method, if it is based on purely geometrical solution to the volume fraction equation, is the fact that all of the developments so far have been done on Carthesian equidistant meshes. I'm trying my best to develop the method using OpenFOAM with support for arbitrary unstructured mesh. It is, well, difficulut, for various reasons.

Adaptive Mesh Refinement has been coupled with the geometrical VoF method in an open source flow solver called Gerris which has remarkable results in this area, but cannot deal with complex geometries because of the nature of the mesh it is built upon. If you need an outstanding DNS solutions for jet breakup, wave breaking, various kinds of interface breakup and coalescence dealt with on very low spatial scales, but for simple geometries, this is the way to go right now. If I manage to implement a true geometrical VoF on top of the polyhedral mesh in OF, then, of course, please do use my code in stead.

If you are interested in killing spurious currents, check out the articles on balanced force algorithm from Francois et al, Renardy and Rendardy and, the best one in my view, Popinet et al. Popinet shows how the combination of the height functions, octree based local AMR, and properly numerically balanced contributions of the pressure gradient and the curvature balanced force), reduce spurious currents to machine tolerance. All of the advancements rely, of course, on the true geometrical captured interface (not numerically compressed + iso-surface).

T.

This is exactly the point!. As i said in my previous post, in my opinion all these techniques work well only with a sharp interface formulation. Im very familiar with the papers you mentioned and i've already tried to adapt the balance force algorithm to OpenFoam but without success in reducing spurious current, due to the fact that my interface is diffused over some cells.
I'm familiar with Gerris as well, but as you said it works only for very simple geometry and mesh, which is not my case (you cannot include obstacle and treat the contact angle as usual for example and you can not use polyhedral mesh). Currently i'm trying to "live" with the problem of spurious velocities and i would like to test how much they affect the results when the dimension of the domain decreases (for simple case and compared to other codes). The major problem is related to the time step, which is controlled by the CFL number and so by the highest velocity in the domain.
It would be great if you want to share your results in case you will be successful!!

Best

andrea
Andrea_85 is offline   Reply With Quote

Old   September 16, 2011, 07:01
Default
  #10
Senior Member
 
Tomislav Maric
Join Date: Mar 2009
Location: Darmstadt, Germany
Posts: 284
Blog Entries: 5
Rep Power: 21
tomislav_maric is on a distinguished road
Quote:
Originally Posted by Andrea_85 View Post
This is exactly the point!. As i said in my previous post, in my opinion all these techniques work well only with a sharp interface formulation. Im very familiar with the papers you mentioned and i've already tried to adapt the balance force algorithm to OpenFoam but without success in reducing spurious current, due to the fact that my interface is diffused over some cells.
I'm familiar with Gerris as well, but as you said it works only for very simple geometry and mesh, which is not my case (you cannot include obstacle and treat the contact angle as usual for example and you can not use polyhedral mesh). Currently i'm trying to "live" with the problem of spurious velocities and i would like to test how much they affect the results when the dimension of the domain decreases (for simple case and compared to other codes). The major problem is related to the time step, which is controlled by the CFL number and so by the highest velocity in the domain.
It would be great if you want to share your results in case you will be successful!!

Best

andrea
Right now I'm concentrated on coding the steps for the volume fraction equation solution, and as soon as the code runs (both interface steps), I'll need to validate it. I'll share the first results that show the code working properly when the time comes for it: there shall be debugging.

Best regards,
Tomislav
tomislav_maric is offline   Reply With Quote

Old   September 19, 2011, 10:58
Default
  #11
ngj
Senior Member
 
Niels Gjoel Jacobsen
Join Date: Mar 2009
Location: Copenhagen, Denmark
Posts: 1,900
Rep Power: 37
ngj will become famous soon enoughngj will become famous soon enough
Hi Tomislav

Congratulations on obtaining your master and also congratulations with obtaining a PhD scholarship. I am glad to see that you will continue to pursuit for a general method for interface tracking.

I need to begin this week to make a hexahedral approach similar to the one we talked about in Gothenburg last year, however, I did in my initial exercises last year experience problems with stability of the momentum equation at the interface. Hence, I would appreciate if you could either elaborate a little more on the references (full author list and title) or potentially send a copy to me.

All of my best wishes,

Niels
ngj is offline   Reply With Quote

Old   September 19, 2011, 17:14
Default
  #12
Senior Member
 
Tomislav Maric
Join Date: Mar 2009
Location: Darmstadt, Germany
Posts: 284
Blog Entries: 5
Rep Power: 21
tomislav_maric is on a distinguished road
Hi Niels,

thanks for the congrats! The master degree I will remember for as long as I live: the acknowledgement process and all that surrounded it was very.... special.

Just send me an e-mail at the CSI email.

Best,
Tomislav

Quote:
Originally Posted by ngj View Post
Hi Tomislav

Congratulations on obtaining your master and also congratulations with obtaining a PhD scholarship. I am glad to see that you will continue to pursuit for a general method for interface tracking.

I need to begin this week to make a hexahedral approach similar to the one we talked about in Gothenburg last year, however, I did in my initial exercises last year experience problems with stability of the momentum equation at the interface. Hence, I would appreciate if you could either elaborate a little more on the references (full author list and title) or potentially send a copy to me.

All of my best wishes,

Niels

Last edited by tomislav_maric; September 22, 2011 at 03:48. Reason: Editing the email...
tomislav_maric is offline   Reply With Quote

Old   October 12, 2011, 05:16
Default
  #13
New Member
 
Stephane Popinet
Join Date: Oct 2011
Posts: 1
Rep Power: 0
popinet is on a distinguished road
Hi Andrea et al,

I just stumbled upon this interesting discussion by chance.

Quote:
Originally Posted by Andrea_85 View Post
I'm familiar with Gerris as well, but as you said it works only for very simple geometry and mesh, which is not my case (you cannot include obstacle and treat the contact angle as usual for example and you can not use polyhedral mesh).
I don't really agree with this. Gerris can be used with complex geometries. With "complex" ranging from relatively simple such as this example:

http://gerris.dalembert.upmc.fr/gerr...ples/ship.html

to multiphase flows in actual injector geometries. For example Figure 22 of the recent paper by Dos Santos and Le Moyne:

http://ogst.ifpenergiesnouvelles.fr/...ogst100079.pdf

It is correct that at the moment contact angle conditions are not treated accurately at the "triple points" between interfaces and complex solid boundaries. On the other hand I would be cautious about methods claiming to solve this accurately while still having substantial problems with spurious currents around stationary droplets.

If what you meant is that Gerris does not (and cannot) use an unstructured polyhedral mesh, this is entirely correct, however this does not mean that problems in complex geometries cannot be solved (or that the mesh itself is necessarily simple).

That said, I think that implementing a geometric VOF model within OpenFOAM is a very interesting project.

cheers

Stephane
kk415 likes this.
popinet is offline   Reply With Quote

Old   October 12, 2011, 05:51
Default
  #14
Senior Member
 
Tomislav Maric
Join Date: Mar 2009
Location: Darmstadt, Germany
Posts: 284
Blog Entries: 5
Rep Power: 21
tomislav_maric is on a distinguished road
Quote:
Originally Posted by popinet View Post

I don't really agree with this. Gerris can be used with complex geometries. With "complex" ranging from relatively simple such as this example:

http://gerris.dalembert.upmc.fr/gerr...ples/ship.html

to multiphase flows in actual injector geometries. For example Figure 22 of the recent paper by Dos Santos and Le Moyne:

http://ogst.ifpenergiesnouvelles.fr/...ogst100079.pdf

(snip)

If what you meant is that Gerris does not (and cannot) use an unstructured polyhedral mesh, this is entirely correct, however this does not mean that problems in complex geometries cannot be solved (or that the mesh itself is necessarily simple).

(snip)
Stephane
Hi Stephane,

there must have been a misunderstanding, please accept my apologies for that. I was under the impression that dealing with complex geometries is a work in progress in Gerris.

Next time, I'll double check what I think I know (sorry again):

"Entirely automatic mesh generation in complex geometries"

http://gfs.sourceforge.net/wiki/inde...&printable=yes


Best regards,
Tomislav
tomislav_maric is offline   Reply With Quote

Old   October 12, 2011, 08:20
Default
  #15
Senior Member
 
Andrea Ferrari
Join Date: Dec 2010
Posts: 319
Rep Power: 16
Andrea_85 is on a distinguished road
Hi Stephan,
I apologize too for that, what i was saying was related to the particular case i want to simulate, for which i'm quite sure i cannot use Gerris (at least at the moment). I did not mean any kind of complex geometry or complex mesh, but it was more related to what I'm trying to do. So i'm sorry again to have explained myself bad.


Best


Andrea
Andrea_85 is offline   Reply With Quote

Old   May 12, 2015, 19:13
Default
  #16
Member
 
james wilson
Join Date: Aug 2014
Location: Orlando, Fl
Posts: 39
Rep Power: 11
jameswilson620 is on a distinguished road
Any progress with your Ph.D., Tomislav?

I am interested in interface reconstruction in Cartesian, structured grids and some cylindrical grids, all using hex cells.

I would like to accurately calculate the interface normal in a manner similar to that described by Kunkelmann: http://tuprints.ulb.tu-darmstadt.de/...Kunkelmann.pdf (pg.25).

James
jameswilson620 is offline   Reply With Quote

Old   January 15, 2016, 14:21
Default
  #17
Senior Member
 
Saideep
Join Date: Apr 2015
Location: INDIA
Posts: 203
Rep Power: 11
Saideep is on a distinguished road
Hi everyone,

I would like to revoke this thread. I am presently facing problems with badly smeared out interface using interfoam.

I tried to refine the domain fine but the smearing seems to reduce but not really satisfactory. Is this caused only due to spurious currents? {else in other words if I try to reduce spurious currents, would the smearing reduce?}

Any new updates over this issue Foamers.

Thanks and happy weekend,
-Saideep
Saideep is offline   Reply With Quote

Old   January 16, 2016, 16:13
Default
  #18
Senior Member
 
kmooney's Avatar
 
Kyle Mooney
Join Date: Jul 2009
Location: San Francisco, CA USA
Posts: 323
Rep Power: 17
kmooney is on a distinguished road
Quote:
Originally Posted by Saideep View Post
Hi everyone,

I would like to revoke this thread. I am presently facing problems with badly smeared out interface using interfoam.

I tried to refine the domain fine but the smearing seems to reduce but not really satisfactory. Is this caused only due to spurious currents? {else in other words if I try to reduce spurious currents, would the smearing reduce?}

Any new updates over this issue Foamers.

Thanks and happy weekend,
-Saideep
In general the OpenFOAM VOF solver just needs lots and lots of mesh to stay resolved. In general if you turn up the interface compression coeffecient you'll get less smearing, however, you end up increasing the magnitude of spurious currents (from what I've seen).

One effective way (I've tested it its awesome) to mitigate spurious currents is to implement this laplacian smoothing sweep:
http://scitation.aip.org/content/aip...1063/1.4756069

http://www.cfd-online.com/Forums/ope...bill-work.html

Cheers!
Kyle
kmooney is offline   Reply With Quote

Old   June 5, 2018, 07:27
Default Sharp interface
  #19
New Member
 
Saicharan
Join Date: Jan 2018
Location: Bangalore, India
Posts: 29
Rep Power: 8
wavefunction is on a distinguished road
Quote:
Originally Posted by kmooney View Post
In general the OpenFOAM VOF solver just needs lots and lots of mesh to stay resolved. In general if you turn up the interface compression coeffecient you'll get less smearing, however, you end up increasing the magnitude of spurious currents (from what I've seen).
Hello! Increasing cAlpha and refining the mesh are the only ways to get a sharper interface then, is it? Is there any other way to get a sharp interface? As I am simulating droplet evaporation, I require a sharp interface for my purpose.

Also, thank you for suggesting the laplacian smoother. I implemented it and it really sped up my simulation.
Elham likes this.
wavefunction is offline   Reply With Quote

Reply

Tags
interfoam, plic

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
Getting a concentration field around a bubble in InterFoam azman OpenFOAM Running, Solving & CFD 3 June 7, 2022 05:21
InterFoam stops after deltaT goes to 1e14 francesco_b OpenFOAM Running, Solving & CFD 9 July 25, 2020 07:36
Slow interFoam compared with other CFD tools? Ralph M OpenFOAM Programming & Development 1 November 17, 2010 07:46
Open Channel Flow using InterFoam type solver sxhdhi OpenFOAM Running, Solving & CFD 3 May 5, 2009 22:58
InterFoam PLIC ntrask OpenFOAM Running, Solving & CFD 1 January 8, 2009 13:32


All times are GMT -4. The time now is 06:39.