 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 2044 | 回復(fù): 4 | ||
| 【懸賞金幣】回答本帖問題，作者cfd菜雞將贈送您 40 個(gè)金幣 | ||
cfd菜雞新蟲 (初入文壇)
|
[求助]
這是我編寫濕空氣冷凝的udf導(dǎo)入fluent后初始化就無法完成，是不是死循環(huán)，求大神指導(dǎo)
|
|
|
#include "udf.h" #include "sg.h" #include "sg_mphase.h" #include "flow.h" #include "mem.h" #define LAT_HT 2.257e6 #define WALL_NUMBER 3 #define scdomain_index 3 #define prim_index 0 #define index_evap_primary 1 int phase_domain_index; real face_center[ND_ND], cell_center[ND_ND], a[ND_ND], distance; real p_op, p, diffusivity; real cell_vapor_pressure, wall_temper, wall_pressure, wall_mass_fraction; real NV_VEC(A); real area, area_density, mass_transfer_coeff; real water_temper, water_satpressure, water_mass_fraction; real vapor_density; int n; Domain *subdomain; cell_t c; Thread *t; Thread **pt; Thread *tp; Thread *ts; Thread *tf, *tb; face_t fa; face_t fb; real source; DEFINE_SOURCE(vap_src, cell, first, dS, eqn) { Thread *mixer, *sec_th; real m_dot_first; real mass_dot; mixer = THREAD_SUPER_THREAD(first); sec_th = THREAD_SUB_THREAD(mixer, 1); real Db; real water_mass_fiction; real area; real volume; begin_c_loop(cell,mixer) { volume = C_VOLUME(cell, mixer);//網(wǎng)格體積// c_face_loop(cell, mixer, n)//對所有cell面循環(huán)// { fa = C_FACE(cell, mixer, n);//獲取面指針// tf = C_FACE_THREAD(cell, mixer, n);//獲取面所在thread指針// if (THREAD_ID(tf) == WALL_NUMBER) //判斷是否為壁面網(wǎng)格// { if (C_VOF(cell, sec_th) == 1.) //判斷液體是否充滿壁面// mass_dot = 0.0; //質(zhì)量源為0// else { C_CENTROID(cell_center, cell, mixer);//獲取cell中心坐標(biāo)// F_CENTROID(face_center, fa, tf);//獲取面中心坐標(biāo)// NV_VV(a, =, cell_center, -, face_center);//體心到面心矢量// distance = NV_MAG(a);//體心到面心距離// F_AREA(A, fa, tf);//面積矢量// area = NV_MAG(A);//面積// wall_temper = F_T(fa, tf) - 273.15;//獲取壁面溫度攝氏度// Db = (4.709 + 0.196 * wall_temper + 0.006122 * wall_temper * wall_temper + 0.0004386 * wall_temper * wall_temper * wall_temper) / 1000;//壁面溫度對應(yīng)飽和含濕量// water_mass_fiction = C_YI(cell, first, 0);//水蒸氣質(zhì)量分?jǐn)?shù)// if (water_mass_fiction > Db) { mass_dot = -C_R(cell, first) * 0.256 * 0.0001 * (water_mass_fiction - Db) / distance * area/volume;//質(zhì)量源求解公式// dS[eqn] = 0; } } } } } end_c_loop(cell,mixer) return mass_dot; } DEFINE_SOURCE(liq_src, cell, second, dS, eqn) { Thread *mix_th, *fir_th; real mass_dot; real Db; real water_mass_fiction; real area; real volume; mix_th = THREAD_SUPER_THREAD(second); fir_th = THREAD_SUB_THREAD(mix_th, 1); begin_c_loop(cell, mix_th) { volume = C_VOLUME(cell, mix_th);//網(wǎng)格體積// c_face_loop(cell, mix_th, n) { fa = C_FACE(cell, mix_th, n); tf = C_FACE_THREAD(cell, mix_th, n); if (THREAD_ID(tf) == WALL_NUMBER) { if (C_VOF(cell, fir_th) == 1.) mass_dot = 0.0; else { C_CENTROID(cell_center, cell, mix_th); F_CENTROID(face_center, fa, tf); NV_VV(a, =, cell_center, -, face_center); distance = NV_MAG(a); F_AREA(A, fa, tf);//面積矢量// area = NV_MAG(A);//面積// wall_temper = F_T(fa, tf) - 273.15; Db = (4.709 + 0.196 * wall_temper + 0.006122 * wall_temper * wall_temper + 0.0004386 * wall_temper * wall_temper * wall_temper) / 1000; water_mass_fiction = C_YI(cell, second, 0); if (water_mass_fiction > Db) { mass_dot = C_R(cell, second) * 0.256 * 0.0001 * (water_mass_fiction - Db) / distance * area / volume; dS[eqn] = 0; } } } } } end_c_loop(cell, mixer) return mass_dot; } DEFINE_SOURCE(enrg_src, cell, mixer, dS, eqn) { Thread* first, * sec_th; real m_dot; real mass_dot; first = THREAD_SUB_THREAD(mixer, 0); sec_th = THREAD_SUB_THREAD(mixer, 1); real Db; real water_mass_fiction; real area; real volume; begin_c_loop(cell, mixer) { c_face_loop(cell, mixer, n) { fa = C_FACE(cell, mixer, n); tf = C_FACE_THREAD(cell, mixer, n); if (THREAD_ID(tf) == WALL_NUMBER) { if (C_VOF(cell, sec_th) == 1.) mass_dot = 0.0; else { C_CENTROID(cell_center, cell, mixer); F_CENTROID(face_center, fa, tf); NV_VV(a, =, cell_center, -, face_center); distance = NV_MAG(a); F_AREA(A, fa, tf);//面積矢量// area = NV_MAG(A);//面積// wall_temper = F_T(fa, tf) - 273.15; Db = (4.709 + 0.196 * wall_temper + 0.006122 * wall_temper * wall_temper + 0.0004386 * wall_temper * wall_temper * wall_temper) / 1000; water_mass_fiction = C_YI(cell, first, 0); if (water_mass_fiction > Db) { mass_dot = -C_R(cell, first) * 0.256 * 0.0001 * (water_mass_fiction - Db) / distance * area / volume; dS[eqn] = 0; } } } } } end_c_loop(cell, mixer) return mass_dot * LAT_HT; |
回復(fù)此樓» 猜你喜歡
A區(qū)一本交叉課題組，低分調(diào)劑，招收機(jī)械電子信息通信等交叉方向
已經(jīng)有15人回復(fù)
-
材料學(xué)碩318求調(diào)劑
已經(jīng)有18人回復(fù)
-
江西師范大學(xué)王昆-盧章輝課題組招化學(xué)方向碩士研究生兩名
已經(jīng)有3人回復(fù)
-
沒上岸的看過來
已經(jīng)有5人回復(fù)
-
290求調(diào)劑
已經(jīng)有5人回復(fù)
-
材料與化工，291，求調(diào)劑
已經(jīng)有4人回復(fù)
-
301求調(diào)劑
已經(jīng)有4人回復(fù)
-
一志愿東北大學(xué)080500，英一數(shù)二343分求調(diào)劑
已經(jīng)有5人回復(fù)
-
274環(huán)境工程求調(diào)劑
已經(jīng)有6人回復(fù)
-
304求調(diào)劑
已經(jīng)有3人回復(fù)
mopsite
木蟲 (著名寫手)
- 應(yīng)助: 2 (幼兒園)
- 金幣: 8834.7
- 散金: 424
- 紅花: 2
- 帖子: 2575
- 在線: 615小時(shí)
- 蟲號: 172296
- 注冊: 2006-01-19
- 專業(yè): 傳熱傳質(zhì)學(xué)
2樓2021-05-12 14:28:33
18438609728
新蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 200.4
- 帖子: 63
- 在線: 3.4小時(shí)
- 蟲號: 24482155
- 注冊: 2020-11-02
- 專業(yè): 傳熱傳質(zhì)學(xué)
3樓2021-09-28 21:22:09
18438609728
新蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 200.4
- 帖子: 63
- 在線: 3.4小時(shí)
- 蟲號: 24482155
- 注冊: 2020-11-02
- 專業(yè): 傳熱傳質(zhì)學(xué)
4樓2021-09-28 21:22:25
18438609728
新蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 200.4
- 帖子: 63
- 在線: 3.4小時(shí)
- 蟲號: 24482155
- 注冊: 2020-11-02
- 專業(yè): 傳熱傳質(zhì)學(xué)
5樓2021-10-01 09:45:33
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 cfd菜雞 的主題更新
| 5 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 301求調(diào)劑 +4 | 李LJR 2026-03-04 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 一志愿天津大學(xué)085600 319分 材料與化工 金屬方向 求調(diào)劑 +6 | 青科11 2026-03-02 | 6/300 |
|
|
|
[考研]
085600 英一數(shù)二272求調(diào)劑
 5+6
