 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 4493 | 回復(fù): 17 | |||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | |||
mopsite木蟲 (著名寫手)
|
[求助]
用uds求解如下的問(wèn)題為什么不收斂啊。謝謝！ 已有1人參與
|
||
|
用uds求解如下的問(wèn)題為什么不收斂啊。謝謝！ 比較簡(jiǎn)單的問(wèn)題，2個(gè)方程，求解兩個(gè)uds，兩外2個(gè)uds儲(chǔ)存數(shù)據(jù)和得出梯度。算例見附件。 當(dāng)不加uds1的源項(xiàng)時(shí)，可以運(yùn)算，但是一加uds1的源項(xiàng)，就發(fā)散。不知道問(wèn)題出在哪，懷疑是在用uds求解梯度時(shí)不太對(duì)，查看可以發(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)督和反饋：小木蟲僅提供交流平臺(tái)，不對(duì)該內(nèi)容負(fù)責(zé)。
本內(nèi)容由用戶自主發(fā)布，如果其內(nèi)容涉及到知識(shí)產(chǎn)權(quán)問(wèn)題，其責(zé)任在于用戶本人，如對(duì)版權(quán)有異議，請(qǐng)聯(lián)系郵箱：xiaomuchong@tal.com - 附件 1 : new.zip
2013-12-17 17:47:22, 227.24 K
» 收錄本帖的淘帖專輯推薦
 傳熱傳質(zhì) |
