mopsite

木蟲 (著名寫手)

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

[求助] 用uds求解如下的問題為什么不收斂啊。謝謝！已有1人參與

用uds求解如下的問題為什么不收斂啊。謝謝！
比較簡單的問題，2個方程，求解兩個uds，兩外2個uds儲存數(shù)據(jù)和得出梯度。算例見附件。
當(dāng)不加uds1的源項時，可以運算，但是一加uds1的源項，就發(fā)散。不知道問題出在哪，懷疑是在用uds求解梯度時不太對，查看可以發(fā)現(xiàn)梯度局部邊界的值較大。

3.jpg#include "udf.h"
#define B_ION 1.5e-4 /* Ion mobility - positive ion estimate */
#define EPS_0 8.85418717e-12

DEFINE_ADJUST(Calc_E_xy,d)
{
Thread *t;
cell_t c;
face_t f;
Thread *t0;
cell_t c0;
real dr0[ND_ND],dr1[ND_ND], xf[ND_ND], dy;

/* Do nothing if gradient isn't allocated yet. */
if (!Data_Valid_P()) return;
thread_loop_c(t,d)
{
if (NULL != THREAD_STORAGE(t,SV_UDS_I(0)) && NULL != T_STORAGE_R_NV(t,SV_UDSI_G(0)))
{
begin_c_loop(c,t)
{
C_UDSI(c,t,2) = -1.*C_UDSI_G(c,t,0)[0];
C_UDSI(c,t,3) = -1.*C_UDSI_G(c,t,0)[1];
}
end_c_loop(c,t)
}
}
thread_loop_f(t,d)
{
if (NULL != THREAD_STORAGE(t,SV_UDS_I(3)) && NULL != T_STORAGE_R_NV(t->t0,SV_UDSI_G(3)))
{
if (BOUNDARY_FACE_THREAD_P(t))
{
begin_f_loop(f,t)
{
c0 = F_C0(f,t);
t0 = THREAD_T0(t);
F_CENTROID(xf,f,t);
C_CENTROID(dr0,c0,t0);
NV_VV(dr1,=,xf,-,dr0);
F_UDSI(f,t,2) = C_UDSI(c0,t0,2) + (C_UDSI_G(c0,t0,2)[0]*dr1[0]);
F_UDSI(f,t,3) = C_UDSI(c0,t0,3) + (C_UDSI_G(c0,t0,3)[1]*dr1[1]);
}
end_f_loop(f,t)
}
}
}
}

/*******************************************************************/
/* UDF_source_phi_rho_11.c for specifying the phi and rho source terms for uds-0 */
/*******************************************************************/
DEFINE_SOURCE(phi_source_temp,c,t,dS,eqn)
{
real phi_source;
dS[eqn] = 0.0;
phi_source = C_UDSI(c,t,1)/EPS_0;
return phi_source;
}
DEFINE_SOURCE(rho_source_temp,c,t,dS,eqn)
{
real rho_source;
real xc[ND_ND];
real dens;
dS[eqn] = 0;
C_CENTROID(xc,c,t);
dens=C_R(c,t);
rho_source = C_UDSI(c,t,1)*B_ION*(C_UDSI_G(c,t,2)[0] + C_UDSI_G(c,t,3)[1]);
rho_source = -rho_source;
return rho_source;
}

DEFINE_UDS_FLUX(rho_ion_flux,f,t,i)
{
cell_t c0, c1 = -1;
Thread *t0, *t1 = NULL;
real NV_VEC(psi_vec), NV_VEC(A), flux = 0.0;
c0 = F_C0(f,t);
t0 = F_C0_THREAD(f,t);
F_AREA(A, f, t);

if (BOUNDARY_FACE_THREAD_P(t))
{

NV_DS(psi_vec, =, F_U(f,t), F_V(f,t), F_W(f,t), *, 1);
flux = NV_DOT(psi_vec, A); /* flux through Face */
}
else
{
c1 = F_C1(f,t); /* Get cell on other side of face */
t1 = F_C1_THREAD(f,t);
NV_DS(psi_vec, =, C_U(c0,t0),C_V(c0,t0),C_W(c0,t0),*,1.0);
NV_DS(psi_vec, +=, C_U(c1,t1),C_V(c1,t1),C_W(c1,t1),*,1.0);
flux = NV_DOT(psi_vec, A)/2.0; /* Average flux through face */
}
return flux;
}

