请教gallongallon，关于源项DEFINE_SOURCE在流场中某一区域温度大于600度加质问题

ruyanzhao

请教gallongallon高手，
  我现在想用DEFINE_SOURCE给流场中添加源项，添加位置是流场中靠近壁面一层网格中温度大于600度的网格加源项
  我用DEFINE_SOURCE循环流畅网格，再通过thread 循环得到靠近壁面的网格，如果温度大于600则加质，程序如下：
       DEFINE_SOURCE（mass,c,t,ds,eqn）
{
        face_tread=lookup_thread(domain,1);
   begin_f_loop(f,face_tread)
{
     T=F_T(f,face_tread);
     c0=F_C0(f,face_tread);
   if((c==c0)&&T>600)
    {
      source=10;
    }
}
end_f_loop(f,face_tread)
ds[eqn]=0;
return source;
}
如果通过解释法进行编译时可以使用，而通过compile法编译连接后
为什么加质时整个靠近壁面的网格全部加质呢？

owencaller

我也遇到相似的问题，有空可以探讨探讨。
我是将靠近壁面的一层流体单元单独作为一个流体域，然后再用源项法。也没算通。
好像DEFINE_SOURCE自动循环域中的所有单元的。

chinafrank

DEFINE_SOURCE是针对整个计算域中的Cell来完成的。可以试试通过Cell的编号或者Cell中心的坐标来筛选。

inneror

这个问题解决了吗？

我也正在困扰着

jinruibin

1. /***************** UDF ***********************/
   
2. /*
   
3. * FLUENT User Defined Function
   
4. *
   
5. * version 4.0
   
6. *
   
7. *
   
8. * Purpose
   
9. * -------
   
10. * Marks cells that contain coordinates in the a specified input file.
    
11. *
    
12. * Author
    
13. * ------
    
14. * William Wangard, Ph.D.
    
15. * Fluent, Inc.
    
16. * Feb 2003
    
17. *
    
18. * Note: Fluent, Inc. does not guarantee the rigor of this User-defined function
    
19. *
    
20. *
    
21. * Usage
    
22. * -----
    
23. *
    
24. * Data input format
    
25. * -----------------
    
26. * The format for the UDF input file (interp-input.dat) is the following:
    
27. *
    
28. * x0, y0, z0
    
29. * x1, y1, z1
    
30. * ...
    
31. *
    
32. *
    
33. * Note, if you are using 2D, then you only enter the x- and y-coordinates.
    
34. *
    
35. *
    
36. * Execution
    
37. * ---------
    
38. * There are 2 functions in this packages: mark, source.
    
39. *
    
40. *
    
41. * Define on Demand function "mark" is used to read the data points from the inpute file
    
42. * and locate the points in the domain. Upon finding each point, a UDM is set
    
43. * to positive value in that cell.
    
44. *
    
45. * In FLUENT, you will separate the cell thread by UDM value. Then, you apply
    
46. * the Mass source to the source-contatining cell thread.
    
47. *
    
48. *
    
49. */
    
50. 
    
51. 
    
52. #include "udf.h"
    
53. #include "surf.h"
    
54. #define MAXPOINTS 20000
    
55. #define UDM_SOURCE 0
    
56. 
    
57. boolean check(void);
    
58. 
    
59. static int np = 0;
    
60. real coordinates[MAXPOINTS][ND_ND] = {{0.}};
    
61. 
    
62. 
    
63. struct interpolation_point{
    
64. cell_t c;
    
65. Thread* t;
    
66. };
    
67. 
    
68. struct interpolation_point point_list[MAXPOINTS];
    
69. 
    
70. 
    
71. boolean check(void)
    
72. {
    
73. boolean ret = TRUE;
    
74. if (sg_udm < 1)
    
75. {
    
76. ret = FALSE;
    
77. Message("Error: You must define at least 1 UDM\n");
    
78. }
    
79. if (!Data_Valid_P())
    
80. {
    
81. ret = FALSE;
    
82. Message("Error: You must initialize the data field\n");
    
83. }
    
84. return ret;
    
85. }
    
86. 
    
87. 
    
88. DEFINE_ON_DEMAND(mark_reset)
    
89. {
    
90. Domain*d = Get_Domain(1);
    
91. Thread*t;
    
92. cell_t c;
    
93. 
    
94. if (!check()) return;
    
95. 
    
96. thread_loop_c(t,d)
    
97. {
    
98. begin_c_loop(c,t)
    
99. {
    
100. C_UDMI(c,t,UDM_SOURCE) = 0.0;
     
101. }
     
102. end_c_loop(c,t);
     
103. }
     
104. }
     
105. 
     
106. DEFINE_ON_DEMAND(mark)
     
