版塊導航: 正在加載中...

登錄注冊

應《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進行實名認證將不得使用互聯(lián)網(wǎng)跟帖服務。為保障您的帳號能夠正常使用，請盡快對帳號進行手機號驗證，感謝您的理解與支持！

24小時熱門版塊排行榜

返回列表

【懸賞金幣】回答本帖問題，作者cfd菜雞將贈送您 40 個金幣

cfd菜雞

新蟲 (初入文壇)

應助: 0 (幼兒園)
金幣: 12
帖子: 15
在線: 6.9小時
蟲號: 26191183
注冊: 2021-04-13
專業(yè): 傳熱傳質(zhì)學

[求助] 這是我編寫濕空氣冷凝的udf導入fluent后初始化就無法完成，是不是死循環(huán)，求大神指導

#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中心坐標//
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);//水蒸氣質(zhì)量分數(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;

回復此樓

» 猜你喜歡

歡迎211本科同學，過A區(qū)國家線，A區(qū)非偏遠一本，交叉學科課題組已經(jīng)有24人回復
新疆大學地質(zhì)與礦業(yè)工程學院招生已經(jīng)有10人回復
招調(diào)劑已經(jīng)有11人回復
0856材料與化工求調(diào)劑！已經(jīng)有11人回復
材料與化工304求調(diào)劑已經(jīng)有10人回復
A區(qū)一本交叉課題組，低分調(diào)劑，招收機械電子信息通信等交叉方向已經(jīng)有31人回復
中科大材料299求調(diào)劑已經(jīng)有7人回復
301求調(diào)劑已經(jīng)有3人回復
材料考研339求調(diào)劑已經(jīng)有3人回復
材料專碩290求調(diào)劑已經(jīng)有5人回復

1樓 2021-04-27 16:39:49

已閱回復此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

mopsite

木蟲 (著名寫手)

應助: 2 (幼兒園)
金幣: 8839.7
散金: 424
紅花: 2
帖子: 2576
在線: 615小時
蟲號: 172296
注冊: 2006-01-19
專業(yè): 傳熱傳質(zhì)學

顯示啥問題啊，是僅計算壁面的冷凝問題嘛，源項里咋沒有時間相關(guān)的項。

贊一下(1人)

回復此樓

2樓2021-05-12 14:28:33

已閱回復此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

18438609728

新蟲 (小有名氣)

應助: 0 (幼兒園)
金幣: 200.4
帖子: 63
在線: 3.4小時
蟲號: 24482155
注冊: 2020-11-02
專業(yè): 傳熱傳質(zhì)學

能發(fā)一份usf

發(fā)自小木蟲Android客戶端

回復此樓

3樓2021-09-28 21:22:09

已閱回復此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

18438609728

新蟲 (小有名氣)

應助: 0 (幼兒園)
金幣: 200.4
帖子: 63
在線: 3.4小時
蟲號: 24482155
注冊: 2020-11-02
專業(yè): 傳熱傳質(zhì)學

能發(fā)一份udf嗎，咱們的研究方向接近

發(fā)自小木蟲Android客戶端

贊一下

回復此樓

4樓2021-09-28 21:22:25

已閱回復此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

18438609728

新蟲 (小有名氣)

應助: 0 (幼兒園)
金幣: 200.4
帖子: 63
在線: 3.4小時
蟲號: 24482155
注冊: 2020-11-02
專業(yè): 傳熱傳質(zhì)學

引用回帖:

2樓: Originally posted by mopsite at 2021-05-12 14:28:33
顯示啥問題啊，是僅計算壁面的冷凝問題嘛，源項里咋沒有時間相關(guān)的項。

我是做水蒸氣冷凝的，但是我模擬出來的熱流密度一直太小。你有沒有相應的udf啊，給我分享一下

發(fā)自小木蟲Android客戶端

贊一下

回復此樓

5樓2021-10-01 09:45:33

已閱回復此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 cfd菜雞的主題更新

返回列表

不應助 確定回帖應助 (注意：應助才可能被獎勵，但不允許灌水，必須填寫15個字符以上)

