| 查看: 534 | 回复: 0 | ||
[求助]
Fluent UDF 求助
|
|
目的：将计算的每一个迭代步的流量导出，并利用欧拉插值得到压力值，并赋回给Fluent作为边界条件，进行下一个迭代步的计算，一直到计算收敛，该时间步结束。 问题：仅仅在第一个时间步的时候发现压力利用欧拉迭代产生了变化，在之后的计算中压力始终等于我赋的初值。求大牛解答这是为什么。 代码如下： DEFINE_PROFILE(DFA_outlet,thread,position) { real NV_VEC(A),NV_VEC(FaceVel); double Q; face_t f; int i; real t=RP_Get_Real("flow-time"  ;real tt; // define the time in one cycle real w=0.8; // period of a cycle real p1[100]; // define the pressure after R1 double R1,R2,C; // define electric modulus int TTS; // total time steps int h,TS; // time interval and time step R1=9.717571*pow(10,8),R2=3.849843*pow(10,9),C=3.712461*pow(10,-10); // set the values of electric modulus TTS=100; // total time steps begin_f_loop(f,thread) { NV_D(FaceVel,=,F_U(f,thread),F_V(f,thread),F_W(f,thread)); F_AREA(A,f,thread); Q=NV_DOT(FaceVel,A); //compute the pressure using the explicit Euler method h=w/TTS; tt=t-floor(t/w)*w; TS=floor(tt/h); p1[0]=12102.65-Q*R1; if(TS=0) F_PROFILE(f,thread,position)=p1[0]+Q*R1; else for(i=1;i<=TS;i++) { p1[TS]=p1[TS-1]+h*(Q-(p1[TS-1]-1333.22)/R2)/C; } F_PROFILE(f,thread,position)=p1[TS]+Q*R1; } end_f_loop(f,thread) } @月只蓝 @beefly 发自小木虫Android客户端 |
回复此楼» 猜你喜欢
欢迎加入课题组
已经有0人回复
-
散金币，求好运，祝面上顺利！
已经有80人回复
-
物理学I论文润色/翻译怎么收费?
已经有209人回复
-
华东师范大学芯片设计徐珑真实水平咋样
已经有1人回复
-
调剂
已经有0人回复
-
【新加坡】纳米电子器件项目组有“联合培养博士生”名额
已经有0人回复
-
中德博士后2026项目5月30日截止- 国内博士后入站、德国两年全职
已经有0人回复
-
面向全球申请人的中德博士后2026项目5月30日截止- 国内博后入站德国两年全职
已经有0人回复
-
【新加坡】纳米电子器件项目组有“联合培养博士生”名额
已经有2人回复
-
重庆交大材料学院26年硕士招生复试实施细则已出。欢迎加入光子学微结构与器件课题组。
已经有0人回复
-
北工大宋晓艳教授团队招收2026年9月入学博士【学博和专博名额均有】！
已经有1人回复
