| 查看: 544 | 回复: 0 | ||
[求助]
用fluent做焊接模拟。程序编译通过了，计算不收敛，求大神看一下啊，谢谢
|
|
#include "udf.h" #define Mu 4*3.1415927e-7 #define T0 293 #define hc 80 #define dsity 1810 #define g 9.81 #define bata 1e-4 DEFINE_PROFILE(heat_flux,t,i) { real x[ND_ND]; face_t f; real r,m,n,k,qmax,a1,b1,c1,a2,b2,c2,s,e; qmax=8e6; a1=7.648e6; b1=0.5133; c1=0.5986; a2=7.647e6; b2=-0.5133; c2=0.5985; s=5.67e-8; e=0.06; begin_f_loop(f,t) { real T=F_T(f,t); F_CENTROID(x,f,t); m=x[0]; n=x[1]; k=RP_Get_Real("flow-time"  ;r=sqrt(pow(m-2*k,2)+pow(n,2)); if(r<=3) F_PROFILE(f,t,i)=a1*exp(-pow((r-b1)/c1,2))+a2*exp(-pow((r-b2)/c2,2))-hc*(T-T0)-s*e*(pow(T,4)-pow(T0,4)); else F_PROFILE(f,t,i)=-hc*(T-T0)-s*e*(pow(T,4)-pow(T0,4)); } end_f_loop(f,t) } DEFINE_SOURCE(x_velocity_source,cell,thread,dS,eqn) { real source; real x[ND_ND]; real r,m,y,z,k,Jr; begin_f_loop(cell,thread) { C_CENTROID(x,cell,thread); m=x[0]; y=x[1]; z=x[2]; k=RP_Get_Real("flow-time" ;r=sqrt(pow(m-2*k,2)+pow(y,2)); Jr=7.473e6*exp(-pow((r-1.211)/1.451,2))+7.473e6*exp(-pow((r+1.211)/1.451,2)); source=-Mu*Jr*(1-Jr/(7.473e6))*(1-z/5)*m/(12*3.1415926*pow(r,2)); dS[eqn]=0; return source; } end_f_loop(cell,cell_thread) } DEFINE_SOURCE(y_velocity_source,cell,thread,dS,eqn) { real source; real x[ND_ND]; real r,m,y,z,k,Jr; begin_f_loop(cell,thread) { C_CENTROID(x,cell,thread); m=x[0]; y=x[1]; z=x[2]; k=RP_Get_Real("flow-time" ;r=sqrt(pow(m-2*k,2)+pow(y,2)); Jr=7.473e6*exp(-pow((r-1.211)/1.451,2))+7.473e6*exp(-pow((r+1.211)/1.451,2)); source=-Mu*Jr*(1-Jr/(7.473e6))*(1-z/5)*y/(12*3.1415926*pow(r,2)); dS[eqn]=0; return source; } end_f_loop(cell,cell_thread) } DEFINE_SOURCE(z_velocity_source,cell,thread,dS,eqn) { real source; real T=C_T(cell,thread); real x[ND_ND]; real r,m,y,z,k,Jr; begin_f_loop(cell,thread) { C_CENTROID(x,cell,thread); m=x[0]; y=x[1]; z=x[2]; k=RP_Get_Real("flow-time" ;r=sqrt(pow(m-2*k,2)+pow(y,2)); Jr=7.473e6*exp(-pow((r-1.211)/1.451,2))+7.473e6*exp(-pow((r+1.211)/1.451,2)); source=Mu*7473000*(1-Jr/7473000)*(1-z/5)/(12*3.1415926*pow(r,2))-dsity*g*bata*(T-800); dS[eqn]=0; return source; } end_f_loop(cell,cell_thread) } 编译完成后，三个方向的动量源项“不是所有的控件路径都返回值” 这是什么情况啊？而且导入fluent中，计算也不收敛 |
回复此楼» 猜你喜欢
欢迎发来filecode的Mz6后的代码验证其规律
已经有118人回复
-
今天放榜吗？
已经有15人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有3人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有4人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有5人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有5人回复
-
filecode，4个jtjc了
已经有11人回复
-
filecode=后面第一个是大写字母
已经有12人回复
-
为什么余额宝的年化利率越来越低？主要原因有哪些？
已经有4人回复
-
今天基金会出结果吗？20260819
已经有8人回复
找到一些相关的精华帖子，希望有用哦~
关于Fluent运算化学反应开始后温度场异常的问题
已经有3人回复
- 求助大神，comsol仿真一直出现不收敛
已经有5人回复
- fluent无法收敛怎么办？怎么判断已经收敛？
已经有12人回复
- fluent非稳态计算两秒后出现不收敛，求指点，谢谢哦
已经有19人回复
- fluent非稳态计算每步都必须收敛吗？偶尔的一步不收敛对最终的计算结果影响大吗？
已经有8人回复
- fluent输出温度场最大值最小值和平均值
已经有4人回复
- fluent14编译UDF时总是出错，求高手啊！！！难道没人遇到过吗？？？
已经有11人回复
- 用fluent模拟时，选择非稳态求解器，收敛以后还继续算会不会对结果有影响？
已经有7人回复
- Fluent连续性方程不收敛
已经有13人回复
- 新人求助，fluent网格优化问题！大神请入
已经有12人回复
- FLUENT计算流场求稳态计算收敛办法
已经有3人回复
- fluent 计算收敛后，看不到结果
已经有4人回复
- fluent中continuity曲线不收敛
已经有8人回复
- fluent 模拟气液两相射流不收敛
已经有3人回复
- 用Fluent中的mixture模型做气固两相流仿真时，粒径的设置及收敛问题
已经有17人回复
- 翼型网格在fluent中收敛性差或不收敛
已经有4人回复
- fluent过来人经验谈之continuity不收敛的问题
已经有36人回复
- fluent新手的福音---fluent不收敛的解决办法
已经有193人回复
- 怎样判定fluent计算是否收敛
已经有4人回复
- fluent温度场分析结果不收敛
已经有6人回复
- FLUENT不收敛通常的解决方式.pdf
已经有115人回复
- 【求助】Fluent 收敛问题
已经有10人回复
- 【求助】fluent中非稳态计算的收敛【已解决】
已经有5人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人
