fluent_udf_velocity_三维模型的例子

sorry

　　昨天分享了fluent udf help文件中关于入口velocity的例子的实现，帮助文件中例子是二维模型，为了更好地理解的udf代码的编写，今天做了个最简单的三维模型的例子。
　　模型如图所示，一个简单管道，坐标系原点位于入口面的圆心。设置入口处速度随圆半径的增大而程抛物线分布，也即靠近管壁出速度为0，圆心出速度最大，设为0.1m/s。



　　udf代码修改如下：

1. 
   
2. 
   
3. 　　#include "udf.h"
   
4. 
   
5. 　　DEFINE_PROFILE(inlet_x_velocity, thread, position)
   
6. 
   
7. 　　{
   
8. 
   
9. 　　real x[ND_ND]; 二维中该数组有两个元素，存储x，y坐标;三维模型，该数组存储x,y,z三个坐标。
   
10. 
    
11. 　　real r,R;
    
12. 
    
13. 　　face_t f;
    
14. 
    
15. 　　R=0.0075; 管道半径为7.5mm
    
16. 
    
17. 　　begin_f_loop(f,thread)
    
18. 
    
19. 　　{
    
20. 
    
21. 　　F_CENTROID(x, f, thread);
    
22. 
    
23. 　　r=sqrt(x[1]*x[1]+x[2]*x[2])/R; x[1],x[2]分别为入口面上face单元质心的y坐标值和z坐标值
    
24. 
    
25. 　　F_PROFILE(f, thread, position)=0.1*(1.0-r*r);
    
26. 
    
27. 　　}
    
28. 
    
29. 　　end_f_loop(f,thread)
    
30. 
    
31. 　　}
    

复制代码

　　设置求解为瞬态求解，液体为液态水，k-e模型，求解时间10s，结果如下所示：

　　速度分布图

　　中心平面上的动态压力分布图

　　经过二维和三维模型的udf实例实现，现在对fluent udf编程有了初步的认识，赶快尝试你的分析求解吧。



转自：http://blog.sina.com.cn/s/blog_15ae0ebcb0102wof8.html