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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » FLUENT & CFX » 能量源项中有温度对时间的微分，需要用uds吗？
北京石油化工学院2026年研究生招生接收调剂公告
51/1返回列表
查看: 5186  |  回复: 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的回帖

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的回帖
查看全部 17 个回答

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的回帖

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的回帖

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的回帖
查看全部 17 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 环境专硕调剂 +6 会说话的肘子 2026-04-06 6/300 2026-04-07 02:00 by Sealedwind
[考研] 070300化学学硕311分求调剂 +12 梁富贵险中求 2026-04-04 12/600 2026-04-06 22:54 by chenzhimin
[考研] 求助 +3 卡卡东88 2026-04-06 4/200 2026-04-06 15:28 by going home
[考研] 专硕0854初试考材科基，求调剂 +6 3220548044 2026-04-06 9/450 2026-04-06 10:26 by barlinike
[考研] 322求调剂 +3 嗯哼哼恒 2026-04-05 3/150 2026-04-05 19:52 by nepu_uu
[考研] 277求调剂 +5 考研调剂lxh 2026-04-05 5/250 2026-04-05 19:03 by chy09050039
[考研] 329求调剂 +17 miaodesi 2026-04-02 20/1000 2026-04-05 18:33 by 蓝云思雨
[考研] 271分求调剂学校 +12 zph158488！ 2026-04-02 13/650 2026-04-05 10:13 by lqwchd
[考研] 353求调剂 +10 MayUxw1 2026-04-03 10/500 2026-04-05 09:23 by 无际的草原
[考研] 085602调剂 初试总分335 +12 19123253302 2026-04-04 12/600 2026-04-05 08:08 by 544594351
[考研] 283分求调剂 +7 小聂爱学习 2026-04-03 7/350 2026-04-04 21:51 by hemengdong
[考研] 316求调剂 +9 墨辰_Orion926 2026-04-04 9/450 2026-04-04 21:35 by lbsjt
[考研] 338求调剂 +7 晟功? 2026-04-03 7/350 2026-04-04 20:37 by 蓝云思雨
[考研] 一志愿南农090401，268，求调剂 +5 一木鸟然 2026-04-04 5/250 2026-04-04 17:07 by babysonlkd
[考研] 286求调剂 +8 lim0922 2026-04-02 8/400 2026-04-03 20:19 by rzh123456
[考研] 314求调剂 +11 1xiaojun23 2026-03-31 12/600 2026-04-02 12:31 by 1xiaojun23
[考研] 材料专硕322分 +11 哈哈哈吼吼吼哈 2026-04-01 11/550 2026-04-02 10:52 by lnilvy
[考研] 0710生物学求调剂 +9 manman511 2026-04-01 9/450 2026-04-02 10:00 by zxl830724
[考研] 材料调剂 +11 一样YWY 2026-03-31 11/550 2026-04-01 11:35 by wangjy2002
[考研] 269求调剂 +4 我想读研11 2026-03-31 4/200 2026-03-31 10:04 by cal0306
信息提示
关闭
请填处理意见
关闭