 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 2 | 1/1 | 返回列表 |
| 查看: 973 | 回復: 1 | ||
[求助]
兩相流UDF編寫求助 已有1人參與
|
|
二維圓柱水滴撞擊問題，編寫UDF讀取壁面第二層網(wǎng)格的水滴速度、密度和體積分數(shù)等，現(xiàn)在已經(jīng)找到第二層網(wǎng)格了，但是讀取值的時候還是有問題，求各位大神指導！ #include "udf.h" DEFINE_ON_DEMAND(test) { Domain *d; d = Get_Domain(3);/*water droplet phase*/ face_t f; cell_t c; cell_t c0, c1= -1; /* 定義相鄰面index,是個整數(shù),內(nèi)部面矢量從c0指向c1 */ Thread *t0, *t1 = NULL; /*pointer to gouping of cell*/ real NV_VEC(A); real NV_VEC(A0); real NV_VEC(A1); real NV_VEC(A2); real NV_VEC(A3); real NV_VEC(A4); int ID = 18; Thread *t = Lookup_Thread(d, ID); Thread *tm = THREAD_SUPER_THREAD(t); /* Thread tm will be a superthread, i.e., mixture */ Thread **pt = THREAD_SUB_THREADS(tm); begin_f_loop(f, t) { F_AREA(A, f, t); c0 = F_C0(f, t); t0 = F_C0_THREAD(f, t); c1 = F_C1(f, t); /* Get cell on other side of face */ t1 = F_C1_THREAD(f, t); //循環(huán)C里面的每一個face 先確定 boundary face, fabf ；讀取法向量nbf；對每一個face，讀法向量；找到對面的face；確定對面的C，c2 計算通量，UC2，// int i = 3; int *n = 0; c_face_loop(c0, t0, i) /* loops over all faces of a cell */ { face_t f0, f1, f2, f3; Thread *tf0, *tf1, *tf2, *tf3; f0 = C_FACE(c0, t0, 0); tf0 = C_FACE_THREAD(c0, t0, 0); f1 = C_FACE(c0, t0, 1); tf1 = C_FACE_THREAD(c0, t0, 1); f2 = C_FACE(c0, t0, 2); tf2 = C_FACE_THREAD(c0, t0, 2); f3 = C_FACE(c0, t0, 3); tf3 = C_FACE_THREAD(c0, t0, 3); F_AREA(A0, f0, tf0); /*get vector A value*/ F_AREA(A1, f1, tf1); /*get vector A value*/ F_AREA(A2, f2, tf2); /*get vector A value*/ F_AREA(A3, f3, tf3); /*get vector A value*/ /*確定f1是第一層網(wǎng)格上平行于壁面的face*/ if (BOUNDARY_FACE_THREAD_P(tf1)==0) /*Get interior face values*/ { cell_t c2, c3 ; Thread *t2, *t3; c2 = F_C0(f1, tf1); t2 = F_C0_THREAD(f1, tf1); c3 = F_C1(f1, tf1); /* Get cell on other side of face */ t3 = F_C1_THREAD(f1, tf1); F_AREA(A4, f1, tf1); /*get vector A value*/ /*確定c3是壁面第2層網(wǎng)格*/ C_UDMI(c0, t0, 0) = C_VOF(c0, t0); C_UDMI(c0, t0, 1) = C_U(c0, t0); C_UDMI(c0, t0, 2) = C_R(c0, t0); /*能正常讀取c0的值，但是c3的讀不出來*/ C_UDMI(c0, t0, 3) = C_U(c3, t3); C_UDMI(c0, t0, 4) = C_R(c3, t3)； C_UDMI(c0, t0, 5) = C_VOF(c3, t3)； } else { } } } end_f_loop(f, t) } |
回復此樓» 猜你喜歡
0856求調(diào)劑
已經(jīng)有13人回復
-
0703化學
已經(jīng)有11人回復
-
329求調(diào)劑
已經(jīng)有7人回復
-
321求調(diào)劑
已經(jīng)有8人回復
-
0856材料化工調(diào)劑 總分330
已經(jīng)有12人回復
-
309求調(diào)劑
已經(jīng)有3人回復
-
070300化學354求調(diào)劑
已經(jīng)有4人回復
-
330分求調(diào)劑
已經(jīng)有3人回復
-
299求調(diào)劑
已經(jīng)有8人回復
-
一志愿北京理工大學本科211材料工程294求調(diào)劑
已經(jīng)有6人回復
2樓2016-05-31 14:01:33
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 嘟嘟張 的主題更新
| 2 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 070300化學354求調(diào)劑 +4 | 101次希望 2026-03-28 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 調(diào)劑310 +7 | 溫柔的晚安 2026-03-25 | 8/400 |
|
|
|
[考研] 材料與化工272求調(diào)劑 +9 | 阿斯蒂芬2004 2026-03-28 | 9/450 |
|
|
|
[考研] 求調(diào)劑 +3 | QiMing7 2026-03-25 | 4/200 |
|
|
|
[考研] 347求調(diào)劑 +3 | 山頂見α 2026-03-25 | 3/150 |
|
|
|
[考研] 085405 考的11408求各位老師帶走 +3 | Qiu學ing 2026-03-28 | 3/150 |
|
|
|
[考研] 0703化學/290求調(diào)劑/本科經(jīng)歷豐富/工科也可 +9 | 丹青奶蓋 2026-03-26 | 10/500 |
|
|
|
[考研] 求調(diào)劑推薦 材料 304 +15 | 荷包蛋hyj 2026-03-26 | 15/750 |
|
|
|
[考研] 考研調(diào)劑 +4 | Sanmu-124 2026-03-26 | 4/200 |
|
|
|
[考研] 266分求材料化工冶金礦業(yè)等專業(yè)的調(diào)劑 +4 | 哇呼哼呼哼 2026-03-26 | 4/200 |
|
|
|
[考研] 305求調(diào)劑 +5 | 哇盧卡庫 2026-03-26 | 5/250 |
|
|
|
[考研] 085601 材料工程 313分 求調(diào)劑 +5 | Ong3 2026-03-27 | 5/250 |
|
|
|
[考研] 0703化學求調(diào)劑 +3 | 丹青奶蓋 2026-03-26 | 5/250 |
|
|
|
[考研] 化學工程085602 305分求調(diào)劑 +17 | RichLi_ 2026-03-25 | 17/850 |
|
|
|
[考研] 總分322求生物學/生化與分子/生物信息學相關(guān)調(diào)劑 +5 | 星沉uu 2026-03-26 | 6/300 |
|
|
|
[考研] 332求調(diào)劑 +6 | 032500 2026-03-25 | 6/300 |
|
|
|
[考研] 【2026考研調(diào)劑】制藥工程 284分 求相關(guān)專業(yè)調(diào)劑名額 +4 | 袁奐奐 2026-03-25 | 8/400 |
|
|
|
[考研] 生物學學碩求調(diào)劑 +7 | 小羊睡著了? 2026-03-23 | 10/500 |
|
|
|
[考研] 一志愿國科過程所081700，274求調(diào)劑 +3 | 三水研0水立方 2026-03-23 | 3/150 |
|
|
|
[考研] 284求調(diào)劑 +3 | yanzhixue111 2026-03-23 | 6/300 |
|