24小时热门版块排行榜    

查看: 403  |  回复: 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-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 X5bHxOULxFqb 2026-10-30 14/700 2026-11-01 11:33 by Equinoxhua
[公派出国] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 1uaZt2kAZj2S 2026-10-31 5/250 2026-11-01 11:18 by Equinoxhua
[基金申请] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 hKe7U1kad5Sg 2026-10-31 6/300 2026-11-01 11:09 by BCdvBAEpt3zm
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 xvJS1BYDxvBS 2026-10-31 4/200 2026-11-01 11:01 by BCdvBAEpt3zm
[论文投稿] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 xvJS1BYDxvBS 2026-10-31 6/300 2026-11-01 10:53 by BCdvBAEpt3zm
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 xvJS1BYDxvBS 2026-10-31 9/450 2026-11-01 10:53 by BCdvBAEpt3zm
[基金申请] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 1uaZt2kAZj2S 2026-10-30 14/700 2026-11-01 10:43 by BCdvBAEpt3zm
[论文投稿] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 1uaZt2kAZj2S 2026-10-30 11/550 2026-11-01 10:35 by BCdvBAEpt3zm
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 1uaZt2kAZj2S 2026-10-30 13/650 2026-11-01 10:27 by BCdvBAEpt3zm
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 1uaZt2kAZj2S 2026-10-30 16/800 2026-11-01 07:01 by YXM1FDCSE8Bn
[基金申请] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 xvJS1BYDxvBS 2026-10-31 8/400 2026-11-01 05:48 by YXM1FDCSE8Bn
[博后之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 1uaZt2kAZj2S 2026-10-30 15/750 2026-11-01 05:22 by YXM1FDCSE8Bn
[考研] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 1uaZt2kAZj2S 2026-10-31 7/350 2026-11-01 04:04 by YXM1FDCSE8Bn
[基金申请] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 1uaZt2kAZj2S 2026-10-31 3/150 2026-11-01 04:04 by YXM1FDCSE8Bn
[找工作] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 1uaZt2kAZj2S 2026-10-31 6/300 2026-11-01 03:56 by YXM1FDCSE8Bn
[基金申请] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 1uaZt2kAZj2S 2026-10-30 10/500 2026-11-01 03:51 by YXM1FDCSE8Bn
[考博] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 1uaZt2kAZj2S 2026-10-31 5/250 2026-11-01 02:25 by YXM1FDCSE8Bn
[考研] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 1uaZt2kAZj2S 2026-10-30 12/600 2026-11-01 02:19 by YXM1FDCSE8Bn
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 X5bHxOULxFqb 2026-10-30 12/600 2026-11-01 00:15 by YXM1FDCSE8Bn
[考博] 申博求助 20+3 Katrinae 2026-10-26 4/200 2026-10-27 09:56 by 北京莱茵编辑
信息提示
请填处理意见