24小时热门版块排行榜    

Znn3bq.jpeg
汕头大学海洋科学接受调剂
查看: 1330  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料工程281还有调剂机会吗 +32 xaw. 2026-04-11 33/1650 2026-04-13 15:21 by 张zhihao
[考研] 一志愿085502,267分求调剂 +17 再忙也要吃饭啊 2026-04-08 18/900 2026-04-13 09:57 by darren_yym
[考研] 一志愿0807 数一英一 313 有没有二轮调剂 +12 emokidd 2026-04-08 13/650 2026-04-13 08:32 by lhj2009
[考研] 366求调剂 +9 不知名的小卅 2026-04-11 9/450 2026-04-13 01:19 by 幸免 ..
[考研] 本科南方医科大学 一志愿985 药学学硕284分 求调剂 +5 弱水听文 2026-04-09 5/250 2026-04-12 13:16 by shengxi123
[考研] 296求调剂 +14 汪!?! 2026-04-10 16/800 2026-04-12 10:48 by zhouyuwinner
[考研] 求助调剂,跨调 +12 X十甫寸Y 2026-04-11 12/600 2026-04-12 10:20 by zhouyuwinner
[考研] 312求调剂 +6 李鸿飞飞 2026-04-06 6/300 2026-04-12 00:34 by 蓝云思雨
[考研] 电子信息279求调剂,有书读就行 +8 wwwooden 2026-04-08 11/550 2026-04-11 20:22 by cq2548
[考研] 359求调剂 +5 胃痉挛累了 2026-04-11 5/250 2026-04-11 19:55 by lbsjt
[考研] 本人女孩 +7 吼吼, 2026-04-10 9/450 2026-04-11 14:45 by ACS Nano——
[考研] 农业管理302分求调剂 +3 xuening1 2026-04-10 3/150 2026-04-11 10:18 by zhq0425
[考研] 本科211 工科085400 280分求调剂 可跨专业 +11 LZH(等待调剂中 2026-04-10 11/550 2026-04-11 08:39 by zhq0425
[考研] 263能源动力专硕求调剂 +4 加大号饭盒袋 2026-04-10 4/200 2026-04-10 20:52 by gong120082
[考研] 一志愿矿大,材料工程专硕314分,0856可调都可以 +15 无懈可击的巨人 2026-04-09 15/750 2026-04-10 18:10 by hmn_wj
[考研] 一志愿211,化学学硕,310分,本科重点双非,求调剂 +27 努力奋斗112 2026-04-07 30/1500 2026-04-10 15:06 by Kilig0317
[考研] 调剂 +19 2261744733 2026-04-08 19/950 2026-04-09 19:11 by vgtyfty
[考研] 085801 总分275 本科新能源 求调剂 +8 bradoner 2026-04-08 9/450 2026-04-09 13:43 by only周
[考研] 机械专硕273请求调剂 +6 庚申壬申 2026-04-07 6/300 2026-04-08 22:41 by bljnqdcc
[考研] 调剂求助(生物与医药) +6 @6952 2026-04-06 6/300 2026-04-07 23:52 by lys0704
信息提示
请填处理意见