|
楼主 |
发表于 2008-3-6 15:45
|
显示全部楼层
我的UDF是这样的
让一端边界(造波板)振动产生波浪,
x=5sin(10*w)
/**********************************************************************/
/* unsteady.c */
/* UDF for 动边界动网格 */
/**********************************************************************/
#include <stdio.h>
#include "udf.h"
DEFINE_CG_MOTION(valve, velocity_inlet, thread, position, t, dtime)
{
real x[ND_ND];
real x;
face_t f;
begin_loop_f(f, thread)
{
F_CENTROID(x,f,thread);
x = x[0];
if(x==-100)
F_PROFILE(f, thread, position) = 5.0*sin(10.*t);
else
F_PROFILE(f, thread, position) = 0;
}
end_loop_f(f, thread)
} |
|