24小时热门版块排行榜    

查看: 4988  |  回复: 16
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

b07988

金虫 (正式写手)

星期六可以休息

[求助] 能量源项中有温度对时间的微分,需要用uds吗? 已有1人参与

能量源项中有dT/dt,就是温度对时间的求导,请问是用
1.c_udsi_m1存储上一步时间步长的温度
2.采用define_uds_unsteady?
希望有好心人帮忙。
不胜感激!
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

先做最重要的事情,再做次重要的事情
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

47425804

木虫 (正式写手)


★ ★
xiegangmai: 金币+2, 鼓励讨论交流! 2012-05-17 00:01:09
DEFINE_ADJUST里面保存温度那一步,换成下面这样试试:
C_UDMI(c,t,0)=C_T(c,t);
DEFINE_SOURCE里的自然就变成:
dT=(C_T(c,t)-C_UDMI(c,t,0))/RP_Get_Real("physical-time-step";
我觉得用UDM比用UDS好理解。

» 本帖已获得的红花(最新10朵)

4楼2012-05-16 22:15:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

47425804

木虫 (正式写手)


★ ★ ★
xiegangmai: 金币+3, 谢谢参与! 2012-05-17 23:09:28
引用回帖:
6楼: Originally posted by b07988 at 2012-05-17 09:24:09:
按你的思路dT会恒等于0.
应该是C_UDMI(c,t,0)=C_T_M1(c,t);吧。但是M1的话这样的话会直接报错。
所以我用了C_UDSI_M1(c,t,0)

恩,你说的没错,这样确实会恒等于0。
但是C_T_M1(c,t)是没用的。
等于零是因为计算源项的时候温度还没有更新,所以两个温度值是相等的。
这个时候其实是需要两个UDM,分别保存上一步和再上一步的温度。
而且可以考虑全部在源项里完成,不用Adjust宏:
source=(C_UDMI(c,t,1)-C_UDMI(c,t,0))/TIME_STEP;
C_UDMI(c,t,0)=C_UDMI(c,t,1);
C_UDMI(c,t,1)=C_T(c,t);
return source;
当然,在刚开始计算的时候要注意下UDM的初始问题,防止出错。
从第三个时间步开始,就是正常的了。

» 本帖已获得的红花(最新10朵)

7楼2012-05-17 11:48:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

47425804

木虫 (正式写手)


★ ★
xiegangmai: 金币+2, 谢谢参与! 2012-05-17 22:59:44
引用回帖:
10楼: Originally posted by b07988 at 2012-05-17 16:45:53:
这个方法挺好的,只要一开始的用DEFINE_INIT,对
C_UDMI(c,t,0)=0.;
C_UDMI(c,t,1)=0.;
就可以了!
我试了下,没有什么问题。谢谢了!

我也用uds试了下,其实也没有什么问题。只要打开uds就可以了。:vi ...

用patch会比较方便,不用DEFINE_INIT。
而且,开始时候给0不是太合适,用个平均温度更好点。
不过也就三个时间步的事,影响很小。

不管UDM还是UDS,能用就是好东西,呵呵~~
11楼2012-05-17 17:26:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 b07988 的主题更新
信息提示
请填处理意见