CFD Online Discussion Forums

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)
-   FLUENT (https://www.cfd-online.com/Forums/fluent/)
-   -   Fatal Error signal in flapping wing analysis (https://www.cfd-online.com/Forums/fluent/146254-fatal-error-signal-flapping-wing-analysis.html)

gopalkrishnan December 22, 2014 09:18

Fatal Error signal in flapping wing analysis
 
Hi everyone

I m trying to 3D flapping wing analysis in fluent. I am using the following UDF.
I got segmentation fault while i preview mesh motion in dynamic mesh setting. My doubts are

1.What is boundary condition i have to give for wing wall.
2. what mesh zone settings?
3. what are dynamic mesh settings i have to assign?

The following procedure i followed

Read Mesh
Define UDF
Material and Boundary condition settings
Dynamic mesh settings

My UDF:

#include "udf.h"
#include "dynamesh_tools.h"
#define degree 180

real w, pi, A_Flap, A_Feat,Phi, dPhi, Theta, dTheta, Psi, dPsi,force,T,Freq,deltat,a;
real Amp,Feather;


DEFINE_CG_MOTION(wing, dt, vel, omega, time, dtime)
{

Thread *t;
face_t f;
real NV_VEC (A);
Freq=0.0;
Freq=20; //dtime is the physical time step defined in the journal file.

pi=3.141592654;

/* define motion variable */

Amp= 50; //This amp variable will be defined in our journal file.
Feather= 60;
A_Flap = Amp*pi/180;
A_Feat = Feather*pi/180;

w=2*pi*Freq; //omega (radians)
T=1.0/Freq;

Phi = -A_Flap*sin(w*time);
dPhi = -A_Flap*w*cos(w*time);
dTheta = -A_Feat*w*sin(w*time)*exp(sin(2*w*time-pi/2)-1)/1.1963;

omega[0] = dTheta*cos(Phi);
omega[1] = dTheta*sin(Phi);
omega[2] = dPhi;

t = DT_THREAD ((Dynamic_Thread *)dt);
force=0.0;
begin_f_loop (f,t)
{
F_AREA (A, f, t);
force += F_P (f, t) * A[1];
}
end_f_loop(f, t)
}

My dynamic Mesh settings i attached here. I given the Time step size 0.001 and Number of Time step 500. Please people help me on this


All times are GMT -4. The time now is 00:47.