24小时热门版块排行榜    

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

b07988

金虫 (正式写手)

星期六可以休息

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

能量源项中有dT/dt,就是温度对时间的求导,请问是用
1.c_udsi_m1存储上一步时间步长的温度
2.采用define_uds_unsteady?
希望有好心人帮忙。
不胜感激!
回复此楼
先做最重要的事情,再做次重要的事情
已阅   回复此楼   关注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的回帖
查看全部 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");

但是没有效果,请帮我看看。
先做最重要的事情,再做次重要的事情
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朵)

4楼2012-05-16 22:15:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
信息提示
请填处理意见