| 查看: 1322 | 回复: 6 | ||||
uu1991金虫 (正式写手)
|
[求助]
fluent多空介质区的动量源相 已有3人参与
|
|
动量源相S=-粘性*速度/多空介质渗透率 多空介质目前只想到通过温度来确定，研究区域采用cell-loop循环来筛选，具体的源相编写如下 /*多空介质之间的动量源相*/ x方向 #include "udf.h" DEFINE_SOURCE(s_x,cell,thread,dS,eqn) { real k0,n,d0,Ts,Tl,k,fs; real T,u; real source; n=0.05; /*粘性*/ d0=1e-6; /*颗粒直径*/ Ts=415; /*固相线温度*/ Tl=425; /*液相线温度*/ k0=d0*d0/180; begin_c_loop(cell,thread) { T=C_T(cell,thread); if(Ts<T<Tl) { fs=(Tl-T)/(Tl-Ts); k=k0*(1-fs)*(1-fs)*(1-fs)/((fs*fs)+0.001); u=C_U(cell,thread); source=-n*u/k; dS[eqn]=0; } else source=dS[eqn]=0; return source; } end_c_loop(cell,thread) } 我研究的激光焊接，采用的是凝固&&融化模型，所以焊接材料应该都是多空介质区fluid,同时y、z方向和类似 在加载运行很短的时间内，软件提示温度发散，错误显示如下： Error: Divergence detected in AMG solver: temperature Error: Divergence detected in AMG solver: temperature Error Object: #f 搞了好长时间了，还是没有结果。 希望大家给点参考意见~~谢谢大家 |
回复此楼» 猜你喜欢
拟解决的关键科学问题还要不要写
已经有11人回复
-
救命帖
已经有5人回复
-
限项规定
已经有5人回复
-
为什么nbs上溴 没有产物点出现呢
已经有9人回复
-
招博士
已经有3人回复
-
存款400万可以在学校里躺平吗
已经有35人回复
-
最失望的一年
已经有18人回复
-
求推荐博导
已经有4人回复
-
求推荐英文EI期刊
已经有5人回复
-
疑惑？
已经有5人回复
» 本主题相关商家推荐: (我也要在这里推广)
匿名
用户注销 (小有名气)
- 应助: 31 (小学生)
- 金币: 1248.2
- 散金: 68
- 红花: 7
- 帖子: 220
- 在线: 78.3小时
- 虫号: 0
- 注册: 2015-01-18
- 专业: 工程热物理与能源利用
2楼2016-10-24 19:05:48
dai26520
专家顾问 (著名写手)
-
专家经验: +219 - 仿真EPI: 1
- 应助: 224 (大学生)
- 金币: 4097
- 散金: 4143
- 红花: 107
- 帖子: 1850
- 在线: 529.6小时
- 虫号: 1619980
- 注册: 2012-02-16
- 专业: 高分子合成化学
- 管辖: 仿真模拟
3楼2016-10-24 19:26:07
uu1991
金虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 1632.4
- 散金: 62
- 红花: 1
- 帖子: 312
- 在线: 107.8小时
- 虫号: 3347621
- 注册: 2014-07-31
- 性别: GG
- 专业: 激光
4楼2016-10-24 19:38:16
dai26520
专家顾问 (著名写手)
-
专家经验: +219 - 仿真EPI: 1
- 应助: 224 (大学生)
- 金币: 4097
- 散金: 4143
- 红花: 107
- 帖子: 1850
- 在线: 529.6小时
- 虫号: 1619980
- 注册: 2012-02-16
- 专业: 高分子合成化学
- 管辖: 仿真模拟
5楼2016-10-24 20:23:10
dai26520
专家顾问 (著名写手)
-
专家经验: +219 - 仿真EPI: 1
- 应助: 224 (大学生)
- 金币: 4097
- 散金: 4143
- 红花: 107
- 帖子: 1850
- 在线: 529.6小时
- 虫号: 1619980
- 注册: 2012-02-16
- 专业: 高分子合成化学
- 管辖: 仿真模拟
【答案】应助回帖
★
uu1991: 金币+1, ★有帮助, 这个程序我之前就看过，但是没有被启发到 2016-10-25 13:49:18
uu1991: 金币+1, ★有帮助, 这个程序我之前就看过，但是没有被启发到 2016-10-25 13:49:18
|
/********************************************************************* UDF for specifying a temperature-dependent viscosity property **********************************************************************/ #include "udf.h" DEFINE_PROPERTY(cell_viscosity, cell, thread) { real mu_lam; real temp = C_T(cell, thread); if (temp > 288.) mu_lam = 5.5e-3; else if (temp > 286.) mu_lam = 143.2135 - 0.49725 * temp; else mu_lam = 1.; return mu_lam; } |
6楼2016-10-24 20:31:32
zhouxman
金虫 (小有名气)
- 应助: 51 (初中生)
- 金币: 1363.4
- 散金: 32
- 红花: 5
- 帖子: 299
- 在线: 115.5小时
- 虫号: 1419015
- 注册: 2011-09-27
- 专业: 零件成形制造
7楼2016-10-25 20:43:36