Fluent-UDF造波源程序集(共7页)

1、精选优质文档-倾情为你奉上法一：边界造波法程序一：inlet.c#include "udf.h" /* #include<stdio.h> (输入输出)*/ #include<stdlib.h>#include<math.h>#define HH 0.5 /*waver hight*/ /*不同波形需要修改的波形参数*/#define LL 0.95 /*波长*/#define g 9.81#define pi 3.#define d 0.5 /*water deepth */DEFINE_PROFILE(x_velocity,threa

2、d,index) real kk = 2.0*pi/LL; real c = sqrt(g*tanh(kk*d)/kk); /*参见一般性公式的推导*/ real TT = LL/c; real xND_ND; /* this will hold the position vector */ real y = 0; real s = 0; real ct = 0; /*相位角*/ face_t f; real t = CURRENT_TIME; real u = 0; t=RP_Get_Real("flow-time"); begin_f_loop(f,thread) /*

3、 loops over all faces in the thread passed in the DEFINE macro argument */ F_CENTROID(x,f,thread); y = x1; s = y+d; ct = kk*(x0-c*t); /*参见一般性公式的推导*/ if(y < 0.5*HH*sin(ct) /*水面以下，其中，0.5*HH*sin(ct)为波面方程*/ u = pi*HH*cosh(kk*s)*sin(ct)/(TT*sinh(kk*d); /*x方向速度分量公式，参见一般性公式的推导*/ else u = 0.0; /*水面以上流体单位

4、速度矢量的x方向分量*/ F_PROFILE(f,thread,index) = u; end_f_loop(f,thread)DEFINE_PROFILE(y_velocity,thread,index) real kk = 2.0*pi/LL; real c = sqrt(g*tanh(kk*d)/kk); real TT = LL/c; real xND_ND; /* this will hold the position vector */ real y = 0; real s = 0; real ct = 0; face_t f; real t = CURRENT_TIME; rea

5、l v = 0; t=RP_Get_Real("flow-time"); begin_f_loop(f,thread) F_CENTROID(x,f,thread); y = x1; s = y+d; ct = kk*(x0-c*t); if(y < 0.5*HH*sin(ct) v = pi*HH*sinh(kk*s)*cos(ct)/(TT*sinh(kk*d); else v = 0.0; F_PROFILE(f,thread,index) = v; end_f_loop(f,thread)DEFINE_PROFILE(voffactor,thread,inde

6、x)real kk = 2.0*pi/LL; real c = sqrt(g*tanh(kk*d)/kk);real TT = LL/c; real xND_ND; /* this will hold the position vector */ real y = 0; real s = 0; real ct = 0; face_t f; real t = CURRENT_TIME; t=RP_Get_Real("flow-time"); begin_f_loop(f,thread) F_CENTROID(x,f,thread); y = x1;s = y+d; ct =

7、kk*(x0-c*t); if(y < 0.5*HH*sin(ct) F_PROFILE(f,thread,index) = 1.0; else F_PROFILE(f,thread,index) = 0.0; end_f_loop(f,thread)程序二：wave.c （此程序同程序一大致相同）#include "udf.h" /* must be at the beginning of every UDF you write new9 case */real AA=0.5; /*waver amplitude*/real LL=0.95; /*不同波形需要修改的

8、波形参数*/real TT=0.78;real pi=3.real kk=2.0*pi/TT;real ww=2.0*pi/LL;real h=0.5; /*water deepth */real ux=1.0; /*此为何变量？*/DEFINE_PROFILE(x_velocity,thread,index) real xND_ND; /* this will hold the position vector */ real y; face_t f; real t = CURRENT_TIME; real u; begin_f_loop(f,thread) /* loops over all

9、 faces in the thread passed in the DEFINE macro argument */ F_CENTROID(x,f,thread); y = x1; if(y<(AA*cos(kk*x0-ww*t) /*(-a*sin(w*t)+2.0为入口的波面随时间的变化*/ /*u=a*w*(exp(k*(y+h)+exp(-k*(y+h)*-sin(w*t)/(exp(k*h)-exp(-k*h)*/ u=9.8*AA*kk/ww*cosh(kk*(y+h)*cos(kk*x0-ww*t)/cosh(kk*h); else u =0.0; /*水面以上流体单位速

10、度矢量的x方向分量*/ F_PROFILE(f,thread,index)=u; end_f_loop(f,thread)DEFINE_PROFILE(y_velocity,thread,index) real xND_ND; /* this will hold the position vector */ real y; face_t f; real t = CURRENT_TIME; real v; t=RP_Get_Real("flow-time"); begin_f_loop(f,thread) /* loops over all faces in the thre

11、ad passed in the DEFINE macro argument */ F_CENTROID(x,f,thread); y = x1; if(y<(AA*cos(kk*x0-ww*t) v=9.8*AA*kk/ww*sinh(kk*(y+h)*sin(kk*x0-ww*t)/cosh(kk*h); else v=0.0; F_PROFILE(f,thread,index)=v; end_f_loop(f,thread)DEFINE_PROFILE(voffactor,thread,index) real xND_ND; /* this will hold the positi

12、on vector */ real y; face_t f; real t = CURRENT_TIME; begin_f_loop(f,thread) /* loops over all faces in the thread passed in the DEFINE macro argument */ F_CENTROID(x,f,thread); y = x1; if(y<(AA*cos(kk*x0-ww*t) F_PROFILE(f,thread,index) = 1.0; else F_PROFILE(f,thread,index) = 0.0; end_f_loop(f,th

13、read)法二：推波板造波法（动网格）程序三：pushboard.c#include<stdio.h>#include"udf.h"#define T 5.8 /*周期*/#define S 1.04 /*冲程*/#define PI 3.DEFINE_CG_MOTION(pushboard, dt, cg_vel, cg_omega, time, dtime) real u = 0; real ww = 0; ww = 2*PI/T; if(time <= 2*T) u = ww*S*time*cos(ww*time)/(4*T); else u = w

14、w*S*cos(ww*time)/2; cg_vel0 = u;法三：摇板造波法程序四：11.c （三维造波）#include"udf.h"#define k 0.77#define w 2.75#define A 0.059 /*振幅*/#define l 1.6 /*波长*/DEFINE_CG_MOTION(yaoban,dt,vel,omega,time,dtime)vel0=0.0; /*x方向的速度*/vel1=0.0;vel2=0.0;omega0=0.0; omega1=0.0;omega2=atan(A/l)*w*cos(w*time); /*z方向的角频率

15、*/程序五：b.c#include "udf.h"static real G=9.81;static real VLA_M=100;static real VLA_F=0.0;static real VLA_V=0.0;DEFINE_CG_MOTION(vla,dt,vel,omega,time,dtime)Thread *t;Domain *d;real dv,CGND_ND,force2,moment2;if(!Data_Valid_P()return;NV_S(vel,=,0);NV_S(omega,=,0);d=THREAD_DOMAIN (DT_THREAD (Dynamic_Thread *)dt); t