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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » FLUENT & CFX » 能量源项中有温度对时间的微分，需要用uds吗？
171/2返回列表12下一页
查看: 5128  |  回复: 16
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

b07988

金虫 (正式写手)

星期六可以休息

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

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

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
先做最重要的事情，再做次重要的事情
1楼 2012-05-16 10:45:49
已阅   回复此楼   关注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，能用就是好东西，呵呵~~
赞一下(3人) 回复此楼
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");

但是没有效果，请帮我看看。
赞一下(1人) 回复此楼
先做最重要的事情，再做次重要的事情
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朵）

b07988
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朵）

b07988
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 的主题更新
171/2返回列表12下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 药学383 求调剂 +3 药学chy 2026-03-15 4/200 2026-03-16 20:51 by 元子^0^
[考研] 326求调剂 +4 诺贝尔化学奖觊� 2026-03-15 7/350 2026-03-16 17:11 by 诺贝尔化学奖觊�
[考研] 304求调剂 +3 曼殊2266 2026-03-14 3/150 2026-03-16 16:39 by houyaoxu
[考研] 283求调剂 +10 小楼。 2026-03-12 14/700 2026-03-16 16:08 by 13811244083
[考研] 309求调剂 +5 花与叶@ 2026-03-10 5/250 2026-03-16 14:13 by 哦哦123
[考博] 欢迎申博同学联系 +3 天道酬勤2026686 2026-03-10 7/350 2026-03-15 19:03 by 天道酬勤2026686
[考研] 288求调剂 +4 奇点0314 2026-03-14 4/200 2026-03-14 23:04 by JourneyLucky
[考研] 255求调剂 +3 李嘉慧， 2026-03-12 4/200 2026-03-14 16:58 by 有只狸奴
[考研] 学硕285求调剂 +13 Wisjxn 2026-03-12 46/2300 2026-03-14 10:33 by JourneyLucky
[考研] 0856材料与化工309分求调剂 +6 ZyZy…… 2026-03-10 6/300 2026-03-14 00:38 by JourneyLucky
[考研] 279求调剂 +3 抓着星星的女孩 2026-03-10 3/150 2026-03-13 23:47 by userper
[考研] 材料工程,326分，求调剂 +6 KRSLSR 2026-03-10 6/300 2026-03-13 23:47 by JourneyLucky
[考研] 【考研调剂求收留】 +3 Ceciilia 2026-03-11 3/150 2026-03-13 20:18 by JourneyLucky
[考研] 302求调剂 +6 负心者当诛 2026-03-11 6/300 2026-03-13 16:11 by JourneyLucky
[考研] 290求调剂 +7 ADT 2026-03-12 7/350 2026-03-13 15:17 by JourneyLucky
[考研] 081200-11408-276学硕求调剂 +3 崔wj 2026-03-12 4/200 2026-03-12 19:33 by 求调剂zz
[考研] 341求调剂 +4 捣蛋猪猪 2026-03-11 4/200 2026-03-12 14:47 by ruiyingmiao
[考研] 420求调剂 +4 莫向外求11 2026-03-10 6/300 2026-03-12 14:41 by ruiyingmiao
[考研] 化工0817调剂 +8 灿若星晨 2026-03-10 8/400 2026-03-10 22:44 by 星空星月
[考研] 0703化学调剂 +3 三dd. 2026-03-10 3/150 2026-03-10 15:45 by peike
信息提示
关闭
请填处理意见
关闭