# 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

 Thread Tools Display Modes Linear Mode

 Posting Rules You may not post new threads You may not post replies You may not post attachments You may not edit your posts BB code is On Smilies are On [IMG] code is On HTML code is OffTrackbacks are On Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post Dante CFX 13 May 31, 2012 08:47 sweetcherry FLUENT 1 May 15, 2008 01:48 Amida Oluwole Main CFD Forum 0 June 12, 2007 06:12 A Student CD-adapco 0 April 18, 2006 11:36 Arpan CD-adapco 1 March 9, 2006 08:52

All times are GMT -4. The time now is 20:41.