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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » FLUENT & CFX » 温度区域能量源项的添加
141/2返回列表12下一页
查看: 2735  |  回复: 13
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

seazsh

金虫 (小有名气)

[求助] 温度区域能量源项的添加

1想在某温度或者某温度区间施加能量源项， 按下面的编写好像加不上去，

2还有如果要在某温区添加动量源项怎么处理？

DEFINE_SOURCE(energy_source,c,t,dS,eqn)
{
   real source;
   real r;
   real x[ND_ND];
   real xi,yi;      
   real cs,temp;
cs=3/pow(R0,2);  
temp=C_T(c,t);   
C_CENTROID(x,c,t);
        xi=x[0];
        yi=x[1];
          r=sqrt(xi*xi);              
  if(temp<=235&&temp>=365&&yi>=0)
      {  
            source=3*cs*P/(3.14159*0.007*(1-1/exp(3)))*exp(-3*r*r*cs/log10(H/(yi+0.0001)));
           dS[eqn]=0;
      }
   else
  source=dS[eqn]=0;
return source;
}
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2013-05-13 23:04:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

shuai_1314

捐助贵宾 (小有名气)

【答案】应助回帖

★ ★
感谢参与，应助指数 +1
xiegangmai: 金币+2, 谢谢参与！ 2013-05-14 21:40:52
你的R0没有定义吧？还有你的return给的不对，要不就在每条条件判断语句后都return，要么在条件判断外部return

