[engrmansoor2534]

I have written a udf for 2D channel wall deformation. Now I want to apply the udf to the wall of a tube. the udf just displaces a certain line of nodes but not all the nodes on the tube wall.what changes should I apply?
the tube has axial direction in the x-axis.

[`e`]

What is your UDF and which different nodes are modified in the 3-D case? For example, if only the nodes in the X-Z plane are affected then you could add a conditional statement for z=0.

[engrmansoor2534]

Originally Posted by `e`

What is your UDF and which different nodes are modified in the 3-D case? For example, if only the nodes in the X-Z plane are affected then you could add a conditional statement for z=0.

Dear 'e'
Here is the udf for 2D channel case. which is deforming the walls according to the specified equation.

DEFINE_GRID_MOTION(grid_motion, domain, dt, time, dtime)
{
 face_t f;
 Thread *tf = DT_THREAD ((Dynamic_Thread *) dt);
 int n;
 Node *v;
 real x;
 real y;
 SET_DEFORMING_THREAD_FLAG(THREAD_T0(tf));

 /* Set flag to indicate that the current node's position
 has been updated, so that it will not be updated in a
 future pass through the loop. */

  begin_f_loop (f, tf)
   {
	f_node_loop (f, tf, n)
	 {
		v = F_NODE (f, tf, n);
		if (NODE_POS_NEED_UPDATE (v))
		  {
			NODE_POS_UPDATED(v);
				x = NODE_X (v);
				y = (D/2)+b*sin((2*pi/lam)*(x-c*CURRENT_TIME));
				NODE_Y (v)=y;

where D is tube diameter
b is wave amplitude and c is wave propagation speed
lam is wavelength

Now I want to apply this equation for a 3D TUBE case. the streamwise direction is same as the channel i.e. x-axis. TUBE diameter is along y and z axis. that is, nodes in both y and z directions need to be deformed in order to produce contraction and expansion in the tube.
what changes do you prefer to the above code so that it can be used to produce deformation in a tube.
I have tried the following changes i,e,
NODE_Y(v), NODE_Z(v)=y;
but this command deforms the the wall of the tube in a rectangular manner. the circular property of the tube is violated.
any suggestions would be highly appreciated.
thank you

[`e`]

The vertical distance, y, in the 2-D case (rectangle) is analogous to the radial distance, r, in the 3-D case (tube/pipe). One approach is to assume the nodes only move radially inwards/outwards from the origin and therefore have a constant angle with the X-Y plane.

For example:

r = (D/2.)+b*sin((2.*pi/lam)*(NODE_X(v)-c*CURRENT_TIME));
theta = atan(NODE_Z(v)/NODE_Y(v));
NODE_Y(v) = r * cos(theta);
NODE_Z(v) = r * sin(theta);