|
[Sponsors] |
November 16, 2012, 18:31 |
3D velocity profile UDF
|
#1 |
New Member
Hadi
Join Date: Nov 2012
Location: Netherlands
Posts: 2
Rep Power: 0 |
Hi guys
I would like to define a 3D velocity profile but I don't know how. I read lots of 2D examples, but I am still confused about 3D. Is there any one who can help me? It is one of the samples for 2D How can I add the third dimension to this code? The desired velocity is Vz = (1-(x^2+y^2)/25) #include "udf.h" DEFINE_PROFILE(inlet_x_velocity, thread, position) { real x[ND_ND]; /* this will hold the position vector */ real y; face_t f; begin_f_loop(f, thread) { F_CENTROID(x,f,thread); y = x[1]; F_PROFILE(f, thread, position) = 20. - y*y/(.0745*.0745)*20.; } end_f_loop(f, thread) } |
|
November 17, 2012, 05:46 |
|
#2 | |
Senior Member
Rick
Join Date: Oct 2010
Posts: 1,016
Rep Power: 27 |
Quote:
I think this should work: Code:
#include "udf.h" DEFINE_PROFILE(my_velocity, thread, position) { face_t f; real x[ND_ND]; begin_f_loop(f, thread) { F_CENTROID(x,f,t); F_PROFILE(f, thread, position) = (1.-(x[0]*x[0]+x[1]*x[1])/25.); } end_f_loop(f, thread) } Last edited by ghost82; November 17, 2012 at 06:40. |
||
November 18, 2012, 08:23 |
|
#3 |
New Member
Hadi
Join Date: Nov 2012
Location: Netherlands
Posts: 2
Rep Power: 0 |
Thank you ghost82
How can I verify the position of my geometry? I upload my geometry. Thank you |
|
November 18, 2012, 10:03 |
|
#4 | |
Senior Member
Rick
Join Date: Oct 2010
Posts: 1,016
Rep Power: 27 |
Quote:
Your Vz velocity is a function of x and y coordinate, so each of the two cases will have different profiles, because x and y coordinates are different; according to this, you have to take into account the position of your inlet face in respect to the origin; Eventually, you can translate your geometry in pre-processing, or alternatively you can modify your Vz formula, considering a translation transformation. Daniele |
||
December 4, 2012, 02:32 |
|
#5 |
New Member
jianglei
Join Date: May 2012
Location: xiamen
Posts: 12
Rep Power: 14 |
ghost82,
I have a question to ask for you .In a non-regulation shape, for example ,ships , how can I get the center of gravity and buoyant coordinates in the Fluent software ? Please help me. Regards |
|
December 4, 2012, 02:32 |
|
#6 | |
New Member
jianglei
Join Date: May 2012
Location: xiamen
Posts: 12
Rep Power: 14 |
Quote:
I have a question to ask for you .In a non-regulation shape, for example ,ships , how can I get the center of gravity and buoyant coordinates in the Fluent software ? Please help me. Regards |
||
December 5, 2012, 04:13 |
|
#7 | |
Senior Member
Rick
Join Date: Oct 2010
Posts: 1,016
Rep Power: 27 |
Quote:
specify in details your problem, you may better open a new thread. |
||
December 5, 2012, 06:39 |
|
#8 |
New Member
jianglei
Join Date: May 2012
Location: xiamen
Posts: 12
Rep Power: 14 |
||
December 5, 2012, 10:16 |
|
#9 |
Senior Member
Rick
Join Date: Oct 2010
Posts: 1,016
Rep Power: 27 |
||
December 5, 2012, 20:29 |
|
#10 |
New Member
jianglei
Join Date: May 2012
Location: xiamen
Posts: 12
Rep Power: 14 |
||
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
UDF problems - velocity profile | Oli | Fluent UDF and Scheme Programming | 6 | October 24, 2016 11:38 |
Velocity profile UDF | egge24 | Fluent UDF and Scheme Programming | 2 | May 10, 2012 14:03 |
Logarithmic velocity profile | cfdworker | Fluent UDF and Scheme Programming | 0 | April 23, 2009 20:09 |
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 |