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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » FLUENT & CFX » fluent多空介质区的动量源相
71/1返回列表
查看: 1541  |  回复: 6
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

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
搞了好长时间了，还是没有结果。 希望大家给点参考意见~~谢谢大家
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)
1楼 2016-10-24 18:53:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

匿名

用户注销 (小有名气)
★ ★
感谢参与，应助指数 +1
uu1991: 金币+2, ★有帮助 2016-10-24 19:34:04
本帖仅楼主可见
赞一下
2楼2016-10-24 19:05:48
已阅   申请仿真EPI   回复此楼   编辑   查看我的主页

dai26520

专家顾问 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与，应助指数 +1
uu1991: 金币+5, ★有帮助 2016-10-24 19:38:29
Define-Source不需要遍历循环，第一个DS≠0，另外是否需要提取多孔介质温度？
赞一下(1人) 回复此楼
个性只是一个代号，做好自己就行
3楼2016-10-24 19:26:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

uu1991

金虫 (正式写手)
引用回帖:
3楼: Originally posted by dai26520 at 2016-10-24 19:26:07
Define-Source不需要遍历循环，第一个DS≠0，另外是否需要提取多孔介质温度？

我想在固液相温度间的区域加载这个动量源相，还有是目前我还没想到其他的想法确定这个区域，所以采取提取温度的这个方法
赞一下 回复此楼
4楼2016-10-24 19:38:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dai26520

专家顾问 (著名写手)
引用回帖:
4楼: Originally posted by uu1991 at 2016-10-24 19:38:16
我想在固液相温度间的区域加载这个动量源相，还有是目前我还没想到其他的想法确定这个区域，所以采取提取温度的这个方法...

我建议你看看官网UDF帮助里面最后Example里面凝固那个udf例子，可能有帮助
赞一下 回复此楼
个性只是一个代号，做好自己就行
5楼2016-10-24 20:23:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dai26520

专家顾问 (著名写手)

【答案】应助回帖

★
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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhouxman

金虫 (小有名气)

【答案】应助回帖

感谢参与，应助指数 +1
fluent 自己就有固液共存的多孔介质区的 选项，没有必要自己写把
赞一下 回复此楼
7楼2016-10-25 20:43:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 uu1991 的主题更新
71/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 重要来源：本周末出结果 +3 瞬息宇宙 2026-08-12 3/150 2026-08-13 06:40 by 玲珑133
[基金申请] FileCode能看出啥？ +8 要乐观耀哥 2026-08-10 23/1150 2026-08-13 00:35 by 要乐观耀哥
[基金申请] 长年满屏的广告，版主太不责任了。基金也等的急 +4 gltch 2026-08-12 4/200 2026-08-12 16:08 by eulota
[基金申请] 奇怪，两个人的filecode固定段从头到尾一模一样 +7 布布和一二 2026-08-10 10/500 2026-08-12 09:12 by 布布和一二
[基金申请] 小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊 +8 Tide man 2026-08-10 9/450 2026-08-11 20:39 by beefly
[硕博家园] 读博的好处 +3 lnee 2026-08-11 3/150 2026-08-11 18:10 by 希望我好好的
[基金申请] 什么时候出结果，有咨询渠道？？？ +3 Tide man 2026-08-11 3/150 2026-08-11 17:54 by kudofaye
[基金申请] 分享一下我之前已中青C的计划书的filecode +3 布布和一二 2026-08-11 4/200 2026-08-11 15:20 by aasahr
[基金申请] 基金中了 +15 laoda193707 2026-08-06 15/750 2026-08-11 00:11 by jiafei2190
[基金申请] 确定了，国自然21号放榜 +6 布布和一二 2026-08-10 7/350 2026-08-10 19:15 by 2000zf36392
[基金申请] 关于代码变化问题，想知道的进来 +17 且听虎啸 2026-08-07 24/1200 2026-08-10 18:35 by zhangduo2008
[基金申请] 国自然结果 +4 Vierhys 2026-08-10 8/400 2026-08-10 15:06 by Vierhys
[基金申请] 据悉今年马上要出结果了 +7 瞬息宇宙 2026-08-10 8/400 2026-08-10 12:42 by Vivilian
[基金申请] 这样的filecode谁见过 +11 布布和一二 2026-08-08 22/1100 2026-08-10 11:10 by wmfsnow
[基金申请] 面上项目filecode邪修 +5 西山十月 2026-08-09 7/350 2026-08-10 07:32 by 仁砚薪传
[基金申请] fileCode有新解读？ +10 Tide man 2026-08-08 18/900 2026-08-09 12:55 by 仁砚薪传
[基金申请] 关于filecode，很负责任的告诉大家 +6 爱看书的可乐 2026-08-08 7/350 2026-08-08 22:13 by a_niu
[基金申请] 国基金的申报应该改成非等额制，评价高的钱多评价低的钱少，但是增加资助率 +7 a089 2026-08-07 7/350 2026-08-08 18:05 by gltch
[基金申请] 固定端突然变了，今天 +6 archvillain 2026-08-06 10/500 2026-08-07 16:03 by 医学老男孩
[基金申请] filecode +8 布布和一二 2026-08-06 11/550 2026-08-06 20:41 by tangpu318
信息提示
关闭
请填处理意见
关闭