107. {
     
108. 
     
109. #if !RP_HOST
     
110. Domain *d = Get_Domain(1);
     
111. cell_t c;
     
112. CX_Cell_Id cx_cell;
     
113. Thread* t;
     
114. int points_found = 0, total_points_found=0;
     
115. #endif
     
116. 
     
117. #if !RP_NODE
     
118. FILE* input;
     
119. #endif
     
120. int n;
     
121. 
     
122. 
     
123. /* Check for UDM */
     
124. if (!check()) return;
     
125. 
     
126. /* READ COORDINATES */
     
127. #if !RP_NODE
     
128. 
     
129. /* Open input file */
     
130. if (!(input = fopen("interp-input.dat","r")))
     
131. {
     
132. Message0("\nWarning: Could not open interpolation input file...\n");
     
133. return;
     
134. }
     
135. 
     
136. /* Initialize */
     
137. for(n=0; n<MAXPOINTS; n++)
     
138. {
     
139. point_list[n].c = 0;
     
140. point_list[n].t = NULL;
     
141. }
     
142. 
     
143. /* Read points from input file */
     
144. n = -1;
     
145. while (!feof(input))
     
146. {
     
147. n++;
     
148. #if RP_DOUBLE
     
149. #if RP_3D
     
150. fscanf(input,"%lg %lg %lg", &coordinates[n][0], &coordinates[n][1], &coordinates[n][2])
     
151. ;
     
152. #else
     
153. fscanf(input,"%lg %lg", &coordinates[n][0], &coordinates[n][1]);
     
154. #endif
     
155. #else
     
156. #if RP_3D
     
157. fscanf(input,"%g %g %g", &coordinates[n][0], &coordinates[n][1], &coordinates[n][2])
     
158. ;
     
159. #else
     
160. fscanf(input,"%g %g", &coordinates[n][0], &coordinates[n][1]);
     
161. #endif
     
162. #endif
     
163. 
     
164. }
     
165. np = n;
     
166. 
     
167. 
     
168. /* Check np with MAXPOINTS */
     
169. if (np > MAXPOINTS)
     
170. {
     
171. Message0("ERROR: You must recompile interpolate UDF with MAXPOINTS at least %i\n", np)
     
172. ;
     
173. return;
     
174. }
     
175. 
     
176. 
     
177. /* Close input file */
     
178. fclose(input);
     
179. 
     
180. 
     
181. 
     
182. #endif
     
183. 
     
184. 
     
185. 
     
186. /* FIND COORDINATES IN CELLS */
     
187. 
     
188. /* Pass coordinates and np to the compute nodes */
     
189. host_to_node_real(&coordinates[0][0],ND_ND*MAXPOINTS);
     
190. host_to_node_int_1(np);
     
191. 
     
192. #if !RP_HOST
     
193. 
     
194. /* Do computations */
     
195. for(n=0;n<np;n++)
     
196. {
     
197. thread_loop_c(t,d)
     
198. {
     
199. begin_c_loop_int(c,t)
     
200. {
     
201. if (SV_is_point_in_cell(&cx_cell, c, t, coordinates[n]))
     
202. {
     
203. point_list[n].c = RP_CELL(&cx_cell);
     
204. point_list[n].t = RP_THREAD(&cx_cell);
     
205. 
     
206. /* SET UDM HERE */
     
207. C_UDMI(point_list[n].c, point_list[n].t, UDM_SOURCE) = 1.0;
     
208. 
     
209. points_found++;
     
210. goto label;
     
211. }
     
212. }
     
213. end_c_loop_int(c,t);
     
214. }
     
215. label: continue;
     
216. }
     
217. 
     
218. total_points_found += points_found;
     
219. 
     
220. #if PARALLEL
     
221. total_points_found = PRF_GISUM1(total_points_found);
     
222. #endif
     
223. 
     
224. 
     
225. /* PRINT MESSAGE */
     
226. if (np != total_points_found)
     
227. Message0("\n\n Warning.... %i points found in domain out of %i points in input file\n",
     
228. total_points_found, np);
     
229. else
     
230. Message0("\n\n Interpolate... all %i points located!\n", np);
     
231. 
     
232. #endif
     
233. 
     
234. }
     
235. 
     
236. 
     
237. 
     
238. DEFINE_SOURCE(mass_source,c,t,dS,eqn)
     
239. {
     
240. 
     
241. real source = 0.0;
     
242. real mdot = 1.0;
     
243. source = (C_UDMI(c,t,UDM_SOURCE) == 1.0) ? (mdot/C_VOLUME(c,t)) : 0.0;
     
244. dS[eqn] = 0.0;
     
245. return source;
     
246. 
     
247. }
     
248. 
     
249. 
     
250. /************************ END OF UDF *****************************/
     

复制代码

CAI

先将代码如帖，拜谢高手指点！
#include "udf.h"
DEFINE_SOURCE(energy_source1,c,t,ds,eqn)
{
real x[ND_ND];
real a,b,cwz;
real source;
cell_t c;
Thread *c_thread;
begin_c_loop(c,c_thread)
{
C_CENTROID(x,c,t);
a=x[0],b=x[1];
cwz=sqrt(a*a+b*b);
source=482252.4/cwz;
ds[eqn]=0;
return source;
}
end_c_loop(c,c_thread);
}

八怪

写的不错 但是还是 有很多的疑惑的啊

luckyhuman1

你试过上面那段程序没？我试了下编译出错，boolean那块通不过，还有啊，您看懂了这段不，我怎么好多看不懂，查UDF帮助很多都没介绍

wangxw0315

#include "udf.h"
DEFINE_SOURCE(energy_source1,c,t,ds,eqn)
{
real x[ND_ND];
real source;
C_CENTROID(x,c,t);
source=482252.4/sqrt(x[0]*x[0]+x[1]*x[1])
ds[eqn]=0;
return source;
}