|
vida_a 2026-03-01 | 47/2350 |
|
|
|
[考博] 26申博 +5 | north， 2026-02-28 | 5/250 |
|
|
|
[考研] 289求調(diào)劑 +8 | yang婷 2026-03-02 | 10/500 |
|
|
|
[考研] 0703 總分319求調(diào)劑 +4 | Xinyuu 2026-03-02 | 4/200 |
|
|
|
[考研] 276求調(diào)劑 +8 | 路lyh123 2026-02-28 | 10/500 |
|
|
|
[考研] 江蘇省農(nóng)科院招調(diào)劑1名 +5 | Qwertyuop 2026-03-01 | 5/250 |
|
|
|
[考研] 一志愿華中科技大學(xué)，化學(xué)專業(yè)344分，求調(diào)劑 +6 | 邢xing1 2026-03-02 | 6/300 |
|
|
|
[考研] 268求調(diào)劑 +6 | 好運(yùn)連綿不絕 2026-03-02 | 6/300 |
|
|
|
[考研] 299求調(diào)劑 +5 | kkcoco25 2026-03-02 | 9/450 |
|
|
|
[考研] 0857調(diào)劑 +6 | 一ll半 2026-02-28 | 8/400 |
|
|
|
[考研] 338求調(diào)劑 +5 | 18162027187 2026-03-02 | 6/300 |
|
|
|
[考研] 321求調(diào)劑一志愿東北林業(yè)大學(xué)材料與化工英二數(shù)二 +5 | 蟲蟲蟲蟲蟲7 2026-03-01 | 9/450 |
|
|
|
[考研] 0856材料調(diào)劑 +5 | 沿岸有貝殼OUC 2026-03-02 | 5/250 |
|
|
|
[考研] 0856求調(diào)劑285 +11 | 呂仔龍 2026-02-28 | 11/550 |
|
|
|
[考研] 一志愿山東大學(xué)材料與化工325求調(diào)劑 +5 | 半截的詩0927 2026-03-02 | 5/250 |
|
|
|
[基金申請] 此成果不能導(dǎo)入原因：元數(shù)據(jù)必填信息不完整，可 進(jìn)行補(bǔ)充。 +4 | Kittylucky 2026-03-02 | 5/250 |
|
|
|
[考研] 274求調(diào)劑 +3 | cgyzqwn 2026-03-01 | 7/350 |
|
|
|
[基金申請]
剛錄用，沒有期刊號，但是在線可看的論文可以放為代表作嗎
10+3
|
arang1 2026-03-01 | 3/150 |
|