24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 5178  |  回复: 16

b07988

金虫 (正式写手)

星期六可以休息

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

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

» 猜你喜欢

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

先做最重要的事情,再做次重要的事情
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

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

金虫 (正式写手)

星期六可以休息

引用回帖:
2楼: Originally posted by ya634984094 at 2012-05-16 15:27:20:
可以用也可以不要用~用的话用C_UDSI_M1和C_UDSI俩个值之差除以时间就行,不用的话用C_UDMI保存温度,再求~

其实我两种都采用了,但不知道是自己程序出问题还是什么,就是没有效果。
以下是部分程序,希望你能给点意见:
CODE:
DEFINE_ADJUST(store_m1,domain)
{
        Thread *t;
        cell_t c;
        real haha=1.;
        domain =Get_Domain(1);
        thread_loop_c(t,domain)
        {
                begin_c_loop(c,t)
                {
                        C_UDSI_M1(c,t,0)=C_T_M1(c,t);
                }
                end_c_loop(c,t)
        }       
}
然后在define_source里面加入
dT=(C_T(c,t)-C_UDSI_M1(c,t,0))/RP_Get_Real("physical-time-step");

但是没有效果,请帮我看看。
先做最重要的事情,再做次重要的事情
3楼2012-05-16 18:43:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

ya634984094

铜虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
b07988: 金币+5, ★★★很有帮助, 谢谢参与 2012-05-16 18:28:33
xiegangmai: 金币+1, 谢谢应助! 2012-05-17 00:00:56
可以用也可以不要用~用的话用C_UDSI_M1和C_UDSI俩个值之差除以时间就行,不用的话用C_UDMI保存温度,再求~
2楼2012-05-16 15:27:20
已阅   回复此楼   关注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的回帖

030301029

捐助贵宾 (正式写手)

【答案】应助回帖

★ ★ ★ ★
感谢参与,应助指数 +1
xiegangmai: 金币+1, 谢谢参与! 2012-05-17 00:01:19
b07988: 金币+3, 有帮助, 用udm未果,所以才用uds。udm个数已经打开,梯度保存也打开了。 2012-05-17 09:37:01
4楼正解。
简单的一个,干嘛要用uds?直接用udm就可以了。
5楼2012-05-16 22:33:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

b07988

金虫 (正式写手)

星期六可以休息

送鲜花一朵
引用回帖:
4楼: Originally posted by 47425804 at 2012-05-16 22:15:38:
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比用 ...

按你的思路dT会恒等于0.
应该是C_UDMI(c,t,0)=C_T_M1(c,t);吧。但是M1的话这样的话会直接报错。
所以我用了C_UDSI_M1(c,t,0)
先做最重要的事情,再做次重要的事情
6楼2012-05-17 09:24:09
已阅   回复此楼   关注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的回帖

hhucy

实习版主 (著名写手)

DOE锅炉工

uds可以使用的吧
人生那么多不确定,你怕什么
8楼2012-05-17 14:53:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

b07988

金虫 (正式写手)

星期六可以休息

引用回帖:
8楼: Originally posted by hhucy at 2012-05-17 14:53:49:
uds可以使用的吧

嗯,uds可以的。
先做最重要的事情,再做次重要的事情
9楼2012-05-17 16:33:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

b07988

金虫 (正式写手)

星期六可以休息

送鲜花一朵
引用回帖:
7楼: Originally posted by 47425804 at 2012-05-17 11:48:27:
恩,你说的没错,这样确实会恒等于0。
但是C_T_M1(c,t)是没用的。
等于零是因为计算源项的时候温度还没有更新,所以两个温度值是相等的。
这个时候其实是需要两个UDM,分别保存上一步和再上一步的温度。
而且 ...

这个方法挺好的,只要一开始的用DEFINE_INIT,对
C_UDMI(c,t,0)=0.;
C_UDMI(c,t,1)=0.;
就可以了!
我试了下,没有什么问题。谢谢了!

我也用uds试了下,其实也没有什么问题。只要打开uds就可以了。
先做最重要的事情,再做次重要的事情
10楼2012-05-17 16:45:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 b07988 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 296求调剂 +3 汪!?! 2026-04-05 4/200 2026-04-05 20:13 by 啵啵啵0119
[考研] 一志愿211生物学280分 求调剂 +4 李rien 2026-04-05 4/200 2026-04-05 18:01 by kk112233
[考研] 275求调剂 +16 waltzh 2026-04-01 16/800 2026-04-05 17:14 by Hdyxbekcb
[考研] 一志愿南航,数一英一学硕317求调剂!! +5 Acaciad 2026-04-04 5/250 2026-04-05 12:31 by 搏击518
[考研] 电子信息调剂交叉学科有推荐吗 +6 jhtfeybgj 2026-04-01 9/450 2026-04-05 11:13 by 猪会飞
[考研] 323求调剂(计算机视觉和大模型项目经历) +3 chaoxiicy 2026-03-31 3/150 2026-04-05 10:33 by zhq0425
[考研] 266分,一志愿电气工程,本科材料,求材料专业调剂 +11 哇呼哼呼哼 2026-04-01 12/600 2026-04-04 23:17 by 永字号
[考研] 材料334求调剂 +15 Eecho# 2026-04-03 15/750 2026-04-04 23:05 by 无际的草原
[考研] 302求调剂一志愿华中师范大学 +8 小江小江江江 2026-04-02 8/400 2026-04-04 19:50 by 蓝云思雨
[考研] 一志愿南农090401,268,求调剂 +5 一木鸟然 2026-04-04 5/250 2026-04-04 17:07 by babysonlkd
[考研] 303求调剂 +10 DLkz1314. 2026-03-30 10/500 2026-04-03 18:03 by Jimmyandyou
[考研] 机械专硕297 +3 Afksy 2026-04-03 3/150 2026-04-03 14:24 by 1753564080
[考研] 315求调剂 +6 顺理成张 2026-04-03 8/400 2026-04-03 14:04 by 百灵童888
[考研] 321求调剂 +17 y-yh 2026-04-01 20/1000 2026-04-03 12:57 by y-yh
[考研] 一志愿华南师范大学-22408计算机-292分-求华南师范大学调剂 +4 爱读书的小鳄鱼 2026-04-02 4/200 2026-04-02 18:35 by 求调剂zz
[考研] 【求调剂】新能源材料本科,一志愿211,初试321 +6 求调剂学校, 2026-04-02 6/300 2026-04-02 09:41 by 晴空210210
[考研] 省双一流重点一本大学招收调剂 +4 wwwwffffff 2026-03-31 7/350 2026-04-01 15:23 by wwwwffffff
[考研] 土木304求调剂 +5 顶级擦擦 2026-03-31 5/250 2026-04-01 08:15 by fdcxdystjk¥
[考研] 一志愿大连理工大学,机械工程学硕,341 +3 西瓜田的守望者 2026-03-30 3/150 2026-03-31 11:08 by asdfzly
[考研] 吉大生物学326分求调剂 +3 sunnyupup 2026-03-31 3/150 2026-03-31 09:28 by longlotian
信息提示
请填处理意见