CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > FLUENT

3D free surface and VOF

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

Reply
 
LinkBack Thread Tools Display Modes
Old   February 5, 2003, 15:48
Default 3D free surface and VOF
  #1
Orome
Guest
 
Posts: n/a
Hi, we are trying to simulate a 3D channel with a free surface. Our problem rise when we want to include the udf function outlet_pressure:

#include "udf.h"

#include "sg_mphase.h"

#define RHO_WATER 998.2

#define RHO_AIR 1.225

#define G 9.81

DEFINE_PROFILE(outlet_pressure,thread,nv)

{ face_t f;

real ymax=-1e4;

real x[ND_ND];

Thread *t_phase2=THREAD_SUB_THREAD(thread,1);

begin_f_loop(f,thread) {

F_CENTROID(x,f,thread);

if ((ymax < x[1]) && (F_VOF(f,t_phase2) > 0.5)) ymax=x[1]; }

end_f_loop(f,thread)

begin_f_loop(f,thread) {

F_CENTROID(x,f,thread);

F_PROFILE(f,thread,nv)=0;

if (x[1] < ymax) F_PROFILE(f,thread,nv)=G*(RHO_WATER-RHO_AIR)*(ymax-x[1]); }

end_f_loop(f,thread) }

Then Fluent gives the message

Error: fluent.6.0.20 received a fatal signal (SEGMENTATION VIOLATION). Error Object: #f

and it exit. We have no idea where is the error.

Thank you in advance for your help.

  Reply With Quote

Old   February 6, 2003, 23:48
Default Re: 3D free surface and VOF
  #2
phd_y2k2001@yahoo.com
Guest
 
Posts: n/a
Just curious, why do you need to write udf function outlet_pressure??
  Reply With Quote

Old   February 7, 2003, 08:12
Default Re: 3D free surface and VOF
  #3
Orome
Guest
 
Posts: n/a
Well, we have a velocity inlet (water), a pressure inlet (air), a pressure outlet (water & air) and walls as boundary conditions. We want to fix the pressure at the outlet depending on hidrostatic pressure.
  Reply With Quote

Old   February 8, 2003, 08:35
Default Re: 3D free surface and VOF
  #4
Heydari
Guest
 
Posts: n/a
I think for your problem , you can use pressure inlet condition for the whole of your domain except your wall b.c.
  Reply With Quote

Old   February 10, 2003, 06:27
Default Re: 3D free surface and VOF
  #5
Orome
Guest
 
Posts: n/a
Maybe, but we think that it is not the poblem, because when we fix the height (and then we don't have to find it), the routine works. It fails when we try to find the position of the free surface.
  Reply With Quote

Old   February 16, 2003, 05:32
Default i think u can help me
  #6
varshavian
Guest
 
Posts: n/a
i want to model an open channel junction i have two problem 1-i have two phases,air and water at the inlet(velocity inlet).depth of water is about 0.3m and total depth of channel is 0.5m and the rest of it is fuul of air.what is the fraction of water at the entrance and outlet?is water primary phase? 2-how can i give velocity profile at the entrance of channel to program?i have experimental data. please tell me how i can model free surface? thanks
  Reply With Quote

Old   February 17, 2003, 13:07
Default Re: i think u can help me
  #7
Orome
Guest
 
Posts: n/a
The primary phase must be always the lighter phase (air) and then water will be the secondary phase. I don't know how to give a velocity profile at the entrance of the channel, I always give an average velocity. I am sorry ( ) In the "Boundary conditions --> Velocity inlet Zone" I fix the Volume Fraction for secondary phase (water) to 1.
  Reply With Quote

Reply

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
VOF free surface vortex simulation mahdiii FLUENT 3 May 14, 2013 15:33
UDF Defining VOF Free Surface at Outlet Alex Fluent UDF and Scheme Programming 13 August 8, 2012 16:50
VOF in free surface alilouu FLUENT 0 September 27, 2011 06:32
Numerical schemes for free surface flows (VOF) botp OpenFOAM 2 March 11, 2011 16:27
free surface of VOF and melting model? wanghong FLUENT 3 March 13, 2006 10:57


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