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

implement level set to compressibleInterFoam

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   January 10, 2018, 10:09
Default implement level set to compressibleInterFoam
  #1
cyw
Member
 
Yuanwei Cao
Join Date: Jun 2017
Location: Munich,Germany
Posts: 54
Rep Power: 8
cyw is on a distinguished road
Dear Foamers,

Does anybody implement level set method into compressibleInterFoam? Could you please tell me some theory guide?

Best Regards
CYW
cyw is offline   Reply With Quote

Old   January 10, 2018, 13:29
Default
  #2
Senior Member
 
Taher Chegini
Join Date: Nov 2014
Location: Houston, Texas
Posts: 125
Rep Power: 12
Taataa is on a distinguished road
I have been working on it in my free time so I haven't had a significant progress yet.

Here is a good starting point. Look for "Coupled Level-Set with VOF interFoam, Sankar Menon (2014)" in the link and you can find the report as well as the files.
Taataa is offline   Reply With Quote

Old   January 11, 2018, 09:50
Default level set in compressibleInterFoam
  #3
cyw
Member
 
Yuanwei Cao
Join Date: Jun 2017
Location: Munich,Germany
Posts: 54
Rep Power: 8
cyw is on a distinguished road
Quote:
Originally Posted by Taataa View Post
I have been working on it in my free time so I haven't had a significant progress yet.

Here is a good starting point. Look for "Coupled Level-Set with VOF interFoam, Sankar Menon (2014)" in the link and you can find the report as well as the files.
Thank you. I have tried what you said. But for compressibleInterFoam, I think it is different. I tried to implement level set in compressibleInterFoam, but the result it not good as the interFoam. The parasitic currents become more obvious...
cyw is offline   Reply With Quote

Old   January 11, 2018, 10:23
Default
  #4
Senior Member
 
Taher Chegini
Join Date: Nov 2014
Location: Houston, Texas
Posts: 125
Rep Power: 12
Taataa is on a distinguished road
I will resume the work next month, will let you know if I have any achievement in this regard and if interested we can work together.
Taataa is offline   Reply With Quote

Old   January 11, 2018, 11:17
Default
  #5
cyw
Member
 
Yuanwei Cao
Join Date: Jun 2017
Location: Munich,Germany
Posts: 54
Rep Power: 8
cyw is on a distinguished road
Quote:
Originally Posted by Taataa View Post
I will resume the work next month, will let you know if I have any achievement in this regard and if interested we can work together.
Thank you, but I cannot wait until next month...I have to do it right now
cyw is offline   Reply With Quote

Old   January 11, 2018, 14:52
Default
  #6
Senior Member
 
Taher Chegini
Join Date: Nov 2014
Location: Houston, Texas
Posts: 125
Rep Power: 12
Taataa is on a distinguished road
Sorry, I've got my hands full this month. If you have any success please let me know. If you have any specific errors regarding implementation also maybe I can help.
Taataa is offline   Reply With Quote

Old   January 12, 2018, 11:28
Default
  #7
cyw
Member
 
Yuanwei Cao
Join Date: Jun 2017
Location: Munich,Germany
Posts: 54
Rep Power: 8
cyw is on a distinguished road
Quote:
Originally Posted by Taataa View Post
Sorry, I've got my hands full this month. If you have any success please let me know. If you have any specific errors regarding implementation also maybe I can help.
No problem. I am simulating a bubble in still water without gravity. For interFoam and sclsVOFFoam, the result is OK and the parasitic current can be controlled. But for the compressibleInterFoam and sclsVOF_compressible, the result is not OK. The parasitic currents exist. Coupled level set method can not solve this problem. I am not sure if the method or my programming is wrong.
cyw is offline   Reply With Quote

Old   January 12, 2018, 12:54
Default
  #8
Senior Member
 
Taher Chegini
Join Date: Nov 2014
Location: Houston, Texas
Posts: 125
Rep Power: 12
Taataa is on a distinguished road
Parasitic currents are a result of the way surface tension is implemented in OF so you have pay extra attention to the way compressibleInterFoam deals with surface tension modeling as it is different than interFoam. So implementing LS is not as straight forward as interFoam.
Taataa is offline   Reply With Quote

Old   January 12, 2018, 13:42
Default surface tension in compressibleInterFoam
  #9
cyw
Member
 
