fluent中MHD模块用UDF添加外部磁场
用udf的形式给MHD模块中外部磁场赋值,赋值的效果想跟自带的MHD模块一样,结果计算的时候发现,外部磁场赋值成功了,计算的时候却发散了我计算的问题是磁流体泵;变化的磁场驱动流体运动。
我的udf:
#include "udf.h"
DEFINE_INIT(my_init_BO_y,domain)
{
cell_t c;
Thread *t;
real xc;
real temp1,temp2;
temp1=0.02;
temp2=0.2;
/* loop over all cell threads in the domain */
thread_loop_c(t,domain)
{
/* loop over all cells */
begin_c_loop_all(c,t)
{
C_CENTROID(xc,c,t);
if(xc<1.0) C_UDMI(c,t,2)=cos(2*3.1415*60.0/360.0);
else if(xc<2.0) C_UDMI(c,t,2)=cos(2*3.1415*120.0/360.0);
else if(xc<3.0) C_UDMI(c,t,2)=cos(2*3.1415*180.0/360.0);
else if(xc<4.0) C_UDMI(c,t,2)=cos(2*3.1415*240.0/360.0);
else if(xc<5.0) C_UDMI(c,t,2)=cos(2*3.1415*300.0/360.0);
elseC_UDMI(c,t,2)=cos(2*3.1415*360.0/360.0);
C_UDMI(c,t,4)=sqrt(C_UDMI(c,t,2)*C_UDMI(c,t,2));
}
end_c_loop_all(c,t)
}
}
DEFINE_ADJUST(adjust_b0_y,d)
{
cell_t c;
Thread *t;
real xc;
real current_time;
current_time=RP_Get_Real("flow-time");
/* current_time=CURRENT_TIME;*/
thread_loop_c(t,d)
{
begin_c_loop_all(c,t)
{
C_CENTROID(xc,c,t);
if(xc<1.0) C_UDMI(c,t,2)=cos(2*3.1415*60.0/360.0+2.0*3.1415*current_time);
else if(xc<2.0) C_UDMI(c,t,2)=cos(2*3.1415*120.0/360.0+2.0*3.1415*current_time);
else if(xc<3.0) C_UDMI(c,t,2)=cos(2*3.1415*180.0/360.0+2.0*3.1415*current_time);
else if(xc<4.0) C_UDMI(c,t,2)=cos(2*3.1415*240.0/360.0+2.0*3.1415*current_time);
else if(xc<5.0) C_UDMI(c,t,2)=cos(2*3.1415*300.0/360.0+2.0*3.1415*current_time);
elseC_UDMI(c,t,2)=cos(2*3.1415*360.0/360.0+2.0*3.1415*current_time);
C_UDMI(c,t,4)=sqrt(C_UDMI(c,t,2)*C_UDMI(c,t,2));
}
end_c_loop_all(c,t)
}
}
现在主要的问题是:计算几步后就会发散,发现感应磁场不正确
这个是附图,供大家参考
请问楼主解决了吗?至于一定要先把电磁场计算出来,然后导入mhd吗?
页:
[1]