
[Sponsors] 
August 27, 2012, 02:43 
free surface modelling using VOF

#1 
New Member
Join Date: Aug 2012
Posts: 9
Rep Power: 5 
I'm trying to write a code in Frotran to simulate free water surface in 2DV using Volume Of Fluid method. My code isn't giving the right answer.
I assume that cells with free surface have known pressures. For the rest of the cells using the momentum equations and free surface continuity equation, I find the new pressure terms. Is it right to use free surface equation for cells that aren't free surface? if not what should I do? Thanks for helping. 

August 28, 2012, 06:21 

#2 
Senior Member
Join Date: Aug 2011
Posts: 251
Rep Power: 7 

August 28, 2012, 08:17 

#3 
New Member
Join Date: Aug 2012
Posts: 9
Rep Power: 5 
I mean the continuity equation that is used for the free surface:


August 28, 2012, 10:24 

#4 
Senior Member
Join Date: Aug 2011
Posts: 251
Rep Power: 7 

August 28, 2012, 12:10 

#5 
New Member
Join Date: Aug 2012
Posts: 22
Rep Power: 5 
If I remember correctly you will always get some spurious currents when using the VoF equation in 2D. I'm not sure if that is what you are experiencing or not, if it is checkout the openFOAM tutorials, I think there is one there which demonstrates what I mean.
Mike 

August 28, 2012, 15:44 

#6 
Senior Member
Join Date: Aug 2011
Posts: 251
Rep Power: 7 
VOF method is based on the solving of this equation:
dF/dt +udF/dx+vdF/dy=0 here d is the partial derivative, F is the volume fraction. When F(i,j)=0 cell (i,j) is fully occupied bythe first fluid, when F(i,j)=1 cell (i,j) is fully occupied by the second fluid. when 0<F(i,j)<1then cell (i,j) is occupied by the two fluids and then the interface cross the cell. Now I don't know what you are talking about with your "free surface continuity equation" 

August 29, 2012, 01:50 

#7 
New Member
Join Date: Aug 2012
Posts: 9
Rep Power: 5 
I've written a code solving this equation you mentioned using young method which works accurately.
what I'm trying to do now is to solve the continuity and momentum equations in order to find the new pressures and velocities in the new time step and then solve the VOF equation with the new velocities. 

August 29, 2012, 01:55 

#8 
New Member
Join Date: Aug 2012
Posts: 9
Rep Power: 5 
Mike I checked it out, but it wasn't related to my problem, thanks anyway.


August 29, 2012, 01:58 

#9 
New Member
Join Date: Aug 2012
Posts: 22
Rep Power: 5 
Okay, no worries, sorry for sending you down a blind alley!
Mike 

August 29, 2012, 02:56 

#10 
Member
Ren/Xingyue
Join Date: Jan 2010
Location: Nagoya , Japan
Posts: 44
Rep Power: 7 
nik,
The surface tension is not very important in single phase model. You don't consider the air phase in your model. Pressure and velocity were calculated by boundary condition at the free surface, so the sharp change across the free surface won't affect your calculation. See SOLAVOF and NASAVOF as a reference. If you want to use free surface tension, the surface force should be changed into body force and only exist in the smooth area based on CSF model, using a delta function. Ren 

August 29, 2012, 07:43 

#11  
Senior Member
Join Date: Aug 2011
Posts: 251
Rep Power: 7 
Quote:
Quote:
so you should solve a poisson equation for pressure which looks like DIV(1/rho GRAD P)= DIV(U*) Many papers exist on this topic. or you can use SIMPLE algorithm alos for variable density flow. 

Thread Tools  
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
UDF Defining VOF Free Surface at Outlet  Alex  Fluent UDF and Scheme Programming  13  August 8, 2012 16:50 
Free surface boundary condition in the VOF metod  ves  Main CFD Forum  2  March 8, 2011 16:23 
Free Surface Modelling and Reversed Flow Error  Javier González  FLUENT  4  August 16, 2004 11:37 
Free surface modelling under CFX 4  Dahong  CFX  2  June 22, 2002 23:07 
CFX 4.4 New free surface option  Viatcheslav Anissimov  CFX  0  April 3, 2002 06:27 