|
这个还是自己动手吧
相关资料在论坛和Ftp都可以找到
下面发个例子:
- /*******************************************************************/
- /* vprofile.c */
- /* UDF for specifying steady-state velocity boundary condition */
- /*****************************************************************/
- #include "udf.h"
- #include "math.h"
- DEFINE_PROFILE(inlet_y_velocity, thread, position)
- { real x[ND_ND];
- real y;
- real z;
- face_t f;
- begin_f_loop(f, thread)!
- {
- F_CENTROID(x,f,thread);
- y= x[1];
- z= x[2];_
- F_PROFILE(f, thread, position) =2-(y*y+z*z)/(0.005*0.005)*2;
- }
- end_f_loop(f, thread)
- }
复制代码 |
|