naca动网格程序

#include "udf.h"#include "unsteady.h"#include "stdio.h" #include "stdlib.h"/************************************************************/real current_time = 0.0 ;Domain * domain ;Thread * thread ;real NV_VEC( origin ),NV_VEC( force ),NV_VEC( moment ) ;/************************************************************/DEFINE_CG_MOTION(AIRFOIL,dt,vel,omega,time,dtime){current_time = CURRENT_TIME ; vel[0] = 0.0 ; vel[1] = 0.0 ; vel[2] = 0.0 ;omega[0] = 0.0 ;omega[1] = 0.0 ;omega[2] = 5.4831*0.5* cos( 31.4159 * current_time ) ;Message("time=%f  omega=%f\n",current_time,omega[2]) ;}