有人帮忙看一下这个UDF吗?

liy02

#include "udf.h"
#define PI 3.14159
#define GRAV 9.81
#define H2.0
#define D6.0
#define T6.0
#define L40.86

DEFINE_PROFILE(x_velocity, thread, position)
{
  real x[ND_ND];    /* this will hold the position vector                */
  real y;
  real AA,BB,CC,DD,EE,FF,GG,LL,ZZ,K;      /* y=position vector, AA,BB,etc are temporary stores */
  face_t f;
  begin_f_loop(f, thread)
    {
      real t = RP_Get_Real("flow-time");
  F_CENTROID(x,f,thread);
      y = x[1];
      ZZ= y - D;
  K = 2*PI/L;
  AA=  cosh (K*(ZZ + D));
  BB=  cosh (K* D);
  CC = sinh (K*D);
  DD=  cos  (PI/2.0 -2.0 *PI*(t/T));
  EE=  cos (2.0*(PI/2.0 -2.0 *PI*(t/T)));
  FF=  cosh (2.0*K*D);
  GG=  cosh (2.0*(K*(ZZ + D)));
  
  LL= D+(H*DD/2.0)+(PI*H*H/(8.0*L))*(BB/pow (CC,3.0))*(2.0+FF)*EE;
  
  if (y <= LL)        /* added ';='; sign after ';<'; to see change               */  
     F_PROFILE(f,thread,position)=H/2.0*(GRAV*T/L)*AA*DD/BB+0.75*PI*PI*H*H/(L*T)*GG*EE/pow(CC,4.0);
      
   else               /* line &#35;35  */
         F_PROFILE(f, thread, position) = 0;
  }
  end_f_loop(f, thread)
}

DEFINE_PROFILE(y_velocity, thread, position)
{
  real x[ND_ND];    /* this will hold the position vector                */
  real y;
  real K,MM,NN,OO,PP,QQ,RR,SS,TT,ZZZ;      /* y=position vector, AA,BB,etc are temporary stores */
  face_t f;
  begin_f_loop(f, thread)
    {
      real t = RP_Get_Real("flow-time");
  F_CENTROID(x,f,thread);
      y = x[1];
      ZZZ= y - D;
  K = 2*PI/L;
  
  MM=  sinh (K*(ZZZ + D));
  NN=  cosh (K* D);
  OO=  cos (PI/2.0 - (2.0 *PI* t/T));
  QQ=  sin (PI/2.0 - (2.0 *PI* t/T));
  RR=sinh (K*D);
  SS=cosh (2.0*K*D);
  TT=sinh (2.0*(K*(ZZZ + D)));
  
  PP= D+(H*OO/2.0)+(PI*H*H/(8.0*L))*(NN/pow (RR,3.0))*(2.0+SS)*OO;
  
  if (y <= PP)          /* added ';='; sign after ';<'; to see change            */
     F_PROFILE(f,thread,position)=H/2.0*(GRAV*T/L)*MM*QQ/NN+0.75*PI*PI*H*H/(L*T)*TT*QQ/pow(RR,4.0);
      
   else               
         F_PROFILE(f, thread, position) = 0;
  }
  end_f_loop(f, thread)
}