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

State of the art concerning "spurious currents" in interFOAM

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

Like Tree5Likes
  • 2 Post By SaufiSangi
  • 1 Post By liquidspoon
  • 1 Post By liquidspoon
  • 1 Post By eric

Reply
 
LinkBack Thread Tools Display Modes
Old   October 3, 2017, 06:36
Default State of the art concerning "spurious currents" in interFOAM
  #1
New Member
 
Abd Essamade Saufi
Join Date: Apr 2017
Posts: 14
Rep Power: 2
SaufiSangi is on a distinguished road
Good morning everybody,
I am working on the simulation of droplets evaporation using the VOF method, modifying the already existing "interFOAM" solver. I started in absence of gravity, thus no surface tension was necessary in my calculations.

Now I am trying to activate gravity. My droplets have a diameter of 0.5-1 mm, and surface tension effects are definitely dominating. As you (who work on surface tension dominating flows in OpenFOAM) know, spurious currents arise and destroy the whole droplet after few time steps.

Now, I read ALL the post and threads concerning this problem in interFOAM, from 2009 to 2016. I found very interesting ideas, which I tried to implement, but none of them is really satisfactory for such small scales. The best one is the one proposed by Ali Q. Raeini et. al. [2012], but it still gives me some problems.

All of the "definitely solved" problems in VOF method concerning parasitic currents are based on the reconstruction of the interface (e.g VOF-PLIC), which I think it's what Fluent is using, or using marker functions ecc.
From what I understood these methods are quite difficult to implement in OpenFOAM, which already uses a compressive scheme for the interface reconstruction. I would be really scared if I had to implement geometrical methods in OpenFOAM, it does not seem to be easy at all, even for skilled programmers.

So, what's the situation now in 2017 about spurious currents in interFOAM? Are there any improvements/news about them? Is there anybody currently working on this major problem of interFOAM? As many people here, my work is not related to this specific problem, but I really need to solve them in order to proceed with my research, that's why I'm asking for your help

Thanks guys for you attention, any suggestion would be highly appreciated!!
Sangi
liquidspoon and sungda like this.
SaufiSangi is offline   Reply With Quote

Old   October 24, 2017, 08:18
Default
  #2
Member
 
Alex
Join Date: Jun 2011
Posts: 32
Rep Power: 8
liquidspoon is on a distinguished road
Hi Sangi,

I find myself in the same boat, and would be very happy to know if there are any robust solutions to this issue.

As a note, I have tried implementing the Raeini surface tension model (not the whole solver, just replacing the surface tension force calculation). It seemed to help reduce spurious currents, but it also seemed to damp out some real interface dynamics.

I have been meaning to try out the isoAdvectorFoam solver on some test cases. It was developed to address other issues, but may help out with the surface tension force calculations as a side effect (no proof of this, just an idea).

I think one part of the issue is the diffuse interface in interFoam. The surface tension forces are distributed across the interface weighted by grad(alpha) (basically a smoothed heaviside function). As a result, a large component of the force is applied to the low-inertia gas side. Typically you see the highest velocities on the gas side. This is still not the complete picture, because large currents can also appear on the liquid side if the surface tension effects are large enough.
SaufiSangi likes this.
liquidspoon is offline   Reply With Quote

Old   October 24, 2017, 15:46
Default
  #3
New Member
 
Abd Essamade Saufi
Join Date: Apr 2017
Posts: 14
Rep Power: 2
SaufiSangi is on a distinguished road
Hi Alex, thank you for your answer
What are you working on exactly?