» 猜你喜歡
求調(diào)劑
已經(jīng)有3人回復(fù)
-
265求調(diào)劑
已經(jīng)有4人回復(fù)
-
085700資源與環(huán)境308求調(diào)劑
已經(jīng)有6人回復(fù)
-
一志愿吉林大學(xué)材料學(xué)碩321求調(diào)劑
已經(jīng)有12人回復(fù)
-
286分人工智能專業(yè)請(qǐng)求調(diào)劑愿意跨考！
已經(jīng)有3人回復(fù)
-
329求調(diào)劑
已經(jīng)有5人回復(fù)
-
申請(qǐng)回稿延期一個(gè)月，編輯同意了。但系統(tǒng)上的時(shí)間沒變，給編輯又寫郵件了，沒回復(fù)
已經(jīng)有4人回復(fù)
-
材料學(xué)碩318求調(diào)劑
已經(jīng)有5人回復(fù)
-
一志愿中國(guó)海洋大學(xué)，生物學(xué)，301分，求調(diào)劑
已經(jīng)有6人回復(fù)
-
081700化工學(xué)碩調(diào)劑
已經(jīng)有3人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- fluent 仿真溫度出現(xiàn)負(fù)值
已經(jīng)有9人回復(fù)
- 求高手幫忙看一段英文，涉及到數(shù)據(jù)分析問(wèn)題，如能給出說(shuō)得通的全文翻譯另有50BB相贈(zèng)
已經(jīng)有5人回復(fù)
- PCG求解器迭代不收斂
已經(jīng)有4人回復(fù)
- ANSYS添加化學(xué)反應(yīng)交流
已經(jīng)有5人回復(fù)
- 用fluent模擬時(shí)，選擇非穩(wěn)態(tài)求解器，收斂以后還繼續(xù)算會(huì)不會(huì)對(duì)結(jié)果有影響？
已經(jīng)有7人回復(fù)
- fluent求解UDS過(guò)程中，出現(xiàn)少數(shù)不合理的值..
已經(jīng)有14人回復(fù)
- fluent 6.3 13.0版本中多相流的比較
已經(jīng)有20人回復(fù)
- UDF中DEFINE_ADJUST宏中使用C_UDSI_G
已經(jīng)有36人回復(fù)
mopsite
木蟲 (著名寫手)
- 應(yīng)助: 2 (幼兒園)
- 金幣: 8891.7
- 散金: 424
- 紅花: 2
- 帖子: 2586
- 在線: 615小時(shí)
- 蟲號(hào): 172296
- 注冊(cè): 2006-01-19
- 專業(yè): 傳熱傳質(zhì)學(xué)
9樓2013-12-19 11:12:55
|
本帖內(nèi)容被屏蔽 |
3樓2013-12-19 08:01:17
mopsite
木蟲 (著名寫手)
- 應(yīng)助: 2 (幼兒園)
- 金幣: 8891.7
- 散金: 424
- 紅花: 2
- 帖子: 2586
- 在線: 615小時(shí)
- 蟲號(hào): 172296
- 注冊(cè): 2006-01-19
- 專業(yè): 傳熱傳質(zhì)學(xué)
4樓2013-12-19 08:19:06
mopsite
木蟲 (著名寫手)
- 應(yīng)助: 2 (幼兒園)
- 金幣: 8891.7
- 散金: 424
- 紅花: 2
- 帖子: 2586
- 在線: 615小時(shí)
- 蟲號(hào): 172296
- 注冊(cè): 2006-01-19
- 專業(yè): 傳熱傳質(zhì)學(xué)
5樓2013-12-19 09:46:18
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿蘇州大學(xué)材料求調(diào)劑，總分315（英一） +3 | sbdksD 2026-03-19 | 3/150 |
|
|---|---|---|---|---|
|
|
[考研] 0703化學(xué)調(diào)劑 +10 | 妮妮ninicgb 2026-03-15 | 14/700 |
|
|
|
[考研] 本人考085602 化學(xué)工程 專碩 +17 | 不知道叫什么！ 2026-03-15 | 19/950 |
|
|
|
[考研] 一志愿福大288有機(jī)化學(xué)，求調(diào)劑 +3 | 小木蟲200408204 2026-03-18 | 3/150 |
|
|
|
[考研] 材料考研調(diào)劑 +3 | xwt。 2026-03-19 | 3/150 |
|
|
[考研] 332求調(diào)劑 +3 | ydfyh 2026-03-17 | 3/150 |
|
|
|
[考研] 354求調(diào)劑 +4 | Tyoumou 2026-03-18 | 7/350 |
|
|
|
[考研] 302求調(diào)劑 +10 | 呼呼呼。。。。 2026-03-17 | 10/500 |
|
|
|
[考研] 299求調(diào)劑 +5 | △小透明* 2026-03-17 | 5/250 |
|
|
|
[考研] 環(huán)境工程調(diào)劑 +8 | 大可digkids 2026-03-16 | 8/400 |
|
|
|
[考研] 265求調(diào)劑 +3 | 梁梁校校 2026-03-17 | 3/150 |
|
|
|
[考研] 268求調(diào)劑 +6 | 簡(jiǎn)單點(diǎn)0 2026-03-17 | 6/300 |
|
|
|
[考研] [導(dǎo)師推薦]西南科技大學(xué)國(guó)防/材料導(dǎo)師推薦 +3 | 尖角小荷 2026-03-16 | 6/300 |
|
|
|
[考研] 318求調(diào)劑 +3 | Yanyali 2026-03-15 | 3/150 |
|
|
|
[考研] 304求調(diào)劑 +3 | 曼殊2266 2026-03-14 | 3/150 |
|
|
|
[考研] 327求調(diào)劑 +6 | 拾光任染 2026-03-15 | 11/550 |
|
|
|
[考研] 求老師收留調(diào)劑 +4 | jiang姜66 2026-03-14 | 5/250 |
|
|
|
[考研] 289求調(diào)劑 +4 | 這么名字咋樣 2026-03-14 | 6/300 |
|
|
|
[考研] 復(fù)試調(diào)劑 +3 | 呼呼？~+123456 2026-03-14 | 3/150 |
|
|
|
[考研] 一志愿哈工大材料324分求調(diào)劑 +5 | 閆旭東 2026-03-14 | 5/250 |
|