24小时热门版块排行榜    

Znn3bq.jpeg
汕头大学海洋科学接受调剂
查看: 1331  |  回复: 0

pjhust

新虫 (初入文坛)

[求助] fluent中UDF设定边界条件,导入出错

最近在做一个fluent模拟,需要设定非稳态边界条件,写了一个UDF程序如下,导入fluent里面时,出现:
Error: C:\Users\\305\313\346\257\Desktop\pexit2.c.c: line 18: parse error.
Error: C:\Users\\305\313\346\257\Desktop\pexit2.c.c: line 19: invalid lvalue in assignment.
的错误,但是鄙人实在找不出错误之处,所以请教各位大神(由于是新虫,总共只有5个金币,对不起各位了!)

UDF程序:


#include"udf.h"
#define r_1 20
#define r_2 15.4
#define a 0.77
#define w 418.879
#define H 18
#define t 3.2
#define m 1.2
#define p_s 587050
#define p_d 1980000
#define V_r pi*(pow(r_1,2)-pow(r_2,2))*H

DEFINE_PROFILE(unsteady_pressure,thread,position)
{
face_t f;
begin_f_loop(f,thread)
{
real t=RP_Get_Real("flow-time";
t=t-0.015*floor(t/0.015);
if(t<0.002083 && t>0)
   F_PROFILE(f,thread,position)=587050;
else if(t<=0.00845256 && t>=0.002083)
   real f_beta,h_beta,V_cbeta;
   f_beta=(1-pow(a,2))*418.879*t-0.5*pow((1-a),2)*sin(2.0*418.879*t)-pow(a,2)*asin((1.0/a-1)*sin(418.879*t))-a*(1-a)*sin(418.879*t)*pow((1-pow((1.0/a-1),2)*pow(sin(418.879*t),2)),0.5);
   h_beta=r_1*(1-(1-a)*cos(418.879*t)-sqrt(pow((1-a),2)*pow(cos(418.879*t),2)+2*a-1));
   V_cbeta=V_r-0.5*H*(pow(r_1,2)*f_beta+h_beta*t);

   F_PROFILE(f,thread,position)=587050*pow((8932/V_cbeta),1.2);

else if(t<=0.01375 && t>0.00854256)
   F_PROFILE(f,thread,position)=2144045.39-45844.5044*418.879*t;
else if(t<=0.0141667 && t>0.01375)
   F_PROFILE(f,thread,position)=44548090.82-3103131857*t;
else if(t<=0.015 && t>0.0146667)
   F_PROFILE(f,thread,position)=587050;
}
end_f_loop(f,thread)
}
DEFINE_PROFILE(unsteady_temperature,thread,position)
{
face_t f;
begin_f_loop(f,thread)
{
real t=RP_Get_Real("flow-time";
t=t-0.015*floor(t/0.015);
if(t<0.002083 && t>0)
   F_PROFILE(f,thread,position)=308.15;
else if(t<=0.00854256 && t>=0.002083)
  
   F_PROFILE(f,thread,position)=308.15*pow(587050*pow((8932/V_cbeta),1.2)/587050,0.16667);

else if(t<=0.01375 && t>0.00854256)
   F_PROFILE(f,thread,position)=308.15*pow((2144045.39-45844.5044*418.879*t)/587050,0.16667);
else if(t<=0.0141667 && t>0.01375)
   F_PROFILE(f,thread,position)=308.15*pow((44548090.82-3103131857*t)/587050,0.16667);
else if(t<=0.015 && t>0.0141667)
   F_PROFILE(f,thread,position)=308.15;
}
end_f_loop(f,thread)
}
回复此楼

» 猜你喜欢

前进
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

智能机器人

Robot (super robot)

我们都爱小木虫

找到一些相关的精华帖子,希望有用哦~

科研从小木虫开始,人人为我,我为人人
相关版块跳转 我要订阅楼主 pjhust 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 求调剂 +10 何气正 2026-04-13 11/550 2026-04-13 19:42 by lbsjt
[考研] 302求调剂 +10 易!? 2026-04-13 10/500 2026-04-13 19:04 by lbsjt
[考研] 复试调剂 +16 积极向上; 2026-04-10 18/900 2026-04-13 16:48 by fenglj492
[考研] 271求调剂 +30 2261744733 2026-04-11 34/1700 2026-04-13 15:24 by Kilig0317
[考研] 一志愿厦大生物学332求调剂 +11 池池池池池池 2026-04-08 11/550 2026-04-13 14:10 by 科研论
[考研] 求调剂 +11 月@163.com 2026-04-07 13/650 2026-04-11 22:55 by BruceLiu320
[考研] 求调剂 +10 璃茉一定上岸 2026-04-10 10/500 2026-04-11 13:31 by 1005715100
[考研] 求调剂 +6 archer.. 2026-04-09 8/400 2026-04-11 10:55 by zhq0425
[考研] 调剂 +12 卷卷卷心菜_ 2026-04-09 13/650 2026-04-10 22:36 by Ftglcn90
[考研] 263能源动力专硕求调剂 +3 加大号饭盒袋 2026-04-10 3/150 2026-04-10 22:23 by 286640313
[考研] 22408 366分,本科211,一志愿西工大 +4 Rubt 2026-04-09 4/200 2026-04-10 19:51 by chemisry
[考研] 0858求调剂 5+5 Gky09300550, 2026-04-10 8/400 2026-04-10 19:13 by chemisry
[考研] 本9 一志愿西工大085601 324求调剂 +5 wysyjs25 2026-04-10 5/250 2026-04-10 16:57 by luoyongfeng
[考研] 301求调剂 +5 149. 2026-04-10 5/250 2026-04-10 15:45 by 柴小白
[考研] 085601初试330分找调剂 +10 流心奶黄包l 2026-04-09 10/500 2026-04-10 08:14 by Sammy2
[考研] 本科211 工科085400 280分求调剂 可跨专业 +3 LZH(等待调剂中 2026-04-09 3/150 2026-04-09 21:29 by wutongshun
[考研] 337求调剂 +4 Gky09300550, 2026-04-09 4/200 2026-04-09 17:18 by 帕尔马拉特
[考研] 085801 总分275 本科新能源 求调剂 +8 bradoner 2026-04-08 9/450 2026-04-09 13:43 by only周
[考研] 287求调剂 +6 Fnhc 2026-04-07 6/300 2026-04-08 10:05 by xingguangj
[考研] 313求调剂 +3 十六拾陆 2026-04-07 3/150 2026-04-07 23:20 by lbsjt
信息提示
请填处理意见