普通表情龍兔虎貓

最具人氣熱帖推薦 [查看全部]		作者	回/看	最后發(fā)表

[考研] 304求調(diào)劑 +3	52hz~~ 2026-03-05	4/200	2026-03-06 00:27 by wutongshun
[考研] 328求調(diào)劑 +6	vuzhdkfjkx 2026-03-04	7/350	2026-03-05 23:28 by Leeding1356
[考研] 271求調(diào)劑 +7	月色c 2026-03-05	8/400	2026-03-05 19:51 by wangjihu
[考研] 0703 東華大學理學化學均過 a 線，總分 281 求調(diào)劑 +4	香香 Lu 2026-03-03	8/400	2026-03-05 17:36 by 香香 Lu
[考研] 紡織、生物、化學、材料等專業(yè) +3	Eember. 2026-03-05	7/350	2026-03-05 16:35 by Eember.
[考研] 材料類考研調(diào)劑 +7	gemmgemm 2026-03-01	8/400	2026-03-05 10:51 by 早安1234
[考研] 302材料工程求調(diào)劑 +7	Doleres 2026-03-01	8/400	2026-03-05 09:59 by oxidpl
[考研] 0703化學學碩理工科均可不區(qū)分研究方向總分279求調(diào)劑 +7	1一11 2026-03-03	7/350	2026-03-05 09:23 by xwxstudy
[考研] 調(diào)劑材料學碩 +6	詞凝Y 2026-03-02	6/300	2026-03-05 07:10 by 昆工985
[考研] 求調(diào)劑院校 +6	云朵452 2026-03-02	12/600	2026-03-04 23:17 by 云朵452
[考研] 293求調(diào)劑 +3	是樂渝哇 2026-03-04	3/150	2026-03-04 23:12 by wutongshun
[考研] 286 +6	ksncj 2026-03-04	6/300	2026-03-04 20:49 by lature00
[考研] 266求調(diào)劑 +7	哇塞王帥 2026-03-03	7/350	2026-03-04 20:30 by 一切OK
[考研] 一志愿鄭州大學，學碩，物理化學， 333求調(diào)劑 +5	李魔女斗篷 2026-03-04	5/250	2026-03-04 18:33 by houyaoxu
[考研] 本科太原理工采礦工程，求調(diào)劑 +3	onlx 2026-03-01	3/150	2026-03-04 15:57 by Stephen_ym
[考研] 一志愿314求調(diào)劑 +7	202111120625 2026-03-03	7/350	2026-03-04 15:56 by zhukairuo
[考研] 347分電子信息專碩求調(diào)劑 +3	咕嘟咕嘟1 2026-03-03	3/150	2026-03-04 12:32 by 熱情沙漠
[考研] 化工335求調(diào)劑 +5	摸摸貓貓頭 2026-03-02	5/250	2026-03-03 09:10 by houyaoxu
[考研] 0856調(diào)劑 +10	劉夢微 2026-02-28	10/500	2026-03-02 23:42 by hcy618
[考研] 一志愿鄭大材料學碩298分，求調(diào)劑 +6	wsl111 2026-03-01	6/300	2026-03-02 11:00 by ydudjddnd

亭亭五月天在线观看,亭亭五月天在线观看,国产最新av一区二区,国产 高清 中文字幕,99re热久久亚洲综合精品成人,熟妇 一区二区三区,一级做a爰片性色毛片武则天,美女的骚穴视频播放,国产美女午夜免费视频

24小時熱門版塊排行榜

cfd菜雞

[求助] 這是我編寫濕空氣冷凝的udf導入fluent后初始化就無法完成，是不是死循環(huán)，求大神指導

» 猜你喜歡

mopsite

18438609728

18438609728

18438609728

亭亭五月天在线观看,亭亭五月天在线观看,国产最新av一区二区,国产高清中文字幕,99re热久久亚洲综合精品成人,熟妇一区二区三区,一级做a爰片性色毛片武则天,美女的骚穴视频播放,国产美女午夜免费视频