[tadevilas]

For analysis of synthetic jet. to give the motion to vibrating membrane(diaphragm,diameter:50mm) for 2d geometry i used these logarithmic c program. now i wants to do same for 3g geometry.
Is this program work for it. Please give me suggestion if any modification is required.
C program:
#include "udf.h"
DEFINE_GRID_MOTION(logarithmic, domain, dt, time, dtime)
{ Thread *tf=DT_THREAD(dt);
face_t f;
Node *v;
real NV_VEC(velocity),NV_VEC(axis);
real NV_VEC(origin),NV_VEC(rvec);
real loco;
real freq;
real amp;
real diam;
int n;
SET_DEFORMING_THREAD_FLAG(THREAD_T0(tf));
freq=200;
amp=0.0002;
diam=0.050;
loco=2*3.141592654*amp*freq*cos(2*3.141592654*freq *time);
Message("time=%f,omeg=%f\n",time,loco);
NV_S(velocity,=,0.0);
NV_D(axis,=,0.0,1.0,0.0);
NV_D(origin,=,0.0,0.0,0.0);
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);
velocity[0]=loco*(1-(4*NODE_Y(v)*NODE_Y(v))/(diam*diam)+((8*NODE_Y(v)*NODE_Y(v))/(diam*diam))*(log((2*NODE_Y(v)+0.000000000001)/diam)));
NV_V_VS(rvec,=,NODE_COORD(v),+,velocity,*,dtime);
NV_V(NODE_COORD(v),=,rvec); } } }
end_f_loop (f, tf); }

[pakk]

If you want to know if the program works, then try it...


Why ask us?

[tadevilas]

sir,
same program not work for 3d.
so what should changes required for 3d

[pakk]

Can you be more specific than "not work"?

- Does it compile? If not, what happens instead?
- Does it give warnings? If so, when and which?
- Does it give errors? If so, when and which?
- Does it run?
- Does it crash? If so, when does it crash?
- Does it give wrong results? If so, which? And which results did you expect?

If you describe your problem in more details, you will get an answer sooner.

[tadevilas]

sir
it give error after pressing RUN
divergence detected in AMG solver :x momentum
and in second try its show :negative cell volume detected

[pakk]

Is the orientation in your 3D-model the same as in your 2D-model?


As in: should the membrane vibrate into the X-direction in both models?

[tadevilas]

sir
please refer the link.

[pakk]

No, but thanks.
I am sorry, I was under the false impression that you wanted to get help in solving your problem. I would have been happy to do that.

But I now understand that you want somebody else to solve your problem. That is not what I want to do. I am deeply sorry for the confusion.

[tadevilas]

thanks sir

[tadevilas]

I change the program for 3d.
But it's show the error msg that is udf library trying to load is not compiled for current platform.
I used visual Studio 2017 for as compiler. And some time it's give error fatal error C1189.
This is code. Please tell me how to avoid this error..


#include "udf.h"
DEFINE_GRID_MOTION(logarithmic, domain, dt, time, dtime)
{ Thread *tf=DT_THREAD(dt);
face_t f;
Node *v;
real NV_VEC(velocity),NV_VEC(axis);
real NV_VEC(origin),NV_VEC(rvec);
real loco;
real freq;
real amp;
real diam;
int n;
SET_DEFORMING_THREAD_FLAG(THREAD_T0(tf));
freq=200;
amp=0.00152;
diam=0.050;
loco=2*3.141592654*amp*freq*cos(2*3.141592654*freq *time);
Message("time=%f,omeg=%f\n",time,loco);
NV_S(velocity,=,0.0);
NV_D(axis,=,1.0,1.0,1.0);
NV_D(origin,=,0.0,0.0,0.0);
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);
velocity[2]=loco*(1-(4*(sqrt((NODE_X(v)*NODE_X(v))+(NODE_Y(v)*NODE_Y(v ))))*(sqrt((NODE_X(v)*NODE_X(v))+(NODE_Y(v)*NODE_Y (v)))))/(diam*diam)+((8*(sqrt((NODE_X(v)*NODE_X(v))+(NODE_ Y(v)*NODE_Y(v))))*(sqrt((NODE_X(v)*NODE_X(v))+(NOD E_Y(v)*NODE_Y(v)))))/(diam*diam))*(log((2*(sqrt((NODE_X(v)*NODE_X(v))+( NODE_Y(v)*NODE_Y(v))))+0.000000000001)/diam)));
NV_V_VS(rvec,=,NODE_COORD(v),+,velocity,*,dtime);
NV_V(NODE_COORD(v),=,rvec); } } }
end_f_loop (f, tf); }