Yuanwei Cao
Join Date: Jun 2017
Location: Munich,Germany
Posts: 54
Rep Power: 8
cyw is on a distinguished road
Quote:
Originally Posted by Taataa View Post
Parasitic currents are a result of the way surface tension is implemented in OF so you have pay extra attention to the way compressibleInterFoam deals with surface tension modeling as it is different than interFoam. So implementing LS is not as straight forward as interFoam.
Thank you for your information. I think they are the same. In compressibleInterFoam, the surface tension is calculated in interfaceProperties, which is the same in interFoam. And the surface tension is needed in UEqn.H,
which is interface.surfaceTensionForce().
cyw is offline   Reply With Quote

Old   January 12, 2018, 14:56
Default
  #10
Senior Member
 
Taher Chegini
Join Date: Nov 2014
Location: Houston, Texas
Posts: 125
Rep Power: 12
Taataa is on a distinguished road
In case of a quiescent flow without gravity, lack of balance in surface tension and pressure gradient forces, and errors in curvature calculations is the source of the currents.

So I would suggest to look into this case and see how compressibility makes a difference in this simple case. At least that's what I was trying to do and understand. Just put a square bubble in a domain where the boundaries are 4 times the radius away from the bubble, remove the gravity and try it with the models that you have and figure out what are the differences.
Taataa is offline   Reply With Quote

Old   January 13, 2018, 15:49
Default
  #11
cyw
Member
 
Yuanwei Cao
Join Date: Jun 2017
Location: Munich,Germany
Posts: 54
Rep Power: 8
cyw is on a distinguished road
Quote:
Originally Posted by Taataa View Post
In case of a quiescent flow without gravity, lack of balance in surface tension and pressure gradient forces, and errors in curvature calculations is the source of the currents.

So I would suggest to look into this case and see how compressibility makes a difference in this simple case. At least that's what I was trying to do and understand. Just put a square bubble in a domain where the boundaries are 4 times the radius away from the bubble, remove the gravity and try it with the models that you have and figure out what are the differences.
Thank you for your information. I have simulating a case similar to you without gravity. The parasitic current exists for the result from compressibleInterFoam and also for that coupled with level set method. I doubt where I am wrong or the compressibleInterFoam coupled with level set cannot solve this parasitic current...
cyw is offline   Reply With Quote

Old   January 13, 2018, 16:03
Default
  #12
Senior Member
 
Taher Chegini
Join Date: Nov 2014
Location: Houston, Texas
Posts: 125
Rep Power: 12
Taataa is on a distinguished road
I see.
In general I don't think coupling with LS provide a significant advantages in this case, otherwise they would have implemented it officially. The CIF solver is complex, adding more complexity wouldn't help I suppose.
There are some alternatives though. If you need sharp interface you can try using isoAdvector solver. It's already a part of OF v1706, here.
There are some nice comparison between interFoam and CLSVOF in this paper. interFoam performs very well.
Also, you can try implementing a force balancing algorithm which are specifically for this purpose, for example this paper.
Taataa is offline   Reply With Quote

Old   January 18, 2018, 07:36
Default
  #13
cyw
Member
 
Yuanwei Cao
Join Date: Jun 2017
Location: Munich,Germany
Posts: 54
Rep Power: 8
cyw is on a distinguished road
Have you simulated the dambreak case in the tutorial with sclsVOFFoam? What is the result?
In the report, the sclsVOFFoam can only handle zerogradient boundary condionts? What does that mean? How to handle other boundary conditions?
cyw is offline   Reply With Quote

Old   January 18, 2021, 06:27
Default
  #14
New Member
 
Armin Alavi
Join Date: May 2019
Location: Tehran
Posts: 22
Rep Power: 7
ArminAlavi is on a distinguished road
Quote:
Originally Posted by cyw View Post
Have you simulated the dambreak case in the tutorial with sclsVOFFoam? What is the result?
In the report, the sclsVOFFoam can only handle zerogradient boundary condionts? What does that mean? How to handle other boundary conditions?
Hello Yuanwei
Did you manage to find out how to handle other boundary conditions in the sclsVOFFoam solver?
ArminAlavi is offline   Reply With Quote

Reply

Tags
compressibleinterfoam, level set, multiphase

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
[ICEM] ICEM Scripting Issues tylerplowright ANSYS Meshing & Geometry 33 September 27, 2021 16:35
OF 1.6 | Ubuntu 9.10 (64bit) | GLIBCXX_3.4.11 not found piprus OpenFOAM Installation 22 February 25, 2010 13:43
Problem on high density ratio in Level Set method Kai Yan Main CFD Forum 10 December 25, 2007 06:12
Env variable not set gruber2 OpenFOAM Installation 5 December 30, 2005 04:27
How to set environment variables kanishka OpenFOAM Installation 1 September 4, 2005 10:15


All times are GMT -4. The time now is 11:56.