From what I know isoAdvector does not solve the problem of spurious currents, there is a thread about that (don't remember where though).
Yes, the highest velocities are in the gas-side, but even the density correction proposed by Brackbill doesn't seem to be enough. This is a major problem for me, I am working with 1 mm liquid droplets, which are extremely difficult to stabilize. I need a smart way to rapidly solve this problem, since this is not the main topic of my work, I don't want to loose too much time on this.
So far I implemented the Raeini model and I am trying to couple it with a Level-Set method, this can be a way.

Keep in touch!

Sangi
SaufiSangi is offline   Reply With Quote

Old   October 24, 2017, 15:51
Default
  #4
New Member
 
Abd Essamade Saufi
Join Date: Apr 2017
Posts: 14
Rep Power: 2
SaufiSangi is on a distinguished road
Here it is by the way: the thread about isoAdvector

IsoAdvector release
SaufiSangi is offline   Reply With Quote

Old   October 24, 2017, 15:58
Default
  #5
Member
 
Alex
Join Date: Jun 2011
Posts: 32
Rep Power: 8
liquidspoon is on a distinguished road
Hi Sangi,
I am also working on capillary scale flows (small droplets, bubbles/Taylor flows in minichannels).

I tried out this coupled LS-VOF solver a while back: https://bitbucket.org/nunuma/public . It was very nice for the developer to release it, but I had mixed experiences. Maybe it could be worth a try for you?

Depending on what you are doing, you may be able to use the extend code interTrackFoam (https://github.com/Unofficial-Extend...interTrackFoam). This uses a "body-fitted mesh" approach. It may not work if you are interested in large deformations, or coalescence, etc.

To be honest, I don't have any great suggestions at this time. I ended up in this thread by asking the same questions you are.
SaufiSangi likes this.
liquidspoon is offline   Reply With Quote

Old   October 24, 2017, 16:05
Default
  #6
New Member
 
Abd Essamade Saufi
Join Date: Apr 2017
Posts: 14
Rep Power: 2
SaufiSangi is on a distinguished road
Yeah I am also having not the best experience with LSM-VOF coupling, but I try to go deeper in that.

It's the first time I see this "interTrackFoam". Is it still a VOF method or it's something like a Front-Tracking method? Did you try it?

Thank you very much anyway

Sangi
SaufiSangi is offline   Reply With Quote

Old   October 24, 2017, 16:16
Default
  #7
Member
 
Alex
Join Date: Jun 2011
Posts: 32
Rep Power: 8
liquidspoon is on a distinguished road
I have not tried interTrackFoam, but it is interesting. It updates the cell vertex positions so that the mesh is coincident with the interface. It seems like this would be a good thing because the pressure jump could be imposed directly across a cell face.
liquidspoon is offline   Reply With Quote

Old   October 30, 2017, 05:21
Default
  #8
Member
 
Knut Erik T. Giljarhus
Join Date: Mar 2009
Location: Norway
Posts: 33
Rep Power: 13
eric is on a distinguished road
It is very difficult to avoid spurious currents on general, unstructured meshes. If you are only interested in the drops/bubbles themselves, I would recommend using Cartesian grids, which makes it much easier to remove the spurious currents. See for instance the level set method with the ghost-fluid method to handle the discontinuities.
SaufiSangi likes this.
eric is offline   Reply With Quote

Old   October 30, 2017, 17:09
Default
  #9
New Member
 
Abd Essamade Saufi
Join Date: Apr 2017
Posts: 14
Rep Power: 2
SaufiSangi is on a distinguished road
Thank you Eric for your answer,

Yeah I noticed that. It would be very nice to use an Adaptive Mesh Refinement at the interface, but most of the times OF gives unstructured meshes, which of course enhance spurious currents.

I've seen some papers about the ghost-fluid method, but I cannot understand it properly, especially in OF environment. Do you have some nice articles about it?


Sangi
SaufiSangi is offline   Reply With Quote

Reply

Thread Tools
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 On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Setting the height of the stream in the free channel kevinmccartin CFX 10 July 9, 2015 21:36
Calculation of the Governing Equations Mihail CFX 7 September 7, 2014 06:27
error message cuteapathy CFX 14 March 20, 2012 07:45
Constant velocity of the material Sas CFX 15 July 13, 2010 08:56
mass flow in is not equal to mass flow out saii CFX 2 September 18, 2009 08:07


All times are GMT -4. The time now is 00:52.