| 查看: 2492 | 回复: 13 | ||
[求助]
温度区域能量源项的添加
|
|
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; } |
回复此楼» 猜你喜欢
（调剂）一志愿报考哈尔滨工业大学0857资源与环境专业378分考生
已经有6人回复
-
262求调剂
已经有7人回复
-
材料专硕322分
已经有10人回复
-
材料调剂
已经有8人回复
-
085600，320分求调剂
已经有8人回复
-
一志愿 江南大学 085602 化工专硕 338分求调剂
已经有11人回复
-
求调剂
已经有6人回复
-
086000调剂
已经有3人回复
-
求调剂
已经有4人回复
-
0703调剂，一志愿天津大学319分
已经有4人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 能量方程用温度计算还是焓计算？
已经有6人回复
- 请教耦合面传热问题
已经有6人回复
- UDS边界条件，求知道的说一下，困惑我Ｎ天了...
已经有5人回复
- 能量源项加载问题
已经有19人回复
- 高斯09计算分子能量时默认温度是多少？
已经有4人回复
- 【求助】FLUENT UDF 气液两相流 汽蚀 汽化潜热(问题详见内容)
已经有21人回复
- UDF 热源随温度变化的处理
已经有4人回复
- UDF添加动能源项，发现温度没什么变化
已经有20人回复
- 能量源项中有温度对时间的微分，需要用uds吗？
已经有16人回复
- fluent能量源项中温度T对时间的t的偏导数，即dT/dt如何用udf编写能量源项
已经有6人回复
- 谁能告诉我一下用分子动力学怎么计算能量和温度，谢谢了啊
已经有7人回复
- 能量方程源项udf
已经有4人回复
- 求助气固两相流的边界条件问题
已经有8人回复
- 水沸腾UDF总是报错
已经有10人回复
- 关于相态转变的模拟
已经有9人回复
- 【求助】vasp 算出的0k下的能量如何扩展到其他温度压力
已经有12人回复
- 【讨论】很奇怪的能量，温度曲线。
已经有7人回复
- 【求助】求助或讨论一下关于能量传递源项中各个焓值的含义
已经有3人回复
- 【其他】vasp结合phonon计算材料不同温度下能量的问题
已经有9人回复
- 【求助】OUT 文件显示怎么不同温度和压力下的能量怎么是相同的啊
已经有4人回复
- 【讨论】请教一个气-液-固三相流模拟可能性的问题
已经有6人回复
- 【求助】冷凝的UDF程序
已经有15人回复
shuai_1314
捐助贵宾 (小有名气)
- 应助: 33 (小学生)
- 金币: 1861.7
- 红花: 6
- 帖子: 167
- 在线: 163.1小时
- 虫号: 2415320
- 注册: 2013-04-12
- 性别: GG
- 专业: 内流流体力学
2楼2013-05-14 11:09:46
3楼2013-05-14 13:24:51
4楼2013-05-14 20:44:11
shuai_1314
捐助贵宾 (小有名气)
- 应助: 33 (小学生)
- 金币: 1861.7
- 红花: 6
- 帖子: 167
- 在线: 163.1小时
- 虫号: 2415320
- 注册: 2013-04-12
- 性别: GG
- 专业: 内流流体力学
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★
xiegangmai: 金币+2, 谢谢参与！ 2013-05-15 11:25:05
seazsh: 金币+5, ★有帮助 2013-05-15 16:52:34
xiegangmai: 金币+2, 谢谢参与！ 2013-05-15 11:25:05
seazsh: 金币+5, ★有帮助 2013-05-15 16:52:34
|
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
shuai_1314
捐助贵宾 (小有名气)
- 应助: 33 (小学生)
- 金币: 1861.7
- 红花: 6
- 帖子: 167
- 在线: 163.1小时
- 虫号: 2415320
- 注册: 2013-04-12
- 性别: GG
- 专业: 内流流体力学
6楼2013-05-14 23:02:44
7楼2013-05-15 08:55:30
8楼2013-05-15 08:57:10
9楼2013-05-15 14:02:03
10楼2013-05-15 14:05:12
