24小时热门版块排行榜    

查看: 402  |  回复: 0

夜独醉12377

新虫 (初入文坛)

[求助] fluent软件中UDF程序报错

#include "udf.h"
#define PI      3.14159
#define GRAV    9.81
#define D       5.0
#define H       1.0
#define T       4.0
#define L       23
#define FL      90
// x direction velocity  //
DEFINE_PROFILE(x_velocity, thread, position)
{
  real x[ND_ND];
  real y;
  real yu,K,W,Lu;
  face_t  f;
  begin_f_loop(f, thread)
    {
        Real t=RP_Get_Real("flow-time";
      F_CENTROID(x,f,thread);
      y=x[1];
      yu=Y-D;
        K=2.0*PI/L;
        W=2.0*PI/T;
      Lu=D+(H/2.0)*cos(K*0.0-W*t)+(PI*H*H/(8.0*L))*(cosh(K*D)/pow
(sinh(K*D),3.0))* (2.0+cosh(2.0*K*D))*cos(2.0*(K*0.0-W*t))};
      if (y<=Lu)
        F_PROFILE(f,thread,position)=(PI*H/T)*(cosh(K*(yu+D))/sinh(K*D))*cos(
K*0.0-W*t)+0.75*(PI*PI*H*H/(T*L))*(cosh(2.0*K*(yu+D))/pow(sinh(K*D),4.0))*co
s(2.0*(K*0.0-W*t));
          else
F_PROFILE(f, thread, position) =0;
  }
end_f_loop(f, thread)
}
// y direction velocity //
DEFINE_PROFILE(y_velocity, thread, position)
{
  real x[ND_ND];
  real y;
  real yu,K,W,PP;
  face_t f;
  begin_f_loop(f, thread)
    {
      real t =RP_Get_Real( "flow-time ";
      F_CENTROID(x,f,thread);
      y=x[1];
      yu=y-D;
        K=2.0*PI/L;
        W=2.0*PI/T;
      PP=D+(H/2.0)*cos (K*0.0-W*t)+(PI*H*H/(8.0*L))*(cosh(K*D)/pow(sinh(K*D),3.0))* (2.0+cosh(2.0*K*D))*cos(2.0*(K*0.0-W*t));
      if   (y<=PP)
        F_PROFILE(f,thread,position)=(PI*H/T)*(sinh(K*(yu+D))/sinh(K*D))*sin(
K*0.0-W*t)+0.75*(PI*PI*H*H/(T*L))*(sinh(2.0*K*(yu+D))/pow(sinh(K*D),4.0))*sin
(2.0*(K*0.0-W*t));
       else
        F_PROFILE(f, thread, position)=0;
      }
  end_f_loop(f, thread)
}
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

智能机器人

Robot (super robot)

我们都爱小木虫

相关版块跳转 我要订阅楼主 夜独醉12377 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 xvJS1BYDxvBS 2026-10-31 6/300 2026-10-31 23:23 by YXM1FDCSE8Bn
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 xvJS1BYDxvBS 2026-10-31 7/350 2026-10-31 23:18 by YXM1FDCSE8Bn
[考研] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 1uaZt2kAZj2S 2026-10-31 6/300 2026-10-31 23:17 by YXM1FDCSE8Bn
[博后之家] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 1uaZt2kAZj2S 2026-10-31 5/250 2026-10-31 23:10 by YXM1FDCSE8Bn
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 X5bHxOULxFqb 2026-10-30 14/700 2026-10-31 22:46 by YXM1FDCSE8Bn
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 X5bHxOULxFqb 2026-10-30 19/950 2026-10-31 22:44 by YXM1FDCSE8Bn
[论文投稿] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 X5bHxOULxFqb 2026-10-30 24/1200 2026-10-31 22:39 by YXM1FDCSE8Bn
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 X5bHxOULxFqb 2026-10-30 11/550 2026-10-31 22:39 by YXM1FDCSE8Bn
[考研] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 xvJS1BYDxvBS 2026-10-31 5/250 2026-10-31 21:59 by YXM1FDCSE8Bn
[找工作] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 1uaZt2kAZj2S 2026-10-31 4/200 2026-10-31 21:54 by YXM1FDCSE8Bn
[考博] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 1uaZt2kAZj2S 2026-10-31 4/200 2026-10-31 21:51 by YXM1FDCSE8Bn
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 1uaZt2kAZj2S 2026-10-31 4/200 2026-10-31 21:51 by YXM1FDCSE8Bn
[基金申请] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 1uaZt2kAZj2S 2026-10-30 9/450 2026-10-31 21:34 by YXM1FDCSE8Bn
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 1uaZt2kAZj2S 2026-10-30 11/550 2026-10-31 21:25 by YXM1FDCSE8Bn
[教师之家] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 X5bHxOULxFqb 2026-10-30 13/650 2026-10-31 17:12 by hKe7U1kad5Sg
[基金申请] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 X5bHxOULxFqb 2026-10-30 18/900 2026-10-31 17:05 by hKe7U1kad5Sg
[考博] 2027博士申请-临床检验诊断学、病原生物学、微生物方向 50+3 yyh@ 2026-10-28 4/200 2026-10-31 12:47 by 可淡不可忘
[硕博家园] 越来越担心硕士毕业后没有竞争力了 +8 zuesli 2026-10-27 14/700 2026-10-30 22:44 by 好运连连公主
[文学芳草园] 闲聊 +3 myrtle 2026-10-25 6/300 2026-10-29 15:28 by myrtle
[考博] 申博求助 20+3 Katrinae 2026-10-26 4/200 2026-10-27 09:56 by 北京莱茵编辑
信息提示
请填处理意见