24小时热门版块排行榜    

Znn3bq.jpeg
查看: 359  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 化学专业申博 +3 赵子羊 2026-05-23 4/200 2026-05-24 18:10 by 工大学长
[硕博家园] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 1rx34o113h 2026-05-23 3/150 2026-05-24 17:41 by 0i3mu4vkjz
[基金申请] 评审有感 +16 popular289 2026-05-18 27/1350 2026-05-24 17:34 by hhs666
[教师之家] 论文撤稿了 +4 bjvtcliu 2026-05-24 7/350 2026-05-24 17:29 by bjvtcliu
[硕博家园] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 hvkbtfonbv 2026-05-23 4/200 2026-05-24 17:21 by 75ui6h7z2t
[考研] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 a2tycdlnq1 2026-05-23 5/250 2026-05-24 16:21 by hhx1yx9evi
[论文投稿] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 a2tycdlnq1 2026-05-23 4/200 2026-05-24 16:16 by hhx1yx9evi
[基金申请] 河北省自然科学基金 +6 Peterchao 2026-05-18 9/450 2026-05-24 16:02 by 130067131
[硕博家园] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 pmo95bazuy 2026-05-23 8/400 2026-05-24 15:56 by 1uy1ht2y9r
[基金申请] 西安交大新媒学院副院长用撤稿论文结题 +3 bjvtcliu 2026-05-24 5/250 2026-05-24 10:16 by kudofaye
[考博] 26/27申博自荐 10+4 ZXW0202 2026-05-22 9/450 2026-05-24 08:47 by bjvtcliu
[基金申请] 青B发送上会通知了吗 +5 chemBioBro 2026-05-22 7/350 2026-05-23 12:35 by zhuifengzhy
[论文投稿] 投稿求助,期刊 +4 希冀,有书读 2026-05-20 8/400 2026-05-22 10:16 by 希冀,有书读
[文学芳草园] 献血感触 +7 呀呀好傻 2026-05-19 13/650 2026-05-21 20:15 by 呀呀好傻
[基金申请] 国自然评分 +4 无名者登山 2026-05-20 5/250 2026-05-21 16:35 by swuq
[基金申请] 国自然上会要求 +7 无名者登山 2026-05-18 11/550 2026-05-21 15:50 by draco1987
[基金申请] 提交了我也来说说感想 +9 fummck 2026-05-20 10/500 2026-05-21 14:17 by draco1987
[有机交流] 反应很差,大量原料没有反应 5+3 Mr.Zot 2026-05-19 8/400 2026-05-20 22:19 by Equinoxhua
[考博] 如果工作了想读博,可以边工作边读全日制嘛? 30+3 铁达火车 2026-05-18 5/250 2026-05-20 09:33 by tfang
[考博] 博士申请 +5 星…… 2026-05-18 6/300 2026-05-18 23:49 by 糊糊涂涂好
信息提示
请填处理意见