[ Last edited by mopsite on 2013-12-17 at 18:38 ]

回復(fù)此樓

» 本帖附件資源列表

歡迎監(jiān)督和反饋：小木蟲僅提供交流平臺，不對該內(nèi)容負(fù)責(zé)。
本內(nèi)容由用戶自主發(fā)布，如果其內(nèi)容涉及到知識產(chǎn)權(quán)問題，其責(zé)任在于用戶本人，如對版權(quán)有異議，請聯(lián)系郵箱：xiaomuchong@tal.com
附件 1 : new.zip

2013-12-17 17:47:22, 227.24 K

1樓 2013-12-17 17:48:09

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

回帖支持 ( 顯示支持度最高的前 50 名 )

gearboy78

禁蟲 (小有名氣)

本帖內(nèi)容被屏蔽

3樓2013-12-19 08:01:17

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

mopsite

木蟲 (著名寫手)

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

引用回帖:

12樓: Originally posted by longma135 at 2013-12-19 21:15:15
你好，我想請教下哈，在我編寫的UDF程序中，連續(xù)性方程，動量方程，組分方程和能量方程都有源項編寫，他們之間也會有聯(lián)系，我都存儲在UDM中，這就有個先后順序問題，fluent怎么識別這個啊

各個方程有個求解順序，你看看udf的幫助文檔中，我記著有個各個方程的求解順序

贊一下(1人)

回復(fù)此樓

13樓2013-12-19 22:26:25

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

普通回帖

mopsite

木蟲 (著名寫手)

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

頂頂誰會啊

回復(fù)此樓

2樓2013-12-18 14:59:45

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

mopsite

木蟲 (著名寫手)

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

引用回帖:

3樓: Originally posted by gearboy78 at 2013-12-19 08:01:17
不收斂的可能性太多，例如控制方程本身不對，松弛因子太大，編程不對等等。

謝謝回復(fù)，您說的控制方程不對，是指哪個地方啊。如果這個地方不對的話，就是根上不對了。
還請高手點撥點撥啊，叩謝了。。。

贊一下

回復(fù)此樓

4樓2013-12-19 08:19:06

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

mopsite

木蟲 (著名寫手)

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

引用回帖:

3樓: Originally posted by gearboy78 at 2013-12-19 08:01:17
不收斂的可能性太多，例如控制方程本身不對，松弛因子太大，編程不對等等。

感謝您的提醒，控制方程可能是這樣的，我修改下程序試下，謝謝。

eq.jpg

贊一下

回復(fù)此樓

5樓2013-12-19 09:46:18

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

longma135

金蟲 (小有名氣)

應(yīng)助: 5 (幼兒園)
金幣: 1186.6
散金: 15
帖子: 199
在線: 207小時
蟲號: 2533271
注冊: 2013-07-05
性別: GG
專業(yè): 能源化工

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

感謝參與，應(yīng)助指數(shù) +1

SV_UDS_I 這個寫錯了吧

贊一下

回復(fù)此樓

6樓2013-12-19 10:00:45

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

mopsite

木蟲 (著名寫手)

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

引用回帖:

6樓: Originally posted by longma135 at 2013-12-19 10:00:45
SV_UDS_I 這個寫錯了吧

這個沒錯吧，您是指？。。。

贊一下

回復(fù)此樓

7樓2013-12-19 10:53:48

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

longma135

金蟲 (小有名氣)

應(yīng)助: 5 (幼兒園)
金幣: 1186.6
散金: 15
帖子: 199
在線: 207小時
蟲號: 2533271
注冊: 2013-07-05
性別: GG
專業(yè): 能源化工

引用回帖:

7樓: Originally posted by mopsite at 2013-12-19 10:53:48
這個沒錯吧，您是指？。。。...

不應(yīng)該寫成SV_UDSI(i)么

贊一下

回復(fù)此樓

8樓2013-12-19 11:06:01

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

mopsite

木蟲 (著名寫手)

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

