版塊導(dǎo)航: 正在加載中...

登錄注冊

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

24小時熱門版塊排行榜

返回列表

tsv249052

新蟲 (著名寫手)

應(yīng)助: 0 (幼兒園)
金幣: 661.4
帖子: 1164
在線: 2小時
蟲號: 20775031
注冊: 2020-02-03
性別: GG
專業(yè): 多相流熱物理學(xué)

[求助] 求助，UDS源項，一計算就發(fā)散是什么問題已有1人參與

求助，我自己編的udf程序，使用了define adjust和uds源項，在fluent中編譯build和加載load都沒有問題，但是計算一開始就顯示uds-0發(fā)散了，不知道是什么原因，部分代碼如下：
#include "udf.h"

#define PI 3.14159 /* 定義π*/
#define K 1.38e-23 /*定義玻爾茲曼常數(shù) */
#define M 2.99e-26 /*定義水分子質(zhì)量*/
#define R 461.5 /*定義水蒸氣氣體常數(shù)*/
#define f_v 0.4 /*相對濕度*/

real S1, G, J, pv1, ps1, T, P, sur_tension, density_c, density_v;/*不飽和度，吉布斯焓，成核率，水蒸氣分壓力，水蒸氣飽和壓力，溫度，壓力*/

/*計算成核率，液滴生長速度，臨界半徑*/
DEFINE_ADJUST(chenghelv, d)
{
Thread *t;
cell_t c;

thread_loop_c (t,d)
{
begin_c_loop_all(c,t)
{
T = C_T(c,t);
P = C_P(c,t);
ps1 = -2.88934 + 0.37966 * T - 0.00925 * pow(T, 2) + 1.54e-4 * pow(T, 3);/*飽和壓力擬合*/
density_v = -0.00134 + 0.00136 * T - 3.27e-5 * pow(T, 2) + 7.76e-7 * pow(T, 3);/*蒸氣密度擬合*/
density_c = 1000.07904 + 0.01231 * T - 0.00586 * pow(T,2) + 1.6e-5 * pow(T,3);/*水密度擬合*/
sur_tension = 0.07566 - 1.41e-4 * T - 2.68e-7 * pow(T,2);/*表面張力擬合*/
S1 = P / ps1;
pv1 = P * f_v;
G = 16 / 3 * PI * pow((M / (density_c * log(S1) * K * T)),2) * pow(sur_tension, 3); /*計算吉布斯自由焓*/
C_UDMI(c,t,0) = sqrt(2 / PI * sur_tension) * pow(M ,-3/2) * pow(density_v, 2) / density_c * exp( -G / K * T);/*成核率*/
end_c_loop_all(c,t)
}
}
}

/*第一個輸送方程Q0源項*/
DEFINE_SOURCE(q0_source,c,t,dS,eqn)
{
return C_UDMI(c,t,0);
}

/*第二個輸送方程Q1源項*/
DEFINE_SOURCE(q1_source,c,t,dS,eqn)
{
return C_UDMI(c,t,0) * pow(C_UDMI(c,t,2),2) + C_R(c,t) * C_UDSI(c,t,0) * C_UDMI(c,t,1);
}

/*第三個輸送方程Q2源項*/
DEFINE_SOURCE(q2_source,c,t,dS,eqn)
{
return C_UDMI(c,t,0) * pow(C_UDMI(c,t,2),2) + 2 * C_R(c,t) * C_UDSI(c,t,1) * C_UDMI(c,t,1);
}

/*第四個輸送方程g的源項*/
DEFINE_SOURCE(g_source,c,t,dS,eqn)
{
return 4 * PI * density_c * (1/3 * C_UDMI(c,t,0) * pow(C_UDMI(c,t,2),3) + C_R(c,t) * C_UDSI(c,t,2) * C_UDMI(c,t,1));
}

/*能量方程源項*/
DEFINE_SOURCE(energy_source,c,t,dS,eqn)
{
real dt, temp, L;
dt = RP_Get_Real("physical-time-step"

;
temp = C_T(c,t);
L = 2502 - 2.42169 * temp;
return L * C_UDSI(c,t,3) / dt;
}

回復(fù)此樓

» 猜你喜歡

西安電子科學(xué)大學(xué)杭州研究院劉麗香教授招收智能多模態(tài)傳感器和微型儲能器件方向博士已經(jīng)有13人回復(fù)
重慶交通大學(xué)光子學(xué)微結(jié)構(gòu)與器件課題組招收2026年碩士研究生信息已經(jīng)有1人回復(fù)
物理學(xué)I論文潤色/翻譯怎么收費? 已經(jīng)有156人回復(fù)
一志愿鄭大材料學(xué)碩298分，求調(diào)劑已經(jīng)有6人回復(fù)
尋合作：應(yīng)力腐蝕多尺度模擬已經(jīng)有3人回復(fù)
考研交流已經(jīng)有0人回復(fù)
【新加坡】納米電子器件項目組有“聯(lián)合培養(yǎng)博士生”名額已經(jīng)有0人回復(fù)
【新加坡】納米電子器件項目組有“聯(lián)合培養(yǎng)博士生”名額已經(jīng)有0人回復(fù)

» 本主題相關(guān)商家推薦: (我也要在這里推廣)

1樓 2020-03-17 19:08:52

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

Diligent雨

新蟲 (初入文壇)

應(yīng)助: 3 (幼兒園)
金幣: 37.6
帖子: 28
在線: 16.2小時
蟲號: 5052080
注冊: 2016-09-23
性別: GG
專業(yè): 金屬材料表面科學(xué)與工程

【答案】應(yīng)助回帖

把udf簡化嘗試，運用排除法試試；可能是fluent內(nèi)部設(shè)置問題，查看與uds0所有的設(shè)置；減小松弛因子。希望有幫助

贊一下

回復(fù)此樓

焊接fluent熔池

2樓2020-03-23 22:01:49

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

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