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

UDF problem : inlet velocity in cyl. coord. system

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

Reply
 
LinkBack Thread Tools Display Modes
Old   June 15, 2004, 11:21
Default UDF problem : inlet velocity in cyl. coord. system
  #1
Jongdae Kim
Guest
 
Posts: n/a
Computational domain of the (tornado) simulation is defined in the 3-dimensional cylindrical coordinate system.

I'd like to change the uniform inlet velocities (radial, swirl and axial velocities) to the boundary layer type velocities.

Please let me know how to define such velocities using UDF.

************************************************** ************************************************** ************************************************** ************************************************** ****************

Why boundary layer type?

From the 2-dim. axisymmetry analysis of a tonado, the uniform inlet boundary condition gives more intensive tornadic core than the case of using the boundary layer type inlet velocity condition. Also, when boundary layer type inlet velocity is used, the velocity profiles show similar trend with the experimental data from tornadic chamber test.

************************************************** ************************************************** ************************************************** ************************************************** ****************

The following is the UDF used for the inlet velocity profiles for a 2-dim axisym. analysis.

/************************************************** **/

/* UDF for specifing the inlet */ /* radial & swirl velocity profile */ /* for inlet boundary condtions */

/************************************************** **/

#include "udf.h"

#define ALFA 0.113

#define ZH 0.025

#define VH 0.3

#define VISC 1.7894e-05

#define DEFU 1.00 /* deffuser effect compesstor */

#define S 0.28 /* swirl ratio */

/* Profile for radial-velocity */

DEFINE_PROFILE(r_velocity, thread, position)

{

real x[ND_ND];

real z;

face_t f;

begin_f_loop(f, thread)

{

F_CENTROID(x, f, thread);

z = x[0];

F_PROFILE(f, thread, position) = (-DEFU*VH)*pow (z/ZH, ALFA);

}

end_f_loop(f, thread)

}

/* Profile for swirl-velocity */

DEFINE_PROFILE(s_velocity, thread, position)

{

real x[ND_ND];

real z;

real V;

face_t f;

begin_f_loop(f, thread)

{

F_CENTROID(x, f, thread);

z = x[0];

V = S*2.0*VH;

F_PROFILE(f, thread, position) = (DEFU*V)*pow(z/ZH, ALFA);

}

end_f_loop(f, thread)

}

  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
How to define a parabolic velocity inlet of particles with UDF zumaqiong FLUENT 2 February 22, 2010 10:46
Velocity inlet BC problem Figd84 FLUENT 5 October 16, 2009 01:46
UDF velocity and temperature Raj FLUENT 3 February 1, 2009 19:29
Variables Definition in CFX Solver 5.6 R P CFX 2 October 26, 2004 02:13
Cannot set correct inlet velocity by UDF?!? Bob FLUENT 5 March 27, 2003 04:17


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