引用回帖:

8樓: Originally posted by longma135 at 2013-12-19 11:06:01
不應(yīng)該寫成SV_UDSI(i)么...

參考了那個溫度梯度的寫法，看這個鏈接 http://combust.hit.edu.cn:8080/f ... tml/udf/node238.htm

贊一下

回復(fù)此樓

9樓2013-12-19 11:12:55

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

longma135

金蟲 (小有名氣)

應(yīng)助: 5 (幼兒園)
金幣: 1186.6
散金: 15
帖子: 199
在線: 207小時
蟲號: 2533271
注冊: 2013-07-05
性別: GG
專業(yè): 能源化工

引用回帖:

9樓: Originally posted by mopsite at 2013-12-19 11:12:55
參考了那個溫度梯度的寫法，看這個鏈接 http://combust.hit.edu.cn:8080/fluent/Fluent60_help/html/udf/node238.htm...

哦，這樣，以后多交流哈

贊一下

回復(fù)此樓

10樓2013-12-19 11:19:09

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

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

[考研] 354求調(diào)劑 +4	Tyoumou 2026-03-18	7/350	2026-03-18 21:45 by Tyoumou
[考研] 328求調(diào)劑，英語六級551，有科研經(jīng)歷 +3	生物工程調(diào)劑 2026-03-17	7/350	2026-03-18 20:41 by Wangjingyue
[教師之家] 焦慮 +8	水冰月月野兔 2026-03-13	12/600	2026-03-18 15:27 by 咪嗚喵嗚
[考研] 0854，計算機(jī)類招收調(diào)劑 +3	胡辣湯放糖 2026-03-15	6/300	2026-03-18 12:09 by 上岸上岸……..
[考研] 280求調(diào)劑 +6	咕嚕曉曉 2026-03-18	7/350	2026-03-18 11:25 by 無際的草原
[考研] 303求調(diào)劑 +4	睿08 2026-03-17	6/300	2026-03-18 11:01 by Iveryant
[考研] 0703化學(xué)求調(diào)劑總分331 +3	ZY-05 2026-03-13	3/150	2026-03-18 10:58 by macy2011
[考研] 268求調(diào)劑 +6	簡單點0 2026-03-17	6/300	2026-03-18 09:04 by 無際的草原
[考研] 一志愿蘇州大學(xué)材料工程（085601）專碩有科研經(jīng)歷三項國獎兩個實用型專利一項省級立項 +6	大火山小火山 2026-03-16	8/400	2026-03-17 15:05 by 無懈可擊111
[考研] 302求調(diào)劑 +4	小賈同學(xué)123 2026-03-15	8/400	2026-03-17 10:33 by 小賈同學(xué)123
[考研] 285化工學(xué)碩求調(diào)劑（081700） +9	柴郡貓_ 2026-03-12	9/450	2026-03-17 10:18 by Sammy2
[基金申請] 國自科面上基金字體 +6	iwuli 2026-03-12	7/350	2026-03-16 21:18 by sculhf
[考研] 333求調(diào)劑 +3	文思客 2026-03-16	7/350	2026-03-16 18:21 by 文思客
[考研] 0703一志愿211 285分求調(diào)劑 +5	ly3471z 2026-03-13	5/250	2026-03-16 16:16 by 哦哦123
[碩博家園] 085600 260分求調(diào)劑 +3	天空還下雨么 2026-03-13	5/250	2026-03-13 18:46 by 天空還下雨么
[考研] 274求調(diào)劑 +3	S.H1 2026-03-12	3/150	2026-03-13 15:15 by JourneyLucky
[考研] 求調(diào)劑 +3	程雨杭 2026-03-12	3/150	2026-03-13 15:06 by JourneyLucky
[考研] 308求調(diào)劑 +3	是Lupa啊 2026-03-12	3/150	2026-03-13 14:30 by 求調(diào)劑zz
[考研] 0817化學(xué)工程與技術(shù)考研312分調(diào)劑 +3	T123 tt 2026-03-12	3/150	2026-03-13 10:49 by houyaoxu
[考研] 333求調(diào)劑 +3	152697 2026-03-12	4/200	2026-03-13 07:08 by Iveryant