24小时热门版块排行榜    

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

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;
}
回复此楼

» 猜你喜欢

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

已阅   回复此楼   关注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的回帖
查看全部 14 个回答

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的回帖
信息提示
请填处理意见