# who can help me about my udf on piston...

 Register Blogs Members List Search Today's Posts Mark Forums Read

 August 11, 2008, 05:25 who can help me about my udf on piston... #1 venus Guest   Posts: n/a i have to simulate a case about stirling engine, it is like the motion of piston,i write a udf to control the piston motion, but the result shows the pressure shakes and increases the pressure value.In normal situlation, the pressure should be shake on steady,the pressure can't always increase, who can help me to solve the problem, please!!! i had think is it my UDF wrong about the "force". can everyone help me!!!>" #include "udf.h" # define UDF_FILENAME "piston.dat" # define PI 3.141592654 # define rp 0.0067 # define l1 0.079 # define rpm 500. # define w 2.*PI*rpm/60. # define q w*time FILE *fp; DEFINE_CG_MOTION(piston1, dt, cg_vel, cg_omega, time, dtime) { Thread *t=DT_THREAD (dt); face_t f; real force, loc; real velo, pressure; real v1, v2, v3; fp= fopen("piston.dat","a"); NV_S(cg_vel,=,0.0); NV_S(cg_omega,=,0.0); if (!Data_Valid_P ()) return; force=0.0; begin_f_loop(f,t) { real *AA; AA = F_AREA_CACHE (f, t); force += F_P (f, t) * AA[0]; if (rp_axi) force *= 2.0 * PI; pressure = force / 0.0099*0.0099*PI; } end_f_loop (f, t) velo = rp*w*sin(PI/2.-q)+1./sqrt(2.*rp*cos(PI/2.-q)* 2.*rp*cos(PI/2.-q)-4.*(rp*rp-l1*l1))*(2.*rp*cos(PI/2.-q)) *rp*sin(PI/2.-q)*(w); loc = 0.01945+0.02+rp*cos(PI/2.-q)+0.5*sqrt ((2.*rp*cos(PI/2.-q))*(2.*rp*cos(PI/2.-q))-4.*(rp*rp-l1*l1)); v1=0.0099*0.0099*PI*0.07; v2=0.0095*0.0095*PI*0.045; v3=v1-v2-(loc*0.0099*0.0099*PI); cg_vel[0]=velo; Message ("\nUDF valve:%f,%f,%f,%f,%f,%f\n",time,loc,velo,dtime ,v3,pressure); fprintf (fp, "%e %e %e %e %e %e\n", time, loc, velo, dtime, v3,pressure); fclose (fp); } vahid-karaj likes this.

 May 28, 2014, 10:37 Question #2 New Member   lakhdar Join Date: May 2014 Posts: 1 Rep Power: 0 Hello I used fluently with Visual Studio in the study of Stirling Engine I have a problem when the engine speed decreases Increase heat spread Knowing that I used the sinusoidal motion