应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » FLUENT & CFX » 用uds求解如下的问题为什么不收敛啊。谢谢！
181/2返回列表12下一页
查看: 3717  |  回复: 17
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
【悬赏金币】回答本帖问题，作者mopsite将赠送您 100 个金币

mopsite

木虫 (著名写手)

[求助] 用uds求解如下的问题为什么不收敛啊。谢谢！已有1人参与

用uds求解如下的问题为什么不收敛啊。谢谢！
比较简单的问题，2个方程，求解两个uds，两外2个uds储存数据和得出梯度。算例见附件。
当不加uds1的源项时，可以运算，但是一加uds1的源项，就发散。不知道问题出在哪，怀疑是在用uds求解梯度时不太对，查看可以发现梯度局部边界的值较大。
用uds求解如下的问题为什么不收敛啊。谢谢！
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 ]
回复此楼

» 本帖附件资源列表

  • 欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
    本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：libolin3@tal.com
  • 附件 1 : new.zip
    • 2013-12-17 17:47:22, 227.24 K

» 收录本帖的淘贴专辑推荐

传热传质

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼2013-12-17 17:48:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

gearboy78

禁虫 (小有名气)
本帖内容被屏蔽
3楼2013-12-19 08:01:17
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mopsite

木虫 (著名写手)
引用回帖:
12楼: Originally posted by longma135 at 2013-12-19 21:15:15
你好，我想请教下哈，在我编写的UDF程序中，连续性方程，动量方程，组分方程和能量方程都有源项编写，他们之间也会有联系，我都存储在UDM中，这就有个先后顺序问题，fluent怎么识别这个啊

各个方程有个求解顺序，你看看udf的帮助文档中，我记着有个各个方程的求解顺序
赞一下(1人) 回复此楼
13楼2013-12-19 22:26:25
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

mopsite

木虫 (著名写手)
顶顶 谁会啊
回复此楼
2楼2013-12-18 14:59:45
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mopsite

木虫 (著名写手)
引用回帖:
3楼: Originally posted by gearboy78 at 2013-12-19 08:01:17
不收敛的可能性太多，例如控制方程本身不对，松弛因子太大，编程不对等等。

谢谢 回复，您说的控制方程不对，是指哪个地方啊。 如果这个地方不对的话，就是根上不对了。
还请高手点拨点拨啊，叩谢了。。。
赞一下 回复此楼
4楼2013-12-19 08:19:06
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mopsite

木虫 (著名写手)
引用回帖:
3楼: Originally posted by gearboy78 at 2013-12-19 08:01:17
不收敛的可能性太多，例如控制方程本身不对，松弛因子太大，编程不对等等。

感谢您的提醒，控制方程可能是这样的，我修改下程序试下，谢谢。
用uds求解如下的问题为什么不收敛啊。谢谢！-1
eq.jpg
赞一下 回复此楼
5楼2013-12-19 09:46:18
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

longma135

金虫 (小有名气)

【答案】应助回帖

感谢参与，应助指数 +1
SV_UDS_I  这个写错了吧
赞一下 回复此楼
6楼2013-12-19 10:00:45
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mopsite

木虫 (著名写手)
引用回帖:
6楼: Originally posted by longma135 at 2013-12-19 10:00:45
SV_UDS_I  这个写错了吧

这个没错吧，您是指？。。。
赞一下 回复此楼
7楼2013-12-19 10:53:48
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

longma135

金虫 (小有名气)
引用回帖:
7楼: Originally posted by mopsite at 2013-12-19 10:53:48
这个没错吧，您是指？。。。...

不应该写成SV_UDSI(i)么
赞一下 回复此楼
8楼2013-12-19 11:06:01
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mopsite

木虫 (著名写手)
引用回帖:
8楼: Originally posted by longma135 at 2013-12-19 11:06:01
不应该写成SV_UDSI(i)么...

参考了那个温度梯度的写法，看这个链接 http://combust.hit.edu.cn:8080/f ... tml/udf/node238.htm
赞一下 回复此楼
9楼2013-12-19 11:12:55
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

