[noa]

I have a 2D channel that its edge is circular. I am moving the walls in y direction and this working. However, when I am not able to move the circular edge. This is the UDF I am using...


#include "udf.h"
#include "mem.h"
#include "dynamesh_tools.h"
/* Constants */
#define b 1.5e-3
#define T 20
#define a 0.0e-3
#define a1 10e-3
#define b1 0.5e-3
#define R 0.5e-3
#define phi 3.1415926535897932384626433832795
#define L 50
#define EdgeHeight 1.0e-3
/************************************************** ***/

DEFINE_GEOM(outlet,domain,dt,position)
{
Thread *tf = DT_THREAD (dt);
face_t f;
Node *node_p;
real x, y;
int n;

SET_DEFORMING_THREAD_FLAG (THREAD_T0 (tf));

begin_f_loop (f, tf)
{
f_node_loop (f, tf, n)
{
node_p = F_NODE (f, tf, n);
if (NODE_POS_NEED_UPDATE (node_p))
{
NODE_POS_UPDATED (node_p);
x = NODE_X (node_p);



y[0]=sqrt(fabs(pow(R,2)-pow((x-a1),2)))+(b1+CURRENT_TIME/T);

y[1]=-sqrt(fabs(pow(R,2)-pow((x-a1),2)))+(b1+CURRENT_TIME/T);

NODE_Y (node_p) = y;
}
}
}
end_f_loop (f, tf);
}