
[Sponsors] 
January 20, 2007, 04:17 
3D UDF Paraboilc Velocity Profile (Can't Maintain)

#1 
Guest
Posts: n/a

I am simulating the air flow over a 3D building in a wind domain using a paraboilc velocity inlet profile.
The inlet profile is written by the following code, which is compiled into Fluent without any errors. After 2000 iterations,the solution is converged, the velocity inlet profile at a point near to the entry is plotted, which agrees well with the udf profile. However, as the velocity profile is plotted at a distance of 1000m away from the inlet, the profile is somehow distorted and does not have the same profile as the inlet profile. The size of the wind domain is 3500(x) x 3500(y) x 150(z). The building is at a distance of 1200m from the inlet. What can i do to maintain the Velocity Profile along the wind domain? Plz Help~~~~~ // UDF code for 3D parabolic Velocity Profile // #include "udf.h" DEFINE_PROFILE(inlet_x_velocity, thread, position) { real x[ND_ND]; real z; // Vertical Distance face_t f; begin_f_loop(f, thread) { F_CENTROID(x,f,thread); z = x[2]; F_PROFILE(f, thread, position) = 24.9 *pow(z,0.27); } end_f_loop(f, thread) } 

January 20, 2007, 05:07 
Re: 3D UDF Paraboilc Velocity Profile (Can't Maint

#2 
Guest
Posts: n/a

Hi,
most probably, this is because the wall function used with the turbulence model is not compatible with the Atmospheric Boundary Layer profile you've implemented in the UDF. The wall functions are valid for small geometries (compared to yours) and relatively small cells. You just tapped a problem tat is very often neglected, but that chouldn't be neglected! You will have to play with the wall functions (I heard that's possible in v6.3) or with the keps model constants... Can you provide some extra details on your project, simulation goals, findings, ... and we can probably help each other. I'll have to do an atmospheric dispersion problem soon, so it's a good time to start cooperating on these issues... talk to you soon! cheers, Laika, still orbiting 

January 20, 2007, 05:29 
Re: 3D UDF Paraboilc Velocity Profile (Can't Maint

#3 
Guest
Posts: n/a

Hi,
Thanks Laika for replying my message~ As i have posted in the previous email, i have to simulate the air flow over an isolated highrise building, in order to observe the flow changes in the wake region in terms of velocity profile and TKE. In my simulation, i am using the standard kepsilon model as the first step to obtain some preliminary results. The simluated results are satisfactory, like obseving recirculation behind the building. However, the problem is that the velocity profile at the front of the building is distorted and is actually quite different from the input parabolic velocity profile. I am afraid that my simulated results would be very different from the results having perfect parabolic profile. Since one of my requirement is to carry out the simulation using various wind speed and profiles, so it is really important for me to keep the velocity profile up to the building. I am using Fluent version 6.1.22 only, how can i modify the kepsilon model constants as Laika mentioned before? Can you tell more information on that? I can supply further information for my project if needed. Thx Sing 

February 25, 2010, 11:19 

#4 
New Member
srinivas
Join Date: Feb 2010
Posts: 9
Rep Power: 9 
dude can u pls help with writing udf .
Im running simulation with a cylinder ,sphere inside it for Non Newtonian fluid. I recently started working with fluent ,hope u could lead ur helping hand Dimensions of cylinder L=20cms D=4cms d(sphere)=2cms Parabolic Velocity equation be V=Vmax(1r/R)^(n+1/n) my code: #include "udf.h" DEFINE_PROFILE(inlet_x_velocity, thread, position) { real x[ND_ND]; real y,r,R; face_t f; begin_f_loop(f, thread) { F_CENTROID(x,f,thread); r = pow(x,2)+pow(y,2); // Vmax=4.1,R=o.o4// F_PROFILE(f, thread, position) = Vmax *pow(1r/R,(n+1)/n); } end_f_loop(f, thread) } Please tell the changes i have to make ,thank u 

September 27, 2010, 06:50 

#5  
Member
User_CFD
Join Date: Mar 2010
Posts: 48
Rep Power: 9 
Quote:


August 27, 2011, 11:36 
velocity inlet for 3d parabolic equation.

#6 
Member

hiii friends i am entirely new to using udf's in fluent. I have a cubical geometry, with two inlets on perpendicular adjacent faces. i have to implement udf to only one of the conditions.
i have two cases 1) parabolic velocity inlet ( space varient) 2) sinusodial velocity velocity inlet ( time variant) can any one pls send me the udf ( codes) ... 

August 31, 2011, 01:09 

#7 
Member

. i have worked out with 2d udf for parabolic velocity inlet.
please help me how this code will change with 3d. i wanna know how will be the equation considering 3rd dimmension. as in 2d we have V(x)=Vmax y^2/(disatnce from centre in y direction )^2 * Vmax. please help me. also give suggestion how to work out on sinusoidal inlet udf. thanks in advance 

September 1, 2011, 01:38 

#8 
Member

can you please explain how does this line fit int o your code
F_PROFILE(f, thread, position) = 24.9 *pow(z,0.27); . i want to use 3d parabolic inlet mine is a cubical cell 1*1*1. inlet i a circular face of o.1 dia. please tell me how to modify according to my problem 

September 17, 2014, 18:40 

#9 
New Member
bhatti
Join Date: Aug 2013
Posts: 2
Rep Power: 0 
could you please tell me
what is the exact equation you are using for developing Parabolic velocity profile in 3D. 

November 26, 2014, 08:58 
udf for inlet error

#10 
New Member
sitasma
Join Date: May 2014
Posts: 8
Rep Power: 4 
i have to simulate a wind turbine of 45 m radius. i have used this udf at the inlet but there are some error when i initialize from inlet as shown in the picture.can anyone please tell me why is this error. what can be done to correct it?
#include "udf.h" #define zo 10 #define uo 5 DEFINE_PROFILE(inlet_x_velocity,thread,index) { real x[ND_ND]; real z; face_t f; begin_f_loop(f,thread) { F_CENTROID(x,f,thread); z= x[2]; F_PROFILE(f,thread,index)=uo*pow(z/zo,0.16); } end_f_loop(f,thread) } 

October 14, 2015, 20:24 

#11 
New Member
Edmundo
Join Date: Jun 2014
Posts: 1
Rep Power: 0 
Hi, did you solved the problem?, if so, did you found an error rate between the perfect parabolic profile and the distorted parabolic profile?
Regards 

October 24, 2016, 04:49 
UDF for inlet temperature

#12 
New Member
mm
Join Date: May 2016
Posts: 24
Rep Power: 1 
Dear all
I have following UDF for inlet temperature, untill 1300s it takes correct values according to equation, but after 1300s values are higher and not accorrding to equation, like at 1301s it should have value of 405C but in simulation inlet temperature is 621C. I could not find the error in my UDF after lot of try. please check this and guide me. help please #include"udf.h" DEFINE_PROFILE(inlet_temperature,thread,position ) { face_t f; begin_f_loop(f,thread) { real t = RP_Get_Real("flowtime"); if (t <=1300.0 ) { F_PROFILE(f,thread,position) = 379.13 + 0.0005*t; } else if (1300.0 < t && t <= 1500.0 ) { F_PROFILE(f,thread,position)= 1.04289036878969*pow(10,10)*pow(t,6.0)+ 8.86126436853789*pow(10,7)*pow(t,5.0)3.13621260398811*pow(10,3)*pow(t,4.0)+5.91804640375908*pow(t,3.0)6.27969461279651*pow(10,3)*pow(t,2.0)+ 3.55273415252714*pow(10,6)*t  8.37223405676245*pow(10,8); } else { F_PROFILE(f,thread,position) = 9.51538261322402*pow(10,23)*pow(t,6) + 8.26192751387975*pow(10,18)*pow(t,5)2.85237398505875*pow(10,13)*pow(t,4)+4.97518353700886*pow(10,9)*pow(t,3)4.58733775886876*pow(10,5)*pow(t,2)+ 2.10251137071757*pow(10,1)*t +3.57252192344954*pow(10,2); } } end_f_loop(f,thread) } 

Thread Tools  
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
UDF problems  velocity profile  Oli  Fluent UDF and Scheme Programming  6  October 24, 2016 10:38 
How to specify the Velocity profile (eqn) by UDF  Anant  FLUENT  1  February 27, 2008 15:54 
Multi step transient UDF velocity profile problem  William177  FLUENT  1  February 3, 2008 07:47 
UDF problem : inlet velocity in cyl. coord. system  Jongdae Kim  FLUENT  0  June 15, 2004 11:21 
particle velocity and velocity profile  HGG  FLUENT  2  June 10, 2001 16:32 