|
[Sponsors] |
UDF problem : inlet velocity in cyl. coord. system |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
June 15, 2004, 11:21 |
UDF problem : inlet velocity in cyl. coord. system
|
#1 |
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) } |
|
Thread Tools | Search this Thread |
Display Modes | |
|
|
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 09:46 |
Velocity inlet BC problem | Figd84 | FLUENT | 5 | October 16, 2009 01:46 |
UDF velocity and temperature | Raj | FLUENT | 3 | February 1, 2009 18: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 03:17 |