[tadevilas]

Hello sir.
I'm trying to do simulation with serial processor and it take to much time for this. So now I'm trying with parallel process (6 core) is there code previous upload code is same apply or any changes required for it to run with parallel processer

[pakk]

According to you, the previous code gave fatal error C1189 in serial mode.

It would be very strange if it would suddenly work in parallel mode.


So: please talk to the tadevilas from October 2018.

[tadevilas]

Hello sir..
Now the code are run smoothly. I resolve this problem hence now I want to do it by parallel mod

[pakk]

Then... Try it!


If you want help, then give information. What you tried, what went wrong, and so on.


We can not look on your screen!

[Spicy]

I searched in the learning forum of Ansys and tried to find the course on 2d analysis using UDF, but I did not find a single piece of information related to it [specifically explaining everything related to the moving of the diaphragm using UDF]. Can you please help me with this, my group project is related to the 2d analysis using UDF [Synthetic Jet]. I even searched it on youtube but did not find information in which I would know how to move the diaphragm. this analysis is related to the first case that is quiescent flow.
Your help would be greatly appreciated.

[Spicy]

Originally Posted by tadevilas

I change the program for 3d.
But it's show the error msg that is udf library trying to load is not compiled for current platform.
I used visual Studio 2017 for as compiler. And some time it's give error fatal error C1189.
This is code. Please tell me how to avoid this error..


#include "udf.h"
DEFINE_GRID_MOTION(logarithmic, domain, dt, time, dtime)
{ Thread *tf=DT_THREAD(dt);
face_t f;
Node *v;
real NV_VEC(velocity),NV_VEC(axis);
real NV_VEC(origin),NV_VEC(rvec);
real loco;
real freq;
real amp;
real diam;
int n;
SET_DEFORMING_THREAD_FLAG(THREAD_T0(tf));
freq=200;
amp=0.00152;
diam=0.050;
loco=2*3.141592654*amp*freq*cos(2*3.141592654*freq *time);
Message("time=%f,omeg=%f\n",time,loco);
NV_S(velocity,=,0.0);
NV_D(axis,=,1.0,1.0,1.0);
NV_D(origin,=,0.0,0.0,0.0);
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);
velocity[2]=loco*(1-(4*(sqrt((NODE_X(v)*NODE_X(v))+(NODE_Y(v)*NODE_Y(v ))))*(sqrt((NODE_X(v)*NODE_X(v))+(NODE_Y(v)*NODE_Y (v)))))/(diam*diam)+((8*(sqrt((NODE_X(v)*NODE_X(v))+(NODE_ Y(v)*NODE_Y(v))))*(sqrt((NODE_X(v)*NODE_X(v))+(NOD E_Y(v)*NODE_Y(v)))))/(diam*diam))*(log((2*(sqrt((NODE_X(v)*NODE_X(v))+( NODE_Y(v)*NODE_Y(v))))+0.000000000001)/diam)));
NV_V_VS(rvec,=,NODE_COORD(v),+,velocity,*,dtime);
NV_V(NODE_COORD(v),=,rvec); } } }
end_f_loop (f, tf); }

i am using same code for 2d analysis in visual studio but it is showing many error like-
1] cannot find open source fie.
2] identifier logarithmic undefined; identifier domain undefined; identifier dt undefined etc
3]expected an expression
and there some more ... can help me solve this issue as I am not good at coding