longma135

金虫 (小有名气)
引用回帖:
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...

哦，这样，以后多交流哈
赞一下 回复此楼
10楼2013-12-19 11:19:09
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 mopsite 的主题更新
181/2返回列表12下一页
不应助 确定回帖应助 (注意：应助才可能被奖励，但不允许灌水，必须填写15个字符以上)
普通表情 龙 兔 虎 猫 SCI英文润色翻译-EnPapers美团队-先修改后付款
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[论文投稿] ACS Nano投稿后分配到副编辑手里12天了，能确定送审了吗？ +6 潇洒怡惜 2024-06-12 11/550 2024-06-18 01:48 by 玩耍绝缘体
[硕博家园] 硕博巨婴，也许才刚刚开始 +33 SNaiL1995 2024-06-12 85/4250 2024-06-18 00:06 by 博士求上上岸
[公派出国] 去英国的小伙伴儿都在哪儿租的房子呀？ +5 65syn 2024-06-14 8/400 2024-06-17 19:57 by 爱吃菜的虫
[论文投稿] 论文已接收，但发现修改稿传的是旧版该怎么办？ +3 zsq54321 2024-06-15 6/300 2024-06-17 19:45 by wjykycg
[基金申请] 博后面上今天有bug可以看到是否资助？ +22 lyfbangong 2024-06-12 33/1650 2024-06-17 16:17 by carolloo
[教师之家] 请问事业编制和年薪制冲突吗？ +10 ZHONGWU_U 2024-06-14 10/500 2024-06-17 14:44 by 周周520
[博后之家] 山东大学(青岛)“天然药物生物智造”课题组 招聘“博士后”(年薪20.4-55.6万元) +5 第二种态度 2024-06-11 11/550 2024-06-17 10:08 by 小懂事k
[找工作] 江西双非一本和四川双一流高校如何选择？ 5+8 寒山敲钟 2024-06-12 25/1250 2024-06-16 22:05 by zhaojiang427
[教师之家] 每次骚扰女学生的都是院系领导，而不是普通教师，小编们要注意措辞正确 +9 zju2000 2024-06-15 11/550 2024-06-16 14:49 by appleapple2
[基金申请] 希望今年自己国自然面上项目和老婆青年项目能中！ +7 恐龙爸爸 2024-06-14 7/350 2024-06-16 14:48 by redfish105
[基金申请] 博士后创新人才支持计划公示 +9 aishida144 2024-06-14 15/750 2024-06-16 09:52 by msjy
[论文投稿] 投稿被一个审稿人恶意评审了怎么样？ +5 1chen 2024-06-14 7/350 2024-06-15 23:15 by xy66xy
[基金申请] 为什么我的博后基金还在流动站审核中？不会是学院给我卡了吧？ +14 王凯12 2024-06-13 26/1300 2024-06-15 15:22 by 好人与坏人
[基金申请] 工材E口JQ有消息了吗 +4 babyduck 2024-06-11 4/200 2024-06-14 17:23 by firepick
[有机交流] ππ堆积会发生在有机溶剂中吗 5+3 zibuyu0420 2024-06-13 4/200 2024-06-14 14:17 by 小肉干
[基金申请] 75批博后基金 +10 kyukitu 2024-06-13 13/650 2024-06-14 10:31 by kyukitu
[论文投稿] 文章proof要求使用机构的邮箱 5+3 不可不信缘 2024-06-11 11/550 2024-06-14 07:00 by 3001160025
[论文投稿] 投稿后发现其他作者的邮箱填错了该怎么办呀 10+4 在飞的猪 2024-06-13 6/300 2024-06-14 04:45 by 小虫子咔咔
[硕博家园] 科研求助 +5 杲www 2024-06-12 6/300 2024-06-13 16:16 by 姓李名明
[论文投稿] 摩擦磨损论文投稿 +3 jmysan 2024-06-12 3/150 2024-06-13 08:36 by 莱茵润色
信息提示
关闭
请填处理意见
关闭