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

Spurious waves and icAlpha

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

Like Tree1Likes
  • 1 Post By Sideshore

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   April 12, 2019, 09:32
Default Spurious waves and icAlpha
  #1
New Member
 
Jan Willem Krijger
Join Date: May 2010
Posts: 8
Rep Power: 16
Sideshore is on a distinguished road
Hello everybody!

Currently I'm working on resistance calculations of a IOM sailing boat (http://www.sailing.org/classesandequipment/IOM.php). As these are remote control boats of 1 meter long the calculations are done at "full" scale at 2 m/s.

I'm using OpenFOAM-v1812, the vessel is free to sink and trim using dynamicMotionSolverFvMesh and the sixDoFRigidBodyMotion, and a sail force is added as a sixDofRigidBodyMotionRestraint linearSpring.

Most of the settings I used are taken from the DTC_Hull tutorial.
I've attached the following files:
dynamicMeshDict
controlDict
fvSchemes
fvSolution

Unfortunately I noticed I have spurious waves in my freesurface and also the water interface along the hull doesn't look very good, see the attached pictures. I traced these problems to the icAlpha option of the alpha.water term in the fvSolution. Setting this to 0 eliminates the spurious waves, as you can see in the other attached pictures where I restarted the calculation and let it run for some timesteps.

The only information I can find about icAlpha is the following:
https://openfoamwiki.net/index.php/InterFoam

and this thread onfortunetely without any answers.
What is scAlpha vs. icAlpha vs interfaceCompression

Can anybody give me any advice which value for icAlpha (and scAlpha) I should use? And what is the range for icAlpha?
Any description would be much appreciated!

Thanks in advance for your comments!

regards,
Jan Willem
Attached Images
File Type: jpg ic_alpha_1.jpg (26.7 KB, 63 views)
File Type: png ic_alpha_1_hull.png (35.5 KB, 57 views)
File Type: png ic_alpha_0.png (188.6 KB, 55 views)
File Type: png ic_alpha_0_hull.png (33.4 KB, 47 views)
Attached Files
File Type: zip files.zip (3.2 KB, 18 views)
Sideshore is offline   Reply With Quote

Old   April 13, 2019, 14:03
Default similar issue but different problem !
  #2
Member
 
Join Date: Apr 2018
Location: UK
Posts: 78
Rep Power: 8
JM27 is on a distinguished road
Hi,


I am also interested in the functionality of icAlpha and scAlpha. It seems that documentation on this is very limited.



Unfortunately, spurious waves are a common problem in VOF simulations. A multitude of things could be triggering the presence of such spurious waves, so I think you need to do some detective work...



Some artificial compression like icAlpha and scAlpha can help but only one you are sure your setup is right



I would try different settings and compare the magnitude of these parasitic currents before trying to implement icAlpha or scAlpha. In any case I think their behaviour is similar to cAlpha, 0 for no compression, 1 for compression and >1 for enhanced compression. This should be explained in the source code.



Have you tried different Co number setting in your simulation?



Goodluck!
JM27 is offline   Reply With Quote

Old   April 18, 2019, 08:46
Default
  #3
New Member
 
Jan Willem Krijger
Join Date: May 2010
Posts: 8
Rep Power: 16
Sideshore is on a distinguished road
Hi,

Thanks for you reply!
For me it's not clear from the source code that the values for icAlpha & scAlpha should be 0 for no compression, 1 for compression and > 1 for enhanced compression. It seems obvious but for me that was/is just an assumption.

So I basically narrowed it down to the icAlpha coefficient. A value of 1 gives me spurious waves and 0 not. However for other cases with other mesh settings I did not have any problems with icAlpha.

This is inline with what I read somewhere (sorry cannot find it anymore) where using an icAlpha value > 0 can give spurious waves if the mesh has high aspect cells around the freesurface.

I tried lowering the maximum courant number in the freesurface interface which did not change the results. I tried lowering the tolerances and changing some numerical schemes without any luck.

Changing icAlpha to 0 was the only solution I found without going to first order methods.

Thanks again,
Jan Willem
lourencosm likes this.
Sideshore is offline   Reply With Quote

Old   April 18, 2019, 11:34
Default
  #4
Member
 
Join Date: Apr 2018
Location: UK
Posts: 78
Rep Power: 8
JM27 is on a distinguished road
Quote:
Originally Posted by Sideshore View Post
Hi,

Thanks for you reply!
For me it's not clear from the source code that the values for icAlpha & scAlpha should be 0 for no compression, 1 for compression and > 1 for enhanced compression. It seems obvious but for me that was/is just an assumption.

So I basically narrowed it down to the icAlpha coefficient. A value of 1 gives me spurious waves and 0 not. However for other cases with other mesh settings I did not have any problems with icAlpha.

This is inline with what I read somewhere (sorry cannot find it anymore) where using an icAlpha value > 0 can give spurious waves if the mesh has high aspect cells around the freesurface.

I tried lowering the maximum courant number in the freesurface interface which did not change the results. I tried lowering the tolerances and changing some numerical schemes without any luck.

Changing icAlpha to 0 was the only solution I found without going to first order methods.

Thanks again,
Jan Willem
Hi Jan,

If you have a look at alphaControls.H you will find some info on icAlpha and scAlpha.

Have you tried changing the cells in the free surface region so that they are close to 1? Maybe that would solve your issue with icAlpha working for some cases and not for others...

What value are you using for cAlpha?
JM27 is offline   Reply With Quote

Reply

Tags
icalpha, vof 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
icAlpha in OF 2.3.0 laurent98 OpenFOAM Running, Solving & CFD 0 June 29, 2014 20:22


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