[ 发自手机版 http://muchong.com/3g ]
赞一下 回复此楼
率真随性的射手
2楼2013-05-14 11:09:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

seazsh

金虫 (小有名气)
引用回帖:
2楼: Originally posted by shuai_1314 at 2013-05-14 11:09:46
你的R0没有定义吧？还有你的return给的不对，要不就在每条条件判断语句后都return，要么在条件判断外部return

R0和H在前面定义了，可以编译，return这样好像也可以吧？
主要的问题是1想在某温度或者某温度区间施加能量源项， 按下面的编写好像加不上去，

2还有如果要在某温区添加动量源项怎么处理？
赞一下 回复此楼
3楼2013-05-14 13:24:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

b07988

金虫 (正式写手)

星期六可以休息

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
xiegangmai: 金币+2, 谢谢参与！ 2013-05-14 21:41:00
seazsh: 金币+5, ★有帮助 2013-05-15 16:52:49
if(temp<=235&&temp>=365&&yi>=0)????
错误其实挺明显的。。。
&&应该是||
还有
exp(-3*r*r*cs/log10(H/(yi+0.0001)));
自己再仔细看看原文献吧，公式不是这样的，亲。
其次，0.0001有必要吗？
自己好好想想吧。
赞一下 回复此楼
先做最重要的事情，再做次重要的事情
4楼2013-05-14 20:44:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

shuai_1314

捐助贵宾 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★
xiegangmai: 金币+2, 谢谢参与！ 2013-05-15 11:25:05
seazsh: 金币+5, ★有帮助 2013-05-15 16:52:34
引用回帖:
3楼: Originally posted by seazsh at 2013-05-14 13:24:51
R0和H在前面定义了，可以编译，return这样好像也可以吧？
主要的问题是1想在某温度或者某温度区间施加能量源项， 按下面的编写好像加不上去，

2还有如果要在某温区添加动量源项怎么处理？...

DEFINE_SOURCE(source_energe, cell, thread, dS, eqn)
{
        real energe;/*热源*/
        real C[ND_ND];
        real x;
        real y;
        real z;
   
       C_CENTROID(C,cell,thread);/*取控制体cell的坐标*/
            x=C[0];
            y=C[1];
            z=C[2];
            
            if(x >= 0.1 && x <= 0.3)
                    {
                            if(sqrt(y*y+z*z) <= 0.014)
                                    {
                                      energe = 500;
                                      dS[eqn] = 0;
                                      return energe;
                                    }
                            else
                                    {
                                      energe = 0;
                                      dS[eqn] = 0;
                                      return energe;
                                    }
                    }
            else
                    {
                            energe = 0;
                            dS[eqn] = 0;
                            return energe;
                    }
这是我以前编过的一个能量源项，判断过后马上接着return，我之前像你那样编过return，但是也加不上去！

动量源项分方向，注意单位！具体的值你要自己算一下嘛
赞一下 回复此楼
率真随性的射手
5楼2013-05-14 23:01:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

shuai_1314

捐助贵宾 (小有名气)
引用回帖:
4楼: Originally posted by b07988 at 2013-05-14 20:44:11
if(temp<=235&&temp>=365&&yi>=0)????
错误其实挺明显的。。。
&&应该是||
还有
exp(-3*r*r*cs/log10(H/(yi+0.0001)));
自己再仔细看看原文献吧，公式不是这样的，亲。
其次 ...

&&也是可以的
赞一下 回复此楼
率真随性的射手
6楼2013-05-14 23:02:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

b07988

金虫 (正式写手)

星期六可以休息
引用回帖:
6楼: Originally posted by shuai_1314 at 2013-05-14 23:02:44
&&也是可以的...

温度可以小于235，又大于365？？？
&&是并且的意思吧。。。。
赞一下 回复此楼
先做最重要的事情，再做次重要的事情
7楼2013-05-15 08:55:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

b07988

金虫 (正式写手)

星期六可以休息
★
xiegangmai: 金币+1, 谢谢参与！ 2013-05-15 11:25:14
引用回帖:
5楼: Originally posted by shuai_1314 at 2013-05-14 23:01:46
DEFINE_SOURCE(source_energe, cell, thread, dS, eqn)
{
        real energe;/*热源*/
        real C;
        real x;
        real y;
        real z;
   
       C_CENTROID(C,cell,thread);/*取 ...

最后再return，语法绝对没有问题。
事实上，就是没有问题。
赞一下 回复此楼
先做最重要的事情，再做次重要的事情
8楼2013-05-15 08:57:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

seazsh

金虫 (小有名气)
引用回帖:
4楼: Originally posted by b07988 at 2013-05-14 20:44:11
if(temp<=235&&temp>=365&&yi>=0)????
错误其实挺明显的。。。
&&应该是||
还有
exp(-3*r*r*cs/log10(H/(yi+0.0001)));
自己再仔细看看原文献吧，公式不是这样的，亲。
其次 ...

温度区间是2350到1365修改时没注意，能源公式随便给的，主要是用温度判据加载不上，用区域（坐标）可以加载。
赞一下 回复此楼
9楼2013-05-15 14:02:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

seazsh

金虫 (小有名气)
引用回帖:
5楼: Originally posted by shuai_1314 at 2013-05-14 23:01:46
DEFINE_SOURCE(source_energe, cell, thread, dS, eqn)
{
        real energe;/*热源*/
        real C;
        real x;
        real y;
        real z;
   
       C_CENTROID(C,cell,thread);/*取 ...

和你一样用区域（xyz坐标）判据可以加载，可是用温度判据就不行，不怎么回事
赞一下 回复此楼
10楼2013-05-15 14:05:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 seazsh 的主题更新
141/2返回列表12下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[教师之家] 某top大学教授说“能够在市场中兑现的能力才是真能力”无比同意！ +8 zju2000 2026-09-26 11/550 2026-09-30 08:04 by xli1984
[博后之家] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +4 IDs3scOF0tjC 2026-09-28 4/200 2026-09-29 19:27 by cqwQDCxMcL3I
[找工作] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +4 ZYXdhzDAy9ZX 2026-09-28 4/200 2026-09-29 16:54 by yCO1Ll7aHtsw
[考博] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +4 ZYXdhzDAy9ZX 2026-09-28 4/200 2026-09-29 16:47 by yCO1Ll7aHtsw
[博后之家] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 ZYXdhzDAy9ZX 2026-09-28 3/150 2026-09-29 16:30 by etmYJ6d2rquH
[硕博家园] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +3 ZYXdhzDAy9ZX 2026-09-28 3/150 2026-09-29 16:03 by etmYJ6d2rquH
[论文投稿] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +4 tqUTOxClQUMF 2026-09-28 4/200 2026-09-29 15:15 by etmYJ6d2rquH
[考博] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +6 GDBe8tDZqE8z 2026-09-28 6/300 2026-09-29 14:41 by etmYJ6d2rquH
[硕博家园] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +3 GDBe8tDZqE8z 2026-09-28 3/150 2026-09-29 14:30 by etmYJ6d2rquH
[考研] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +3 CfXuS1rDhLYN 2026-09-28 3/150 2026-09-29 14:24 by etmYJ6d2rquH
[找工作] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +6 IDs3scOF0tjC 2026-09-28 6/300 2026-09-29 14:13 by etmYJ6d2rquH
[博后之家] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +6 IDs3scOF0tjC 2026-09-28 7/350 2026-09-29 14:07 by etmYJ6d2rquH
[博后之家] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +3 IDs3scOF0tjC 2026-09-28 3/150 2026-09-29 13:56 by etmYJ6d2rquH
[公派出国] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +4 IDs3scOF0tjC 2026-09-28 5/250 2026-09-29 13:50 by etmYJ6d2rquH
[找工作] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +3 ZYXdhzDAy9ZX 2026-09-28 3/150 2026-09-29 06:28 by ZvyPK8n6Nfki
[硕博家园] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +3 CfXuS1rDhLYN 2026-09-28 4/200 2026-09-28 22:55 by ez6fGg9abYaj
[公派出国] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +4 IDs3scOF0tjC 2026-09-28 4/200 2026-09-28 22:55 by ez6fGg9abYaj
[有机交流] 同一个分子，一条来自文献，一条来自AI——不告诉你答案，你会选哪条？ +3 tianxiaoxian 2026-09-28 6/300 2026-09-28 22:44 by tianxiaoxian
[论文投稿] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +3 CfXuS1rDhLYN 2026-09-28 3/150 2026-09-28 19:08 by 9lS3ad5oOymn
[考研] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +3 IDs3scOF0tjC 2026-09-28 4/200 2026-09-28 18:51 by 9lS3ad5oOymn
信息提示
关闭
请填处理意见
关闭