CFD Online Logo CFD Online URL
Home > Forums > Fluent UDF and Scheme Programming

modified SurfaceTension (with filmthickness)

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

LinkBack Thread Tools Display Modes
Old   March 31, 2010, 11:13
Default modified SurfaceTension (with filmthickness)
New Member
Join Date: Mar 2010
Posts: 2
Rep Power: 0
Martin D! is on a distinguished road

I'm very new to UDF programming. The problem I'm working with is as follows: water should rinse down a vertical plate (with gravity) while gas comes from the bottom in a countercurrent way. It's know that the surface of the water will be covered with waves in a chaotic many. These waves are variable in time and place.
To simulate this behaviour i'm using a 2D-VOF model (yet the problem is 3D - but this cames later on) in a transient calculation. Everything works fine, but sometimes the film got split apart. This isn't observed at my Reynoldsnumbers used!
Now I try to increase the Surfacetension with an UDF when the film becames very small after a wavecrest. So I will account for the cohesion-forces which depend on the inter-molecular forces, and prevent naturally the film from getting ripped apart.
I can load the UDF in, hook it, but it doesn't seem to work.

Please have a look at the code. Maybe u can say me what's wrong with it ...
The formula used is from an article, I will give it a try ...

/*UDF (UserDefiniedFunction) for modified SurfaceTension with filmthickness*/

#include <udf.h>

DEFINE_PROPERTY(User_SurfaceTension,c,t) /*!!!Loops automatically over all cells!!!*/

real k = 2.3e-5;
real vof = C_VOF(c,t);
real delta;
real surfacetension;
real sigma0 = 0.07199; /*SurfaceTension at 25C (standard)*/

if (vof == 0.05)
delta = C_WALL_DIST(c,t); /*get filmthickness at vof=0.05*/
surfacetension = pow(k/delta,4) + sigma0; /*formula to increase Surfacetension esp. when delta is low*/
return surfacetension;

Thanks for every answer!

Martin D! is offline   Reply With Quote


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
Modified scalarTransportFoam, but result is not right , need help panda60 OpenFOAM 2 December 2, 2009 20:50
writing modified U, but U was not change at all ? panda60 OpenFOAM 2 November 18, 2009 03:30
LES and modified pressure clrsgfn CFX 0 November 2, 2009 12:21
Modified Equation for CFX algorithm Craig Johansen CFX 0 August 27, 2004 23:02

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