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

Density averaging the surface tension term in interFoam

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   March 30, 2022, 12:32
Default Density averaging the surface tension term in interFoam
  #1
Member
 
Venkat Ganesh
Join Date: May 2020
Location: Cincinnati, Ohio
Posts: 49
Rep Power: 4
Venky_94 is on a distinguished road
Hey multiphase researchers, I'm trying a standard modification to the surface tension term in interFoam which has been suggested in Brackbill's paper on CSF method, which is to include a density averaging term in the surface tension force term calculation. (I've shared an image of the equation)

In code within openfoam, the current surface tension force term without the density averaging is implemented as
Code:
fvc::interpolate(sigmaK())*fvc::snGrad(alpha1)
I obtained the code for density averaged surface tension term from the post (In attempt to decrease spurious currents in VOF) through Kevin's github adaption (https://github.com/floquation/OF-kva...faceProperties). It is written as
Code:
fvc::interpolate(sigmaK()*rho)*fvc::snGrad(alpha1)*2/(rho1+rho2)
He has shared a detailed explanation of how the equation is arrived and while it all made perfect sense, adopting the change to interFoam is giving me very poor results in terms of bubble shapes, bubble departure diameters and times. I'm finding that that adaption is causing a severely elongated bubble instead of the typical spherical shapes and that the departure diameters and times are overpredicted. I was wondering if someone else had any luck adopting the density averaging approach from Brackbill's paper and if there is an alternate way of expressing that in code?

I'm sharing an image of the bubble shape I obtained using the current code for reference. It would be great if someone could confirm if they've faced similar issues with adopting the change as well.

"Brackbill, J. U., Kothe, D. B., & Zemach, C. (1992). A continuum method for modeling surface tension. Journal of computational physics, 100(2), 335-354."
Attached Images
File Type: jpg Surface tension term.JPG (15.4 KB, 7 views)
File Type: jpg Bubble Shape.jpg (61.2 KB, 8 views)
Venky_94 is offline   Reply With Quote

Reply

Tags
csf model, interfoam, surface tension

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
Plotting Surface tension force with interfoam Mahmoud_aboukhedr OpenFOAM Programming & Development 5 December 1, 2016 09:37
Surface tension Interfoam nb977 OpenFOAM Running, Solving & CFD 1 March 9, 2016 04:02
[Gmsh] Error : Self intersecting surface mesh, computing intersections & Error : Impossible velan OpenFOAM Meshing & Mesh Conversion 3 October 22, 2015 12:05
VOF +surface tension force modeling+ open channel flow+cyclic region= fatal error? SJSW Fluent Multiphase 2 November 18, 2014 05:15
Flame Surface Density. Burn Ing Main CFD Forum 3 December 5, 2004 18:12


All times are GMT -4. The time now is 19:28.